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Abstract 

A boil-off flow calorimeter designed for use with 
corrosive gaS mixtures at temperatures up to 400 °C and 
pressures up to 20 MPa is described. Results of tests with 
N, and CO,, up to 10 MPa and 200 °C, indicate a continuing 
problem with heat leaks, despite changes in the construction 
and Ssadaexorn of the calorimeter in orddnts reduce these 
heat leaks. Some test runs were made for a gas mixture of 
CO., H., and CH,, with the results comparing reasonably well 
with predictions based on the Peng-Robinson equation of 
State. 

A survey of heat leak correction calculations made by 
other workers suggests that a heat leak correction 
correlation dependent on reciprocal process fluid flow rate 
may be appropriate if further improvements are realized as a 


result of the recommendations given here. 
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is tintroduction 

Calorimetry involves the measurement of energy property 
changes or energy interactions which occur in some process 
of interest. While some of these thermophysical properties 
can be predicted accurately for an increaSing number of 
Substances and their mixtures, their accurate measurement 
remains procteae: either for successful process design and 
operation, or for assessment and improvement of the quality 
of the available predictive methods. Today's need for 
enthalpy data, as well as other data, on coal gasification 
plant streams, synthetic natural gas streams, and processes 
involving steam and hydrocarbon mixtures at the associated 
higher temperatures and pressures (Yen et al., 1977; 
Wormald, 1982), demands that facilities be built for the 
purpose of determining these properties. 

The work described here was initiated in order to meet 
some of these perceived needs, and consequently a 
calorimeter was designed and constructed for use with 
corrosive mixtures at high temperatures and pressures. A 
Survey of previous work, a discussion of the relevant 
thermodynamics and the various modes of fluid calorimetry 
lead to a consideration of the important design criteria for 
comparative boil-off calorimeters. A description of the 
apparatus as currently operating is given. What is not 
contained here is an extensive review of the prior 
development work that was invested into the achievement of 


an operational device. The design of the calorimeter is 
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Similar in some respects to that of Sahgal et al. (1965). 
Freon-ii was used as the reference fluid in the device 
described here. Nitrogen and carbon dioxide, two gases 
whose thermal properties are well known, were used as 
Standards in order to calculate the heat leaks. A third 
fimrcd=.. composed: of) amma xtures of. CH, (61.7%), COs. (285.6%), 
and H, (9.7%), was also tested. When it became apparent 
that heat leaks were a major problem, efforts were made to 
modify the construction and operation of the apparatus in 
order to improve its performance. 

The method and precision of the calculations involved 
are summarised, and certain of the processed results are 
compared in an attempt to note trends, as well as to _ find 
probable causes of any variation in the results of similar 
runs. The heat leak correction methods of other workers 
were surveyed and their applicability to this work is 
examined. 

It is concluded that correlation of the measured 
property (in this case, enthalpy change from ambient), and 
hence heat leak, with the reciprocal flow rate would be 
appropriate if the present Peer cmance could be sufficiently 
improved. Recommendations which should result in some 


improvements are given. 
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2. Context 


2.1 Survey of the Literature 

A literature survey in the general area of calorimetry 
readily reveals the existence of previous surveys which, 
combined, form a rather extensive history and context for 
present work. 

These surveys are those of Partington and Shilling 
(1924), reviewing the calorimetric work done and reported 
prior to 1924; Masi (1954), who summarized "all of the 
literature reports of determinations of heat capacity of 
gases in the period 1924-1952, inclusive"; Faulkner (1959), 
who reviewed methods of constant volume calorimetry and 
isenthalpic expansion in addition to ae survey of flow 
calorimetric methods up to that time; Yesavage et al. 
(1967), who briefly discussed experimental methods, ren 
Surveyed much of the data for fluid mixtures with examples 
of how it may be analysed, and finally treated methods of 
prediction of certain thermal properties; and Mather (1978), 
who reviewed recent experimental and predictive methods for 
the determination of fluid enthalpies. 

Since the 1978 survey by Mather, the work of Barry et 
al. (1982), Christensen et al. GiO8-h) Clarke et: ali. 
(1979), Cunningham and Wilson (1978), Eubank et al. (1982), 
Juda et al. (1980), Kidnay and Yesavage and co-workers 


(McConnell (1976), Sharma (1977), Yesavage et al. (197-7): 
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Andrew (1978), Omid (1978) ), Lammers and co-workers 
(Lammers et al. (1978), van Kasteren and Zeldenrust 
(1979) ), Miyazaki et al. (1980), and Wormald (1979, 1982) 
and co-workers (Hutchings et al. (1978), Wormald et al. 
(1979), Richards and Wormald (1981), Richards et al. 
(1981) ) has been documented. 

Focusing the survey on comparative calorimetric 
methods, one finds at least two means of employing a 
reference fluid. The first involves the use of a reference 
liquid which undergoes a temperature change. Subsequent 
calculations require knowledge of the heat capacity of this 
liquid in order to determine the enthalpy change occurring 
in the process under study. 

The other main technique is the use of a vapourising 
reference fluid (at constant temperature). The external 
information required in this case is the data on enthalpy of 
vapourisation of the reference fluid. The first reported 
use of this technique appears to be the work of Nelson and 
Holcomb (1953), who used Freon-11 in the determination of 
enthalpies of hydrocarbon mixtures. Since then, as 
indicated in Table 2.1, many other workers have constructed 
equipment using the same principle in other designs and 
applications. (Some of the information in the table is 
extracted from Mather's survey (1978) and from McConnell's 


thesis (1976).) 
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2.2 Thermodynamic Considerations 

Because the enthalpy of a pure substance or mixture in 
a single homogeneous phase is a state function, changes in 
enthalpy are dependent only on the initial and final states 
and independent of the process or path along which the 
change was effected. Thus, the change in enthalpy is an 


exact differential and may be expressed as follows: 


n 
dh dh oh 
an = (38) ar + (23) dps + Rice) ax, (SY 
oT Pex a yx 2, 7 Ete * 
For pure fluids or mixtures of constant composition, the 
thermal effects of composition changes, i.e. mixing, are not 
present and equation (2.1) may be more simply written: 
dh = (= at + (38 dP (22:2,) 
x 


Ply 


For convenience, the partial derivatives above are, of 
course, need respectively as Cp, the constant pressure 
heat capacity, and @, the isothermal expansion coefficient. 
Sufficient knowledge of these quantities for a certain 
Substance or mixture is necessary for the construction of a 
pressure - enthalpy diagram' and allows calculation of an 
(integral) change of enthalpy corresponding to a change of 


State specified by changes in temperature and/or pressure. 


‘Necessary and sufficient would be knowing either Cp at one 
pressure and all temperatures, with @ at all temperatures 
and pressures, or @ at one temperature and all pressures, 
with Cp at all temperatures and pressures. 
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Since temperature and pressure are more eaSily observed 
in many situations than are other properties, equation (2.2) 
is a convenient form. However, the more fundamental 
relationship for the enthalpy property is the following, 
derived from the first and second laws of thermodynamics and 


the definition of enthalpy: 


dh = T ds + v dP (2.3) 


One can show that @, the isothermal expansion coefficient, 


is fundamentally expressed as follows: 


dh ds 
¢= (3) = (55 + (2.4) 
2 aP } ie 


Use of the fourth Maxwell relation, derived from the 


definition of the Gibbs free energy, 


-(35) : ty) (2.5) 
oP + oT P 
allows the expression of @ as a function of pressure, 


volume, and temperature only: 


dh ov 
oe 


Thus, it is seen that precise P-v-T data (physical property 
measurements) might be used as a check on this type of 


calorimetric data (thermal property measurements). 
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Precision in the P-v-T data is necessary since the 
derivative involved will tend to emphasize any errors which 
are present. 

The foregoing is in contrast to the case of the 
constant pressure heat capacity, Cp. It is not possible by 
Similar methods to express Cp as a function of P-v-T only 
and hence calorimetric methods for its determination are 
necesSary. 

In these methods, it is necessary to exchange energy 
with the substance in some way, while noting the resultant 
change in state of the substance by accurately observing 
temperatures and pressures. The first law of thermodynamics 
may be written as an energy balance: 


2 2 Z 


u u u 
[h+>+ ez]ém, - [h+5+ gz]om, + 6Q- SW= d([e +> + ez]m,_.) 


C7) 
For a steady state flow situation, changes in overall system 
mass and energy are zero. In flow calorimetry, no. shaft 
work is done on or by the system (any electrical energy 
being considered as a heat input) and potential and kinetic 
energy changes as the material moves through are made 
negligible by proper design and the use of small pressure 
drops, so that velocity changes are not large. Given all of 
this, equation (2.7) can be greatly Simplified to the form 


below: 
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where g is the amount of heat transfer per mole of process 
Fluids Further elaboration on this flow calorimetry 


equation appears in the following section. 
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2.3 Methodology 
Major Methods. As has been previously pointed out, 
there are published several extensive surveys of the general 
mfineld (of calorimetry which span at least the last 80 years 
of work in the area. What follows is a brief review of main 
flow calorimetry methods and their relative merits. 
Despite the comments of Masi (1954): 
"Obtaining and measuring a constant flow rate 
to the desired precision, measuring 
temperatures which accurately represent the 
state of the gas and reducing the heat leak of 
the calorimeter to manageable proportions are 
only three of the problems which are 
encountered, and which tend to make extremely 
complicated any calorimeter buat for 
precision work"; 
and the words of Barieau (1965), who decided against the 
direct measurement of enthalpies, concluding that 
"methods that require no mass flow rate or 
energy measurements are to be preferred [over 
direct methods, since]... it appears that the 
most Ci etieuit Quantity to determine in 
reducing data so that an enthalpy diagram can 
be constructed is the determination of the 
flow rate", 
it seems that the basic methods of flow calorimetry have 
Survived. The reasons for this include the directness with 
which Cp or other thermal properties can be calculated 
(since, for steady-state operation, the heat capacity of the 
apparatus itself does not enter the calculation), the wide 
range of temperatures and pressures for which the methods 
can be applied, and their applicability to mixtures. With 


modifications to increase its accuracy and applicability, 


flow calorimetry will continue as an important technique for 
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gathering reliable thermal property data. Under the general 
category of flow calorimetry, there are several main types 
used and discussed in the literature: 

1. differential flow micro-calorimetry 

2. "resistive heating" (direct) 

3. boil-off (comparative) 

4. matching boil-off 

The first of these is not subject to the criticisms of 
Masi, Barieau, and others about the precision of flow 
measurements, since it is actually a comparative method 
using the fluid at one set of conditions as its own 
reference for another set of conditions, all at the same 
flow rate. A ratio of heat capacities is obtained from the 
data and, therefore, a knowledge of ideal gas heat capacity, 
attainable from calorimetry or from spectroscopic 
determinations, is needed in addition to the experimental 
data to calculate actual values for Cp. 

The second type, "resistive heating" flow calorimetry, 
has the longest history of the remaining three. It involves 
the use of a carefully meaSured amount of electrical energy 
to restore or maintain certain flow conditions in the 
calorimeter. This energy is transferred to the fluid 
through dissipation from some resistive element either 
insulated from or in contact with the fluid. Thus the heat 


exchanged can be calculated as 


< 
tH 
tH 
Ww 


(2,89) 
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and, since q equals the measured heat transfer plus the 
transfer due to heat leaks, all per mole of process fluid, 


equation (2.8) can be rewritten: 


wt + ma | (2.10) 
Me i 

Note that this method requires accurate measurement of the 
process fluid flow rate. If this is done volumetrically, it 
is necessary to know the density of the fluid precisely in 
order to reduce the uncertainty in the resulting mass flow 
rate. This basic method has been used extensively, as 
Surveyed by Mather (1978) with more recent examples being 
Cunningham and Wilson (1978) and Clarke et al. GH982)i. to 
determine enthalpies and heat capacities of many fluids at a 
wide range of conditions and also for excess enthalpy (heat 
of mixing) experiments, as recently documented by Barry et 
al. (1982), Christensen et al. (1981), and Wormald et al. 
(1977) and many of his later co-workers. 

The third method, somewhat more recent than the second, 
is the comparative method of boil-off flow calorimetry, in 
which the fluid of interest undergoes a change of state and 
cools. The energy released is transferred to the reference 
Plu idira tee ives) Fibo ings pownt . Liquid nitrogen, Freon-11 
(trichlorofluoromethane), and water are the most commonly 
used reference fluids. The choice depends on operating 
conditions and also on the precision with which the enthalpy 


of vapourisation is known~ at the conditions used. 
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Generally, the fluid is caused to boil at ambient pressure. 


In this technique, the measured heat exchange is given by: 


dm 
sede Ld = = C214) 
Me. 


and the resulting operating equation becomes: 


A_m Q 
h aot a _L C2712) 
out in ° e 
ME Me 


Here it is seen that two flow rates must be measured. 
However, in the work that has been reported to date, the 
accuracies appear to be comparable to those of the resistive 
heating method. | 

McConnell (1976) has summarised the reasons which 
justify the choice of this method for certain applications. 
These include the ability to measure changes in enthalpy 
over larger ranges of conditions, the possibility of simpler 
construction, especially for higher temperature and pressure 
applications, and the adequacy of its accuracy for 
engineering design purposes. Another advantage may be that, 
rather than heating the process fluid as is done in the 
direct method, possibly causing its thermal decomposition, 
in the boil-off method the fluid is cooled, preventing 
further decomposition. 

The fourth method listed is a recent modification of 
the boil-off method which seeks to minimise possible 


problems of reference flow rate meaSurements, accuracy of 
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value of enthalpy of vapourisation, and heat leaks by 
attempting to duplicate the boil-off rate observed in 
"regular" operation by electrically heating the reference 
fluid under conditions which match the expected heat leaks. 
Eubank et al. (1982) recently presented details of this 
method. 

This method still requires knowledge of the process 
fluid flow rate, but the actual measurement of the reference 
fluid boil-off rate is replaced by a measurement of the 
amount of electrical energy required to produce an 
equivalent boil-off rate. Thus, precision and repeatability 
of reference fluid flow rates is required, but the need for 
accurate calibration is eliminated. The operating equation 


is therefore as follows: 
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Modes of Calorimeter Operation. A general 


pressure-enthalpy diagram for some mixture of constant 
composition is given in Figure (2.1). A number of process 
paths are also_ shown, indicating several modes of 
calorimetric investigation. 

Mode (a) is an expansion in which energy must be added 
to the fluid in order to maintain isothermal conditions. 
Any deviation from the temperature of interest must be 
corrected by the term on the far right of equation (2.14a) 


if the results are to be reported as isothermal. The values 
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of @ may be calculated from these data if the pressure drops 
involved are not too great. 

Mode (b) is the isobaric cooling or heating of the 
process fluid, possibly involving phase change. A small 
correction must usually be made for pressure drop since flow 
of the fluid cannot occur if inlet pressure exactly equals 
outlet pressure, but this adjustment may be negligible, 
especially where the isotherm is nearly vertical. The 
"resistive heating” method is generally used to increase the 
enthalpy of the incoming fluid, but could be used to restore 
the enthalpy of a fluid having cooled from its initial 
temperature. Boil-off methods involve only cooling of the 
process fluid. If sufficient data for small increments in 
temperature are taken, determination of the values of Cp is 
possible. 

Mode (c) is the Joule-Thomson experiment where a 
pressure drop under adiabatic flow conditions results in a 
temperature change. This change is often a decrease, but 


not necessarily, depending on the local slope of the 


isotherms. The Joule-Thomson coefficient, defined as 
follows: 
_ (a7 
— ( P) (2.15) 


may be determined from the data if sufficiently small 


pressure drops are used. 
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The thermodynamic consistency of the data gathered in 
the above three modes of experimentation may be examined 


uSing the relation 


6 te Ber (2.16) 


Te Ey, bi these properties have been measured at the states 
of interest. 

A fourth mode, not shown, might be some path in which 
changes in both pressure and temperature are significant. 
Such a process might be a simple expansion from some higher 
temperature and pressure to ambient conditions. The work on 
hydrocarbon enthalpies reported by Wiener (1966), who used a 
boil-off method, is an example of this. The data thus 
produced are enthalpy values based on a zero reference value 
at the ambient or near ambient conditions (eg., 25 °C and 
100 kPa) and might represent overall heating or cooling 
loads although they would be of little use in calculation of 


local values of Cp, @, and u. 
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3. Description of Equipment and Procedure 


3.1 Design Considerations 
The basic boil-off flow calorimeter equation was 


previously shown to be: 


(252) 


+ 
ae es e 


An examination of the terms of Rice equation provides the 
most direct indications of what the primary design 
considerations for this calorimeter should be. 

ies Heat Leaks. Minimisation of heat leaks by 
radiation, convection, and conduction is of great importance 
in the accuracy of the experimental results for enthalpy 
changes. As can be seen upon rearrangement of equation 
(2.12), the heat leak term is positive when an excess of 
reference fluid boils off: 
Qo. i‘ Mg Hout > hin? - dm (3.1) 
If the magnitude of the (positive) boil-off term ie greater 
than the magnitude of the (negative) “cooling” term, then 
heat is leaking into the boil-off chamber and Q is positive. 
Q is negative when too little reference fluid boils off. 

If the heat leak cannot be eliminated by modification 


of design and/or operation, it must be evaluated and 
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corrected, if possible, as a function of flow rate and/or 
temperature. This evaluation involves the use of a 
"Standard" process fluid and use of the known values of 
outlet and inlet enthalpies in equation (3.1). If such a 
correction correlation for heat leak is attainable, the 
calorimeter may be profitably used for enthalpy 
determinations, despite the existence of a measurable heat 
leak. 

2. Steady-state Operation and Capacity. The rapid 
achievement of steady-state operation requires that the 
"thermal inertia" or heat capacity of the calorimeter vessel 
be relatively low while the volume still contains enough of 
the boil-off (reference) fluid to handle anticipated heat 
flows. Thus the rapid "rise-time" may be seen in a 
trade-off relationship to "capacity". 

The mass of the reference fluid vessel and the 
exchanger (whether that be a simple coil of tubing or some 
more complex apparatus) will be of importance in keeping the 
time to reach steady-state to a minimum. The hot process 
fluid must heat at least the tubing walls to their 
steady-state temperature before,  conscvancy ene inlet 
temperature is achieved. Shortly thereafter, given a 
constant flow rate, the reference fluid boil-off rate will 
"level off". Despite the variation of flows that occur 
during experimental work, it 1S necessary to design around 
some figure for capacity, given by a residence time 


multiplied by an "average" flow rate of the reference fluid. 
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39 Adequate Flow Rates. Theoretically, the 


Significance of the heat leak term can be reduced by the use 
of an adequate flow rate of the process fluid. As well, the 
amount of heat lost by the fluid due to any of the heat 
transfer mechanisms should be lower per unit mass since the 
fluid has less total contact time with the inlet tubing. 
From a practical POTHEPOEM View; = CLE WSePalso: desirable tS 
measure a Significant ats eye flow rate and thus a larger 
process fluid flow rate is better. A constraint on this is, 
of course, that the higher the flow rate, the higher the 
pressure drop for a given restriction in the exchanger coil. 
This pressure drop effect must be considered if one wishes 
to report “isobaric” results. 

4, Accurate Pressure, Temperature, and Flow Rate 
MeaSurements. Pressures and temperatures must be accurately 
determined in order to accurately specify the state of the 
fluid at inlet and outlet. Pressure is perhaps the easiest 
to determine, uSing calibrated Heise gauges or pressure 
Gifferential cells with the associated dead weight 
equipment. The measurement of temperature in ae flow 
Situation in small tubes presents a larger possibility of 
error. One must decide either to measure temperature 
SEn-Tinewe Oro“ons tche-GSGexterior® of4ethe Wline If the 
thermocouple or platinum resistance thermometer being used 
is in contact with the tubing wall, there may be ae small 
conduction effect cooling the device below the bulk fluid 


temperature. Also, depending on the way the device is 
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inserted into the line, there will be some conduction down 
the probe itself, tending to heat or cool the operative part 
of the probe. For in-line measurement, there is some 
influence of flow rate due to convection. If the 
temperature is not taken in-line, there may also be 
radiation effects in addition to the above. 

Flow rate measurement errors probably comprise the 
second largest sources of inaccuracy after the heat leak 
problems, and, with the heat leaks minimised by proper 
design and/or correlated and corrected for, may actually 
represent the largest source of error in the value of 
enthalpy change for most calorimetric methods requiring 
values of flow rate. Thus, accurately calibrated flow 
meterS are necessary for both process fluid and reference 
fluids 

ae Variety of Process Fluids. The design should 
include consideration of a variety of fluids at a reasonably 
wide range of temperatures and pressures. As new processes 
become important in such areas as coal gasification and 
other SNG techniques, for example, increasingly more extreme 
process conditions of temperature, pressure, and corrosivity 
are being encountered. A calorimetric facility for handling 
these same process fluids must be able to withstand the same 
conditions using tough, corrosion resistant alloys in parts 
and tubing with which the process fluid comes into contact. 

6; Reference Fluid. Other design considerations 


relate to the handling and type of reference fluid. A fluid 
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with a boiling point near the ambient temperature should 
make possible significant reductions of heat leak in from 
the surroundings. The desired endpoint temperature for the 
process fluid enthalpy determinations will also be an 
important factor in deciding upon the best boiling 
temperature. Operation complexity will be reduced if an 
appropriate fluid with a normal boiling point near the 
temperature desired can be found. Increasing the boiling 
pressure above ambient will, of course, allow greater 
boiling temperatures for the same fluid, but this may bring 
further difficulties. The values of the enthalpy of 
vapourisation for the fluid should be well-known and 
documented for various conditions, since the boiling 
pressures will vary somewhat, even from day to day. A 
decision must be made as to how the flow rate will be 
measured, and whether the fluid will be vented or condensed 
and recirculated. The cost and availability of the fluid is 
important here, as well as its toxicity in the environment. 
7. Miscellaneous. Adequate provision must be made for 
preheating the process fluid and for flow and pressure 
control. Ease of repetition of the experiments will depend 
to some extent on the manner in which especially flow rates 


are controlled and adjusted. 
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3.2 Description of Apparatus 

The basic design of the calorimeter is similar to that 
described by Sahgal et al. (1965) used for low temperature 
work. Others have used modifications of this design for 
other work. (See the summarization of work by Sagara et 
al., Sood and Haselden, and Laverman and Selcukoglu in Table 
2.1.) One basic modification used in this work was the use 
of Freon-11 as the reference fluid and the consequent use of 
a higher range of temperatures than those used by Sahgal et 
al. Another was the recirculation of the F-11 rather than 
condensing and collecting the liquid for weighing, as other 
workers have generally done. 

The following description begins by tracing the path of 
the process fluid through the system, and will conclude with 
a description of the F-11 circulation system. A schematic 
illustration is given- in Figure 3.1 and a detailed layout of 
the calorimeter vessel and internals appears in Figure 3.2. 
Appendix 1 contains some of the design calculations that 
were done prior to construction and reconstruction of 
certain components of the system. 

From the supply cylinder, the process fluid was passed 
through a high pressure Tescom regulator and into the 
preheating section (not shown in Figure 3.1). In the case 
of CO7:; the supply cylinder was wrapped with 
"Heat-By-The-Yard" (HBTY) heating tape to raise the vapour 
pressure and so achieve higher inlet pressures for the CO, 


runs. Preheating of the process fluid was achieved simply 
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PD - pressure differential cooling fluid 
cell > out 
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MFM = mass flowmeter = 
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Figure 3.1 Schematic Diagram of Calorimeter 
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Table 3.1 : Legend for Figure 3.2 (Calorimeter Detail) 


A - Inner vessel vapour exit line 

B =)’Vacuum line 

C - Process *fintiduinier line 

D - Bolt holes (8) 

E - Inlet line radiation shield 

F - Vapourising Freon-11 

G - Approx. position of upper inlet T/C 
H 


- Approx. position of lower inlet T7¢ 


I - Inner radiation shield 

J - Guard radiation shield 

K = Approx. position of ourlecw © 

L - Heat exchanger core and coil 

M - Outer calorimeter vessel (evacuated) 
N - Guard vessel 

O - Inner vessel 

P - Screen vapour demister 

Q - Teflon wrapped aluminum ring seal 
R - Calorimeter lid (support plate) 

S - Process fluid outlet line 


T >= Liquid, Freon cecurne tine 
(glass float level indicator inserted here) 
U - Guard vessel vapour exit line 


Scale: O.D. of Outer Vessel - 5 in (12.7 cm) 


&) 
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¢ (Legend in Table 3.1) 


Figure 3.2 Calorimeter Internal Detail 
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by tightly wrapping the section of 1/8 in (3.18 mm) O.D. 
Stainless steel tubing? with HBTY heating tape. This 
section also included a Nupro sintered steel 60 micron 
filter element as a precaution against particulate build-up 
in downstream locations, and an in-line thermocouple. 

After preheating, the fluid flowed past the capped 
connection (not shown) for future attachment of a gas 
chromatograph (G.C.) sample line, the inlet pressure gauge 
connection, and through the top of the calorimeter into the 
evacuated outer calorimeter shell where the line size is 
reduced from 1/8 in (3.18 mm) to 1/16 in (1.59 mm). Two 
Thermoelectric 0.020 in (0.51 mm) dia. Type J thermocouples 
are inserted into the inlet line through a Conax ceramic 
multihole sealing gland which is teed into the line just 
after the preheating section. The junctions of these 
thermocouples are 2 in (51 mm) apart in order to provide an 
indication of the temperature gradient, if any, along the 
inlet line. The upper thermocouple was later used as an 
in-line resistance heater to raise the inlet temperature 
above that achievable by the preheater. An aluminum foil 
radiation shield was placed around the inlet line inside the 
calorimeter in an attempt to decrease any radiative heat 
transfer from the inlet line to the calorimeter internals. 

The fluid next passed into the exchanger composed of 


some 30 feet (approx. 9 m) of 0.018 in O.D. by 0.00425 in 


*Type 316 stainless steel was used for all process fluid 
tubing. 
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wall (0.46 mm O.D. by 0.11 mm wall) stainless steel tubing 
coiled around a 1.25 in (31.8 mm) copper core which is in 
thermal contact with the F-11 in the inner chamber of the 
calorimeter vessel. A radiation shield assembly was also 
constructed in order to minimise radiative heat losses from 
the exchanger coil to the outer vessel wall. 

Upon leaving the exchanger, the fluid had undergone a 
small pressure drop and was 2 to 3 °C above the boiling 
temperature of the F-11. This temperature was measured by a 
Op020arinnd (025 linmm) dediac Thermoelectric iron-Constantan 
thermocouple inserted into the line in a manner similar to 
that of the inlet line. The outlet line leaves’ the 
evacuated outer vessel of the calorimeter through the top 
plate. At this point, a second Heise gauge was connected 
for outlet pressure observation. A capped connection for 
future attachment of aG.C. sample line was also provided 
here. 

The fluid then passed through a Nupro fine metering 
valve and a Whitey shut-off valve to a knock-out drum from 
which any entrained condensed fluid may be collected*, and 
on tobayPrecisions Scientific: wetptestometer for volumetric 
measurement. A mercury thermometer located in the wet test 
meter indicated the meter temperature. The fluid leaving 


7It should be noted that passage through the needle valve 
into the knock-out drum is, in effect, an atmospheric flash. 
For runs in which a mixture undergoes some condensation, the 
amounts and compositions of the phases at this point will be 
different from the amounts and compositions of the phases 
leaving the calorimeter. The mixture used in this work, 
however, did not change phase in the calorimeter. 
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the wet test meter was vented to the air purge system of the 
laboratory. 

Freon-11 was chosen as the reference fluid Since its 
boiling point is near room temperature at the ambient 
pressure, which is near 93 kPa in this location. As such, 
the heat leak was reduced from that possible for fluids of 
lower boiling points. The vapourising F-11 was contained in 
the inner and guard boil-off chambers located within the 
calorimeter vessel. The guard vessel surrounds the inner 
chamber and was intended to be the sink for any heat leaking 
into the calorimeter from the surroundings. Sanceneliquid 
F-11 flowed into the inner chamber from the guard chamber 
and the pressure differential between the chambers was 
monitored and maintained at a small value, the levels of 
F-11 in each were about the same and the height of the glass 
float in the the guard vessel was thus an indication of the 
level in the inner vessel. This float assembly replaced the 
unsuccessful thermistor-based level sensing system which was 
initially installed. Level changes during a run can cause 
changes in relative heat leak and small changes in enthalpy 
of vapourisation and boiling temperature. Because the F-11 
was condensed and circulated back to the calorimeter, the 
liquid level did not change appreciably once steady 
operation was achieved, but the float was necessary to gauge 
the extent of filling. 

The F-11 vapourised in the inner cell passed through a 


simple demister and left the inner chamber, flowing by the 
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pressure differential cell, in-line thermocouple, and 
shut-off valve to a calibrated Brooks linear mass flow meter 
assembly and then to the reflux condenser through the vapour 
Space in the liquid reservoir. The vapour boiled off in the 
guard vessel normally passed by the pressure differential 
cell, through a needle valve and into the reflux condenser, 
also by way of rhe liquid reservoir. After some unsuccesful 
attempts to operate an in-line condenser, an all glass 
reflux condenser was installed to operate at ambient 
pressure and use cooling water at 2 to 3 °C. At this 
temperature the vapour pressure of F-11 is 45 kPa, and very 
little was lost to the surroundings. Approximately 100 in 
Wemoum) #Ote0 + 2ouineOLr Debye 0.0G25einewall (Orson) OFD. by 
1.59 mm wall) glass tubing was formed into a coil of about 
1.5 in (38 mm) in diameter and 20 in (0.5 m) high to provide 
the surface area for condensation. A calcium chloride 
dessicant was used to trap any water vapour that otherwise 
may have entered the top of the condenser. Vapour entered 
and liquid left the condenser through the connection to’ the 
reservoir immediately below the condenser. 

From the reservoir the liquid F-11 returned to the 
guard vessel of the. calorimeter through a _ liquid leg, 
shut-off valve, and past an in-iine thermocouple for 
observation of the return temperature. The liquid leg was 
meant to serve aS a vapour lock forcing all vapour flow 
through the lines described earlier and thus preventing the 


"short-circuiting" of guard vapour (and possibly some of the 
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3.3 Instrumentation 

Thermocouples. The thermocouples used were 
iron-Constantan (type J) grounded thermocouples. The three 
of these which were located in the fluid inlet and outlet 
lines inside the calorimeter were 0.020 in (0.51 mm) in 
diameter. The two thermocouples in the inlet line were 
manufactured by Omega, and the one in the outlet line by 
Thermoelectric. All others were 0.0625 in (1.59 mm) in 
diameter, also manufactured by Thermoelectric. Each 
thermocouple was connected to a separate ice-point reference 
junction and the resulting e.m.f£. was read through a 
Thermoelectric multipoint switch by a Hewlett-Packard Model 
3450A Digital Multi-Function Meter, which was also used to 
monitor signal voltages generated by other instruments 
described in the following sections. 

All thermocouples were calibrated in the Chemical 
Engineering Instrument Shop against a platinum resistance 
thermometer. The process fluid thermocouples were 
calibrated for the temperatures from room temperature 
(approx. 22 °C) to about 400 °C. For calculation purposes, 
theseacalibuationjdata were correlated by» a least» squares 


fit of equations of the following form: 
Y = A + B*M + C/M + D&ln(M) h3<2) 


where M is the thermocouple e.m.f. in millivolts and Y is 


one of the following functions of temperature: 
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19 WoT (TA a ne SC (32.248) 
aoe (Tas nueK) (3535) 
By)» Yesee ant) (Maine Tce (23c) 
4) Y = ln(T) CTAGNOER:) (2034) 


From a comparison of the above four function types for each 
thermocouple, type (1) generally appeared to be the most 
accurate, based on the sum of the sQuared residual errors in 
temperature. The type (1) correlation was used for all 
thermocouples except those manufactured by Omega, for which 
the type (3) correlations were more accurate. It was 
estimated that the resulting precision of the temperature 
observations was +0.05 °C. The accuracy of any temperature 
determined by these correlations in the range used in this 
work was +0.1 °C. 

The calibration charts, accompanied by the correlations 
and corresponding parameter values, appear in Appendix 2. 

Pressure Measurements. The calorimeter inlet and 
outlet pressures were observed using two Heise gauges with 
316 stainless steel Bourdon tubes. The specifications of 


each are given below: 


Inlet gauge Outlet gauge 
Full scale 3000 psig (20.7 MPa) 10 MPa (1450 psig) 
Accuracy better than 0.1% F.S. better than 0.1% F.S. 


Hysteresis minimal below 1500 psig minimal 


Precision 1 psig (6.9 kPa) 5 kPa (0.7 psig) 
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These gauges could be individually isolated from the flow 
line and compared to check mutual accuracies, if desired. 
If necessary, they could be vented rapidly to the laboratory 
air purge System via a shut-off valve in the event of 
accidental over-pressuring. The absolute pressures were 
calculated by adding these gauge pressures to the barometric 
pressure, which was read froma Fisher Scientific mercury 
barometer. Calibration details for the Heise gauges are 
provided in Appendix 2. 

The calibrations of the 3 in (76 mm) Matheson gauges 
used in the regulator assembly were not extensively checked 
Since these gauges provided approximate indications of the 
supply and upstream feed pressures. 

A Celesco Transducer Products differential pressure 
cell was situated between the guard vessel and inner’ vessel 
vapour lines. It was calibrated to provide a 0 to 10 volt 
Signal tcorrespondingiato am0rtomi0leus psig (CO toed szatkPa)) 
differential pressure range. This range was adequate for 
the operating conditions encountered. 

Flow Meters. After some initial false starts involving 
dry test meters for the measurement of the F-11 vapour flow, 
a Brooks linear mass flow meter (Model No. 5810 B1F2E5A) was 
installed and was successfully employed for measurements of 
low to moderate vapour flow rates. For moderate to high 
vapour rates, a Brooks Model No. 5810 BiJ2E5A flow meter was 
used. An attempt was made to shield the flow meter in use 


from room ‘air currents; «but © itis was> «not: insulated >= and 
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operated near room temperature. The actual vapour 
temperature wasS measured by an in-line thermocouple just 
prior to entry into the flow meter. 

These flow meters produce a full scale linear signal 
ranging from 0 to 5 volts corresponding to 0 to 500 standard 
cubic centimeters of N, per minute (SCCM N,) and 0 to 5000 
SCCM N; for Models 5810 BIF2E5A and 5810 BIJ2E5A, 
respectively. This signal was recorded on a Hewlett Packard 
Model No. 7127A Strip Chart Recorder and monitored through 
the Thermoelectric multipoint switch by the H.P. digital 
multimeter described earlier. Their calibration for F-11 
was accomplished by vapourising the F-11 at certain rates, 
passing this vapour through the meters, condensing the 
vapour, and weighing the condensed liquid collected ina 
certain period of time. Efforts were made to prevent vapour 
losses. The data thus gathered for each meter were 
correlated by a linear least squares fit. The correlations 
and the figures comparing the fit with the calibration data 
appear in Appendix 2. The calibration accuracies were 
better than 25 mg/min for the smaller capacity meter (flow 
eee A) and 100 mg/min for the larger capacity meter (flow 
meter B). For each, the meaSurement precision was +0.01 
Volt, or +0.2)%°0f full seale, ‘corresponding®=<to: 2) mg/min 
(0.015 mmol/min) for flow meter A and 16.9 mg/min (0.12 
mmol/min) for flow meter B. 

ABPrecision,Scientificm0..mcubicy foot. (2.83) 1)vwetatest 


meter uSing distilled water was used to measure the volume 
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of process fluid passing through the calorimeter during the 
period of the run. The meter could be read to the nearest 
0.0005 cubic foot (14.1 ml) and the accuracy was +0.5%. 
This meter operated at the room temperature, indicated by an 
internal mercury thermometer, and at the ambient barometric 
pressure. The process fluid mass flow rate during a run was 
calculated from this volume, temperature, and ambient 
pressure using appropriate values of the compressibility 
factor ganda ithereduratid onoy ofhe thest run: The meter was 
calibrated with a Precision Scientific 0.1 cubic foot (2.83 
1) prover at several different flow rates. 

Power Supplies. For all applications of the 
Heat-By-The-Yard heating tape, the 120 volt Variac power 
Supplies were used. A Heathkit Model IP-27 regulated power 
Supply was used to maintain a 300 mA current through the 
iron-Constantan in-line process fluid heater, which provided 


a resistance of approximately 35 ohms. 
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3.4 Experimental Procedure 

Three different fluids were used in a series of test 
runs for the approximate inlet pressures of 1400, 1000, and 
500 psig (9.7, 6.9, and 3.4 MPa). At each of these nominal 
test pressures, three or four flow rates were used, 
resulting in different inlet temperatures for a given 
preheater duty. The fluids were nitrogen, carbon dioxide, 
End aimixturevdotitethS .gXde CO.mb28:6%) + handh€HyodiGjn 7X « 

A general procedure for the test runs may be itemised 


as follows: 


1. The refrigeration unit for condenser cooling water was 


Started:and maintained at! ‘Isto 2 °C. 


2. The coolant circulation pump was started to cool the 


F-11 condenser. 


3. F-11 was charged to the calorimeter and the liquid 


reservoir through the chilled reflux condenser. 


4, The process fluid preheater power supply was turned on. 


=) The vacuum in the calorimeter outer vessel was checked 


and the vessel was re-evacuated if necessary. 


6. The Heise gauges were zeroed and the thermocouple 


reference ice-water dewar was refilled. 
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7. The supply cylinder valve was opened and the regulator 


was set to provide the desired inlet pressure. 


8. The power supply for the in-line process fluid heater in 


the calorimeter inlet line was activated. 


9. The F-11 vapour flow rate recorder was turned on, as was 
the digital multimeter used to monitor the thermocouples and 


the mass flow meter. 


70). The outlet needle valve was adjusted concurrently with 
the adjustment of the regulator to set the process fluid 
flow rate by setting the inlet pressure and the pressure 
drop across the exchanger. This was continued as necessary 
during startup to enhance achievement of steady operation, 
which was assumed when the inlet temperature, inlet 
pressure, outlet pressure, and F-11 vapour flow rate were 


constant. 


is The timer was started simultaneously with observation 
of the volume indicated on the wet test meter. Shortly 
thereafter, the F-11 vapour flow meter signal, and all 


thermocouple and pressure gauge readings were recorded. 


125 Readings were repeated and recorded every ten to 


fifteen minutes. 
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13. Some gradual pressure changes, if any, were corrected 


by minor adjustments of the regulator and the needle valve. 


14. After a sufficient period of steady operation was 
maintained and recorded, the process fluid flow rate was 
altered by adjustments of the outlet needle valve to cause a 
different ‘pressure drop. This occasionally required some 


adjustment of the regulator as well. 


LS: Steps 10 through 14 were repeated as often as required 


to generate the desired data, then step 16 was undertaken. 


or The process fluid heater power supplies, F-11 vapour 
flow rate recorder, and digital multimeter were shut off. 
The supply cylinder valve was closed. The refrigeration 
unit and the condenser coolant circulation pump remained 
running to prevent excessive losses of F-1i, especially if 


the series of runs was not yet completed. 


There were several deviations from this procedure for 
runs with CO.: | 

- Prior to step 7, the CO, supply cylinder was heated 
with HBTY heating tape to increase the supply pressure above 
the normal room temperature vapour pressure Ore cO-mtOlabouc 
ig00rpsig 011.7 MPa), 

- There was a large cooling effect at the outlet needle 


valve which tended to alter the original setting of the 
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pressure drop as the run continued. A short section of HBTY 
heating tape was wrapped around the valve and activated to 


warm it. This resulted in a much more steady operation. 
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4. Results and Discussion 


4.1 Calculation Procedure 

The experimental data collected during a series of runs 
were entered into data files similar to the sample given in 
Appendix 5.. The program CALOR, written specifically for the 
purpose, and given in Appendix 3, read these files and 
converted the data to the more useful units for temperature, 
pressure, and flow rate. With these values and a data file 
providing F-11 enthalpy of vapourisation values dependent on 
the Saturation pressure, the program calculated the 
experimental results for the process fluid enthalpy change. 

Further, for N, and CO,, which were used as standards, 
the program CALOR calculated predicted values of enthalpy 
change by using the inlet. “and outlet. conditions of 
temperature and pressure. This was done by taking the 
difference of the values predicted for each state uSing a 
two-dimensional linear interpolation of enthalpy data 
extracted from the I1.U.P.A.C. tables on N, (Angus et al., 
1979) and on CO, (Angus et al., 1973), and provided to -the 
program in appropriate data files of which copies appear in 
Appendix 4. In order to check the validity of linear 
interpolation for these data, the fitting and interpolation 
algorithms of Forsythe et al. (1977) were later used in a 
two-dimensional cubic spline interpolation in the same data 


files. The accuracy of these interpolations will be 
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discussed later. 

The program compared the experimental values with their 
associated predicted values of enthalpy change to calculate 
the heat leak at the particular experimental conditions, 
after the manner of equation (3.1). 

As previously mentioned, the third process fluid used 
in the calorimeter was a mixture of COpe en Hava ONG Chie NO 
Standard thermodynamic data are available for this mixture 
anGmmtnuSse ites es typicaler or sf lulds eutor which accurate 
calorimetric data may yet be needed. The experimental 
results in this case were compared to the predictions of the 
EQUI-PHASE program. This program 1s based on the 
"two-constant" Peng-Robinson equation of state (Peng and 
Robinson, 1976)) and is used extensively to predict the 
behaviour and properties of many mixtures of interest to the 
hydrocarbon processing industries. 

The experimental conditions and results, which comprise 


the output of the program CALOR, are listed in Appendix 6. 
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4.2 Precision Analysis 

Before the discussion of the accuracy of the 
calorimeter as a function of heat leaks, operational 
methods, and design decisions, it is appropriate to examine 
the effects of the precision of the experimental 
measurements on the certainty of the calculated results for 
enthalpy change. An Bxamination of the precision of the 
corresponding predicted values of enthalpy change is also 
necessary in order to determine ne limits of uncertainty 
for the estimation of heat leaks in the calorimeter. 

The total derivative of the measured heat exchange, or, 


more specifically, the derivative of equation (2.11) is: 


fez iiioc Uieittenss (4.1) 
: Aleda Zi sles meer da. Teg ao ities 2 sins ; 
Me Me (m-) 


The precision of the experimental value of enthalpy change 
for a given run may be calculated conservatively by 
approximation of the differentials in equation (4.1) by the 
uncertainties in the corresponding variables. Since some of 
these experimental variables were themselves determined from 
other measured values, calculations of their precision were 
done in a similar fashion. Since especially the reference 
and process fluid flow rates change from run to run, a fair 
estimate of the precision of the measured enthalpy change 
must consider several different sets of experimental 


conditions. 
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Calculations for a number of representative runs with 
both N, and CO,, are given in Appendix 7. The precision of 
the experimental values of enthalpy change for this work was 
thus determined to be better than +485 J/mol for N, at the 
lower flow rate (approx 10 mmol/min), changing to +140 J/mol 
as flow rate increased to 35 mmol/min, and +485 J/mol for 
CO, changing to as much as +575 J/mol when the higher 
capacity mass flow meter was used at the larger flow rates. 

The corresponding predicted value of the enthalpy 
change was calculated from a difference of outlet and inlet 
values Of enthalpy which were determined using a 
two-dimensional linear interpolation as previously 
indicated. This interpolation scheme is summarised in the 
equations below to clearly indicate how the calculation of 


precision proceeded. 


bh ib os = 
Dy p = * Li |: ¥ i) >A Doane (42) 
fe oe 
te | ge + hp- ot (4.3) 
: Po —aPe JT : 
ee tbr lje en 
Dy p = PI. - T ) + Arm p (4.4) 


The question of the accuracy of linear interpolation is 
momentarily left aside. Calculation of the precision was 
done by evaluating the total derivative of equation (4.4), 
and this involved similar evaluations of the derivatives of 


equations (4.2) and (4.3). If it can be assumed that the 
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values of enthalpy given in the I.U.P.A.C. tables are exact, 
or at least much more precise than _ the experimental 
measurements of temperature and pressure, then the precision 
of a predicted value of enthalpy depends upon the precisions 
of temperature and pressure in the following way, derived in 


Appendix 7: 
dh = (j+f£P)aT + (k+£T)AP (A7.4) 


where the constants are calculated from only the tabulated 
values. 

The parallel of equation (A7.4) with equation (2.2) is 
obvious and can, indeed, be used to judge whether the linear 
method of interpolation is appropriate for the data set. If 
Comparisons of Gp “with” “(j+tP) “"and*® O wren CkKrET) ae any 
condition (T,P) do not indicate near-equality, then the 
assumptions of linearity of enthalpy with temperature and 
pressure are not good in the vicinity of that state and the 
linear interpolation method is unsuitable for use at that 
point. A sample of these comparisons is given in Appendix 7 
and is representative of the buik of the points tested. The 
conclusion drawn was that two-dimensional linear 
interpolation was a sufficiently accurate method for the N, 
data set and much of the CO, data set. 

This conclusion was tested further by comparing the 
results of a more versatile interpolation method which does 


a better fit of any curvature present in the data. A run by 
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run comparison of results using a two-dimensional cubic 
spline interpolation with the results of the two-dimensional 
linear interpolation was done for both N, and CO,. For the 
N, runs, the enthalpy changes calculated from the _ spline 
interpolations were, on average, only 0.034 % larger than 
the changes calculated using linearly interpolated values of 
enthalpy. Similarly for CO, runs, the calculated enthalpy 
changes were only an average of 0.01490 larger when 
calculated using the spline interpolation. It seems clear, 
then, that the major contribution to ene uncertainty of the 
predicted enthalpy change was due to the precision of the 
experimental measurements of temperature and pressure, and 
not, at least in this work, to the method of representation 
of the tabular data. This uncertainty in the enthalpy 
change values was +7.0 J/mol for all N, runs and from +55 to 
+18 J/mol, decreasing as run pressure decreased, for CO,. 
Finally, the precision of the values of the heat leak 


for a certain run may be determined as follows: 


= d(Ah oa cp Lei eee (4.5) 


The precision of the calculated heat Teak values ranged from 
+490 J/mol to +150 J/mol for N,, as flow rate increased, and 
from +5170 J/mol to +630 J/mol for 'cO;, ‘depending on which 
flow meter was used. Since the magnitude of uncertainty in 
the measured enthalpy change is larger than that of the 


predicted value, it appears that the first attempts to 
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improve the precision must be in the experimental 
measurements, especially those of flow rate. The words of 


Masi (1954), quoted earlier, appear to be at least partially 


borne out. 


ye z 


lstneadveqes ‘9h? AE 


30 ebtow «af? 46987 volt ty a 


yiisizieq tenet ts ad°oF 


48 


4.3 Heat Leaks and Calorimeter Accuracy 

Examination of the results of the experiments, given in 
Appendix 6, and as illustrated in Figures 4.1, 4.2, and 4.3, 
reveals that the performance of the calorimeter was rather 
severely hampered by the presence of heat leaks which ranged 
from approximately 0 to +35 % of the measured boil-off heat 
exchange. In other words, the predicted enthalpy changes 
were anywhere from 35 % below to 35 % above the 
experimentally determined values. 
Effects of desi decisions The 


and change 


on 


discussion on precision in the preceding section indicates 
that there was a Sizeable uncertainty in the values of the 
heat leak, but given the size of the heat leaks calculated, 
it is quite clear that these results were not merely 
artifacts of the precision of the experimental measurements. 
Attempts have been made throughout the experimental work to 
modify the operation and also the construction of the 
calorimeter to reduce these heat leaks, and these attempts 
have met some, albeit limited, success. These attempts, 
briefly described, involved the redesign and reconstruction 
of the alee line to the calorimeter in order to increase 
the resistance to conductive heat transfer from the hotter 
inlet line to the calorimeter internals prior to _ the 
exchanger core, the discarding of the heated air fluidised 
bed thermostat as a Surrounding medium for the calorimeter 
vessel in favour of insulating the vessel and operating it 


at room temperature with the process fluid preheating being 
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done with HBTY heating tape, and installation of an aluminum 
foil radiation shield around the inlet line in the 
calorimeter vessel to reduce possible radiative heat 
transfer. 

The primary effect of these modifications was a 
Sizeable decrease in the amount of "background" heat leak, 
Peemarareductionrinevthe sramount-ncof. heaty leakss into: the 
calorimeter boil-off chamber under non-flow conditions. 
Determinations of background heat leaks were done and an 
example is included in Appendix 6 (Run CO2.1 #1). However, 
Since the temperature profile changed from the non-flow 
profile when fluid was flowing through the calorimeter, 
these background values were of limited usefulness in 
correlating the actual heat leak. Eubank et al. (1982), 
whose work has been discussed earlier, did use non-flow 
values of heat leak to correct for heat leaks during flow 
conditions, since their calorimeter was appropriately 
equipped to adequately reproduce, during non-flow 
conditions, the temperature profile of the corresponding run 
with fluid flow. 

A further design or development decision which appears 
to be detrimental to the eventual correlation of heat leaks 
is the decision not to include sufficiently automatic 
control of flow or inlet temperature of the process fluid. 
While such equipment may not be necessary if heat leaks were 
completely negligible, and may not be of any significance if 


the heat leaks are large, analysis of the mechanisms of heat 
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leak may be eaSier if, for example, the temperature change 
could be independently controlled while flow rates and 
pressures were varied through a series of runs. 
Additionally, the reproducibility of the experimental 
conditions and results would be enhanced. 

Because the reproduction of experimental conditions of 
_ temperature and pressure and, hence, flow rates was not 
Ee achieved, it was not possible to directly investigate 
the effects of probable changes in other variables, such as 
small changes in reference fluid level from one run to 
another, due to the necessary periodic recharging of fluid 
to the system. Other examples are changes in atmospheric 
pressure with the corresponding small changes in saturation 
temperature, enthalpy of vapourisation, and resultant 
changes in temperature profile. Tables 4.1 and 4.2 
summarise experimental work, the parameters of which, though 
not exactly identical, are perhaps similar enough to permit 
some comparison and analysis of difference in the results. 
A first impression was that calorimeter performance changed 
after recharging of F-11. Since F-11 levels were carefully 
monitored during the runs, these performance differences 
should not be due to differences in reference fluid levels. 
Casual checks of the F-11 for discoloration or dissolved 
oils indicated minimal change from the state of the original 
F-11 charged to the system, though it may be argued that 
even minimal changes in composition over time could cause 


substantial differences in the enthalpy of vapourisation of 
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Run # 


N24 #5 
IN232 #1 
Basis of 
External 


N2Z.4 #6 
N22 #2 
Basis of 
External 


N2.3 #4 
N23a #5 
N2@2n#4 
Basis of 
External 


NZS, #5 
N273. #6 
Basis of 
External 


N2Z65 #2 
N232 #2 
Basis of 
External 
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Table 4.1 : Comparison of Run Results for N, 
m P P Zk Ak %L Ah 
[ in out in out calc 
22.4 6.975 16.694 S900) 300.5 -14> <-2793.5 
245501 750061R6e7 040U395%. Cor s00zs #19 =2975.9 
Comparison: inlet, outlet conditions similar 
Changes: see next set 
29.2 6.988 6.594 403.3 302.3 -29 -3141.5 
296A FL 6OSSSSE wee Soot es O05. 7aens0de5 -Ame—-3242.9 
Comparison: inlet, outlet conditions similar 


Changes: reconfiguration of vapour lines 
Calorimeter recharged with F-11 


23.6 6.S989aR6rve94 
as above 
as above 

Comparison: inlet, outlet conditions similar 


Changes: possible loss of F-11 over several days 


Soa Te SOs aes = Zon OO 


36.4 6.896806%687 418.1 9304.9 °-=38" -3505.3 
SGDor e600 96a CGRSc9, PS O4R ees 0s ose 64 = 283)1'..0 
Comparison: flow rates and pressures Similar 
Changes: in-line htr current reduced from 
300 mA to 200 mA 
1909) 6.06 Or O04 665095018 299).01 = Ceara aie 
as above 


Comparison: outlet pressures similar 
Changes: inlet pressure reduced 
in-line heater current reduced to 200 mA 
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-4@ 
c2@.2 @.@r See 
a Shs 
neeiteqmae?) to > 


afr 
‘e2pnen: 


Run # 


Cog ein #2 
CO2,2 8:1 
Basis of 
External 


COZ. 1° #35 
CO2..2.43 
Basis of 
External 


COZ. Fil 
COgs2 Ho 
Basis of 
External 


CO2.0. #2 
[LOPS bate BH 
COZe2, 45 
Basis of 
External 


CO seo 
COZ.3 #1 
COZ. 2.43 
Basis of 
External 


CO2.3 #4 


COZ teh): 


Basis of 
External 


COZ. i3. 00 
COZ v2 #1 
CO 2s 2 
Basis of 
External 


Table 4.2 ;: 
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Comparison of Run Results for CO, 


m P Pp gu 7 %L Ah 

f in out in out calc 
Tie) S.409 Sells 398.5 300.2 -3 -4496.3 
Toes 2 aes Oo 404. 7 S06.) = 1.20 2 4680. 2 


Comparison: inlet, outlet pressures similar 
Changes: calorimeter recharged with F-11 


0p 3) Oe eo oe SO Oe 2) 20, Oo aah 
(See Seems oO 50 On 99D a0 
Comparison: flow rates,pressures Similar 
Changes: as above 
EPolimreveluslightlysnighy ineCOZ. 2.23 


-4018.0 
=4£037.8 


<3 hoya Bs pe PRLS TS 5 /fag bls Losi 9 Oe Cts Ne ey hnemalle’ AS I> TI ORD Ge 5G Paik IF es beg | 
as above 
Comparison: pressures Similar 
Changes: as above 
1a Oe te COM ae Ol es SOO poe oOo une see se | O65 
as above 
as above 
Comparison: attempt to repeat C0O2.3 #1 
Changes: none 
Tae Ue Oe OOl. Se 5205 Coes wen UU we =5 *=4315.5 
as above 
. as above 
Comparison: pressures Similar 
Changes: calorimeter recharged with F-i1 
DrLOreLOBCO2. se5 
Oils Sp OOD umes 2 Om eS OO shl 6200. Ue e 5, 3295.6 
65.0. ci oe Ome oS Ole eco ask eae! Ob el hea ool 2 


Comparison: outlet pressures, flowrates similar 


Changes: inlet pressure set higher in CO2.3 #4 
one run prior to C02.3 #4, F-11 was added 

Tie 2 oes OO. gi Sienieo ue Oi euler U2 ecg metal 473.60 

To. oS De) oe OCm ests 7 SO, Omen 12) 2680.2 

TO 3 sO 20) oe oumee OC Our 0d oa ee Gi 96: 6 


Comparison: flowrates and pressures similar 
Changes: probable higher F-11 levels in CO2.2 runs 
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Run # 


CO2.5 #4 
CO2.5* #5 
Basis of 
External 


CO2.5 #6 
CO2.5: #5 
CO2.5. #4 
Basis of 


External 


COZ. 5a 87 
CO2.5 #3 
Basis of 
External 


CO2.5 #8 
CO2.'S* #7 
CO2Z.5: #3 
Basis of 
External 


CO2.6 #1 
CO2.3781 
Basis of 
External 


CO2.6 #2 
CO2.3. 42 
Basis of 
External 


COZ 63 #3 
C02.3 #4 
Basis of 
External 


m Pp P 4g a %L Ah 

£ in out in out calc 
28.4 9.745 9.563 438.1 310.4 -5 -12576.6 
26.3 9.752 9.570 433.1 309.7 -3 -12502.8 


Comparison: #5 is attempt to repeat #4 
Changes: none 


26.0 9.755 9.577 432.6 309.6 -10 -12486.4 


as above 
as above 
Comparison: attempt to repeat after add'n of 
; some F-11 


Changes: addition of some "make-up" F-11 


34.4 9.727 9.405 451.9 313.7 

37.4 9.759 9.403 458.4 314.6 
Comparison: pressures similar 
Changes: F-11 added one run prior to C02.5 #7 


eS Svi2 O02 od. 
+S5te12025%12 


24.8 9.728 9.408 
as above 
as above 

Comparison: pressures similar 
Changes: possible operational or reading errors? 


430.4 309.2 Og 20209 


TAR 20S. 536983  2OSTESSINS ,¥ 29904 ON7S |43005 

1Seule s+ Oo ues 20 OU meee (+25. —4 4a 
Comparison: flowrates, pressures Similar 
Changes: none 

14,019 3). 54 Sem 2601 PSS0R 5S o2994ne+ 1918-42464 5 

14,6 3.526. 93.260 5 396.9 »§ 298.9) 9+24 54116 —5 
Comparison: flowrates, pressures similar 
Changes: none 

INR ODOR 628° IS5e67en 379068 i296.5 baeeo ses 96n:2 

Oe oOo ers Sem SOO lec SO. Uy cou 3293.06 


Comparison: pressures similar 
Changes: none 
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the reference Plo. Other runs showed performance 
variations which cannot be linked to possible variations in 
F-11 level or composition. External causes for these are 
difficult to postulate. 

Corrections for heat leaks. Keeping in mind that the 
reproducibility of these experimental results is limited and 
hence may lead to a correlation, or heat leak correction 


function, that is suspect of some error or at least 


incompleteness, it remains appropriate to attempt a 
correlation or partial explanation, in at least a 
qualitative sense, of the apparent behaviour of the 


calorimeter. In the literature, one finds a variety of 
techniques employed, depending on the type of calorimeter or 
the author's own preference. 

One method, used successfully with glass calorimeters 
measuring Cp of organic vapours (Montgomery and DeVries, 
1942), is a linear correlation of the measured property at 
some condition with the reciprocal of either the flow rate 
or the square of the flow rate, and subsequent extrapolation 
to zero reciprocal flow rate. This implies a similarly 
linear relationship between the heat leak per mole of 
process fluid at a given condition araetnt reciprocal of the 
flow rate, or of the flow rate squared. The method has an 
intuitive appeal in terms of equation 2.12 in which this 
heat leak term approaches zero as the process fluid flow 
rate increases, so long as the actual heat leak itself does 


not increase faster than the flow rate. Of course, 
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variation in flow rates dictates variation of at least the 
pressure drop and these variations from the nominal 
conditions must also be corrected for in the determination 
of the measured value of the property at a given flow rate. 

Judd et al. (1980) used a resistive heating method to 
measure enthalpies of mixing of condensable vapours. A 
correction for heat exchange with the surroundings was made 
in a manner similar to that of Montgomery and DeVries. 
Since the calorimetric method was completely different, it 
is interesting to note that a linear correlation of the 
measured property value, at a given temperature and 
pressure, with reciprocal flow rate was also successfully 
used here for extrapolation to the "corrected" value at zero 
reciprocal flow rate. MThus, thas method appears to have 
broader applicability. 

A correction for heat losses is also extensively 
discussed by van Kasteren and Zeldenrust (1979) who used a 
resistive heating method to measure enthalpies of 
condensable natural gaS mixtures. These workers analysed 
heat losses more explicitly, in terms of conduction and 
radiation due to temperature differences in the calorimeter. 
It can be shown that this correlation of theshéeatceleakomper 
mole of process fluid for a given temperature profile 
involves a proportionality with reciprocal flow rate similar 
to those discussed earlier. Additionally, this 
proportionality can be calculated from a knowledge of the 


temperature profile. There is, however, a certain non-zero 
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value of the heat leak at infinite flow rate for which a 
correction must also be made. 

Sagara et al. (1977) used a calorimeter also patterned 
after that of Sahgal et al. and thus similar in_ some 
respects to the one developed and tested in this work. The 
correlation used is one in which the heat leak is a function 
of the flow rate of the reference fluid and not of the 
process fluid as was used in the techniques discussed above. 
The reference fluid used was liquid N,. This correlation, 
given in equation (4.6), can be shown to be similar to those 
described above which relate the heat leak to the reciprocal 


of the process fluid flow rate. 


10.75 (4.6) 


: = -0,.04416 + 
m A m A 
see pe ge 


If the heat leak is relatively small compared to the total 
measured heat exchange, then the reference fluid flow rate 
is approximately proportional to the process fluid flow 
rate. Thus the heat leak 1S approximately proportional to 
the process fluid flow rate and, following equation (4.6), 
the percentage heat leak is proportional to the reciprocal 
process fluid flow rate. The limits to this correlation are 
that it was said to hold only for boil-off heat exchanges 
below 243 cal/min (1017 J/min corresponding to approx. 180 
mmol N,/min), where the heat leak percentage is equal to 
Zen0,eandvalsopthatywi« can holdsionivsabovesl0.3 cal/min «(43 


J/min» icorresponding,.to -approx.s, 7..7 mmol N»/min) at ywhich 
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point the heat leak percentage would be 100. It is to be 
noted that there was no mention of dependence on inlet 
temperature, though the data for which this correlation was 
used included inlet temperatures ranging from 110 K to 280 
Re 

Applicability of the correlations to this work. In the 
experimental work presented here, significant changes in 
flow rate at a certain nominal pressure could only be 
accomplished with sizeable changes in pressure drop, and 
large changes in the temperature profile were also observed. 
While the process fluid heaters could have been manually 
adjusted to minimise the change in temperature profile as 
flow rates changed, this was not done. The heat leak 
analysis of a similar calorimeter by Sagara et al. appears 
to imply that this is not necessary. In addition, the major 
concern through most of the work was to determine the size 
of the heat leak and to make changes to the calorimeter or 
recommendations for such changes in order to reduce the heat 
leaks to the point where questions of reproducibility really 
become important. As well, the run time that would be 
expended and the change of the base setting of the process 
fluid heater that would be required, if restoration to the 
original inlet temperature were attempted, were 
considerations that discouraged this. 

Hence, the new flow resulted in new conditions so that 
the heat leak could be determined there. Since the 


conditions were not restored to the original nominal values, 
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correlation of the measured property values or heat leak 
with only flow rate varying is not possible here. However, 
following the lead of such a correlating technique, the 
experimental heat leak percentages are plotted versus flow 
rates of N.,, CO,, and the CH,-H,-CO, mixture, respectively, 
in Figures 4.1, 4.2, 4.3. It must be pointed out here that, 
as flow rate increased, so did the inlet temperature (due to 
better heat transfer from the in-line heater as well as in 
the preheating section) in a roughly linear way (Figure 
4.4), regardless of the pressure. Thus temperature was a 
"hidden" variable here, and the value of these plots is not 
to be found in a directly extracted heat leak versus flow 
rate relationship. They do, however, serve to illustrate 
the magnitude of the problem and the trend of the problem 
with flow rate, given the operating strategy used. 

It is encouraging to note in Figure 4.3 that, for the 
gas mixture, the agreement of the experimental results with 
the predictions of the Peng-Robinson equation of state is 
somewhat better than what could be expected based on the 
experience with the N. and CO, tests. 

Most of the neeneds described above deal with positive 
heat leaks which approach zero at high flow rates. Though, 
as previously stated, the functionality of heat leaks with 
flow rate cannot be determined from Figures 4.1, 4.2, and 
4.3 alone, the decreasing trend of heat leak with flow rate 
appears to be borne out in this work too. However, as_ seen 


in these figures, rather than approaching zero at higher 
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flow rates, the heat leak passes through zero and becomes 
negative for many runs. The reasons for this behaviour are 
the subject of the following discussion and of possible 
future correction by design and reconstruction. 

Phenomenology. Obviously, there must be an operating 
mechanism for heat loss from the fluad. after the: inlet 
temperature is measured, which increasingly negates the 
(positive) heat leak into the calorimeter boil-off chamber 
as flow (and with it, inlet temperature) increases. It may 
be useful then to regard the total leak as a combination of 
a relatively constant positive heat leak into the 
calorimeter and a negative term (heat loss from fluid not 
transferred to the vapourising reference fluid) which 
increases with flow rate, at least in the range of 
conditions studied. (It may be that at much larger flow 
rates, some steady negative value will be attained.) Here 
again it should be stated that, had there been more data in 
which the temperature was controlled independently of flow 
rate, the mechanism might have been easier to deduce and 
quantify. However, some discussion which may serve to. shed 
light on the problem follows. 

As the process fluid flow rate increases, the positive 
contribution to heat leak probably increases slightly, since 
the temperature gradient down the inlet line coincidentally 
increases. Thus, conduction down the line to the heat 
exchanger core and eventually to the inner vessel may be 


increased, even though the tubing walls are very thin and 
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provide minimal conduction of heat flow. This is a common 
heat leak path to which others have also referred (eg. 
Bubanktietvaial’., Wwi982)% Obviously though, if there is an 
increase in heat leak into the boil-off chamber, it appears 
to be more than compensated for by some larger increase in 
the heat loss terms. 

It is felt that one of these terms is actually more 
properly described as a matter of accuracy of the 
measurement of the bulk fluid temperature at the inlet. 
Because of the small diameter of the tubing relative to the 
diameter of the thermocouple, lit may be that the 
thermocouple is in partial contact with the wall which is at 
a Slightly lower temperature than the bulk fluid at that 
point. Hence both the inlet temperature used and the 
prediction of the enthalpy change will be too low. Because 
the actual fluid enthalpy at inlet is higher, the actual 
measured enthalpy change is higher and a positive heat leak 
is calculated. As flow rates increase, however, the wall 
temperature will approach the bulk fluid temperature more 
closely, and the thermocouple temperature will more closely 
approximate the bulk temperature so that the positive "heat 
leak" calculated due to this effect will be reduced. This 
contribution will approach zero as flow rate increases and 
will never be negative. 

Perhaps more important is the increase of heat loss 
from the heat exchange core due to radiation and conduction 


as the inlet temperature (and with it, flow rate) increases. 
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If the radiation shields used are not effective, the 
radiative heat transfer to the surroundings increases 
proportionally with the fourth power of the absolute 
temperature, and this could account for some of the 
increasingly negative heat leak with flow rate. 
Calculations done for the inlet line above the inlet 
thermocouple, given in Appendix 1, indicated that radiative 
losses were more significant than conduction down the line, 
but the installation of a foil radiation shield around the 
inlet line produced little noticeable increase in inlet 
temperature. The inlet temperature was obviously not so 
high that radiation had a large effect, so the radiation 
shields may not be that important at the temperatures used 
in this work, except in so far aS improper construction 
could actually increase the heat loss. 

A second loss mechanism could be the conductive heat 
transfer through the dead air space, if the vacuum was not 
sufficiently high to be effective. However, this effect is 
smaller than the radiation effect, especially at the 
temperatures used here. 

To summarise the previous discussion, it suffices to 
says “thate -thes accuracy, of, the calorimeter is severely 
affected by heat leaks, the primary mechanisms of which 
appear to change with flow rate and temperature. For: = this 
calorimeter, then, the correction correlation does not 
depend on flow rate alone. This is not the implication of 


the correlation used by Sagara et al. and thought to be 
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applicable to this work as well. As stated very early in 
the section on design considerations, a calorimeter with 
measurable heat leak may still be usefully accurate if an 
adequate correlation for the correction of this heat leak is 
attainable. At this point, such a correction correlation 
does not appear possible based on flow rate alone and 
further work would have to be done to examine independently 
the SErebt of inlet temperature and of flow rate on the heat 
leak. The advisability of this is doubtful, however, given 
the magnitude of the Heat leaks still occurring, since in 
severe cases the correction term might be as high as one 
third of the measured value of the enthalpy change. rt 
would be best to further improve the performance of the 
calorimeter to the point where, for example, the heat leaks 
are only +5 % of total measured heat exchange, and once that 
has been achieved, to seek to correlate the remaining 
"error". Recommendations for further improvements are given 


in the next section. 
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4.4 Recommendations for Further Work 

Since the calorimeter continues to suffer from sizeable 
heat leaks, further changes will be required for 
Satisfactory performance. The size of the heat leaks should 
be reduced before installation of equipment for better 
overall control would be advisable. Should this not be 
achievable, work on a complete redesign of the apparatus 
should be begun, assuming that it is still deemed desirable 
to obtain this type of data in this department. The 
recommendations, then, which follow are in these three 
areaS: improvement of performance of the existing apparatus; 
improvement of control and range of operation; and design 
and construction of a new device, if necessary. 

Improvement of performance. 

i. The insulation around the outer vessel may be 
removed and the calorimeter submerged ina larger dewar 
filled with F-11 at the same conditions aS are present in 
the inner vessel. This will reduce heat leaks from the 
Surroundings, which can vary Slightly due to variations in 
room temperature. This will involve some changes to the 
inlet line preheating system as well, but these changes are 
discussed further in the next section. 

a Professional advice on radiation shielding 
construction should “be obtained prior to reconstruction of 
the shields. Reconstruction will become increasingly 


necessary as the inlet temperatures increase. 
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3. A vacuum diffusion pump should be installed in 
order to determine if very high vacuum in the outer vessel 
causes any change. 

4, A thinner top supporting plate may reduce heat leak 
by conduction from the hot inlet line. This would involve a 
great deal of alteration and would not be necessary if 
recommendation (1) above were acted upon. 

Improvement of control and _ The 


of) ,0peration. 


range 
following steps should be taken if the actions recommended 
above are successful. They should also be considered during 
the design of anew device if that is deemed necessary or 
desirable. 

Hig Brocesspantluidegheating) should) be done by a 
preheating section and a final trim heating section composed 
of an in-line heater controlled to a desired temperature by 
an automatic temperature controller. (A small diameter 
iron-Constantan thermocouple may be used as an in-line trim 
heater.) Heating capability should be sufficient to reach 
at least 400 °C for all flow rates desired. The inlet line 
should be well insulated, particularly where it runs through 
the surrounding F-11 bath recommended earlier. 

2. Tighter control of process fluid flow rate would be 
possible with either a back pressure controller or a flow 
controller. A larger range of flow rates, for pressure 
drops below a given maximum, would be possible if larger 
diameter tubing were used in the exchanger coil. This would 


allow for more nearly isobaric results at higher flow rates. 
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Sik A mass flow meter could be installed in place of 
the wet test meter for the process fluid flow rate 
determination. This would give instantaneous flow rate 
readings and may be more accurate in some cases but would be 
more difficult to calibrate for some of the mixtures for 
which data are desired. 

4. A Ruska pump, providing Sactietee displacement of 
fluids at variable speeds, might be installed in the system 
along with an appropriate heated mixing manifold with 
pressure ballast tank to allow preparation of mixtures of 
gases and steam. In addition, another pump may be used, in 
lieu of high pressure supply cylinders, to achieve the 
higher pressures when desired. 

5. The existing sampling tees may be connected to a 
G.C. in order to analyse for any composition changes between 
inlet and outlet due to coking as a result of high inlet 
temperatures or due to condensation during cooling in the 
colle 
Design and 


COnStructionm ory 


a new apparatus. It may 
well be found that even after further work to improve 
performance of the existing apparatus, the behaviour remains 
unacceptable for high quality work. If so, the following 
Suggestions relative to the calorimeter redesign might be 
considered in addition to those above. 

1. The boil-off matching idea of Eubank et al. (1982) 


for comparative calorimetry appears to be successful, and 


its application to a future apparatus should be examined. 
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2 The volume of reference fluid in the calorimeter 
should be increased above that of the existing device, 
especially if the higher operating conditions are to be 
investigated accurately. A more directly visual method of 
liquid level assessment would be helpful, but not absolutely 
necessary. 

She The F-11 circulation System, with its’ reflux 
condenser and liquid reservoir, is effective as currently 
operated. For larger vapour flow rates, the vapour line 
diameter should be increased to 1/2 inch (12.7 mm) to 
minimise any pressure drops. 

4. Further radiation losses may be avoided by 
immersing the exchanger coil directly in the reference fluid 
within the boil-off chamber. However, this would result in 
the loss of ability to easily control the outlet temperature 
at a value other than the reference fluid boiling point. 
This may be of concern if heat capacity values, for example, 


rather than enthalpy changes, were desired. 
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Nomenclature 


A, B, C, D - correlation coefficients for a thermocouple 

Con constant pressure heat capacity, J/mol K 

e - internal energy, J/mol 

f, j, k - interpolation constants determined solely from the 
values of h, T, and P given in data tables for N, and 
CO, 

g - acceleration of gravity, 9.8 m/s? 

h - specific enthalpy, J/mol 

I - electric current, amp 

m<- mass, kg 

M - molar flow rate of process fluid, mol/min 

my - molar flow rate of reference fluid, mol/min 

M - thermocouple e.m.f., millivolts 

P - pressure, MPa 

P*, P” - pressure entries in the data tables, just greater 
than and just less than the pressure of interpolation, 
MPa 

q ~ heat transferred per mole of process fluid, J/mol 

QO. - heat leak rate, J/min 

R - electrical resistance, ohm 

S$ - specific entropy, J/mol K 

Ts temperature, K 

T*, T”- - temperature entries in the data tables, just 
greater than and just less than the temperature of 


interpolation, K 
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u - velocity, m/s 

v - specific volume, m?/mol 

V - electrical potential, volt 

xX - composition, mol fraction 

Y - a correlative function of temperature 


z - displacement, m 
%*L - heat leak as a percent of the total heat exchange 


6m - change of mass, kg 

6Q - heat flow into system, J 

6t - change in time, s 

6W - work out of system, J 

Ah - difference in enthalpies of the outlet stream and inlet 
stream, denoted by subscript as either predicted by 
interpolation of each in tables or by P-R equation of 
State, or measured by experiment, J/mol 

\- ~ enthalpy of vapourisation of reference fluid 
(for F-11 at 93 kPa, = 24872 J/mol) 

u =- Joule-Thomson coefficient, K/MPa 


@ - isothermal expansion coefficient, J/mol MPa 
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Appendix 1 - Design Calculations 


Some design calculations are included here as examples 
of the types of calculations necessary for the construction 
of this device. 

A. Calorimeter Size. 

The calorimeter inner vessel was sized to allow a 
"residence time" of about 30 minutes, based on a heat 
exchange rate of approximately 2000 J/min, corresponding to 
a F-11 circulation rate of 11 g/min. The actual heat 
exchange rates ranged from 2% to 25% of the design heat 


exchange rate. 


Q = 2000 J/min ma = Q/r, = oe mol/min = 0.080 mol/min 


Capacity m x 30 min residence time 


2°5"mnol*F-11" (sat. “liq. at 0. 093°MPa) 
Volume required = 2.5 mol x 92.5 cc/mol = 230 cc 


A cylindrical volume of at least 2.5 in (6.4 ecm) high and 
1.35 in (3.4 cm) in radius will contain this amount of F-11. 
Surrounding this inner vessel is a "guard" vessel of 1.75 in 


(4.4 cm) in radius containing F-11 at similar conditions. 
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B. Process fluid line size. 

Based on some initial calculations using a high process 
fluid flow rate which would produce a (high) velocity of 
approx. 100 ft/sec (30.5 m/sec) at 750 °F (400 °C) and 2000 
psia (13.9 MPa), the heat exchange coil tubing size selected 
WaSmUierTG "Il "Oro. DY Oss conn Wald lo. coenm O. De Dy 0. 1 
mm wall). The wall thickness is sufficient for the 
conditions as seen from the following calculation as 
outlined in the Chemical Engineers’ Handbook (Perry and 
Chileon, 197s)" chapters 6. 


maiy.. wall thickness =t = Pise Or De 


m 2*(S * E + P * Y) 


where P is the pressure, O.D. is the outside diameter, and 
S is the maximum allowable stress as per A.S.M.E. 
Code at 750 °F : S = 9100 psi 
E is the weld joint efficiency (Perry, Table 6-5) 
E = 0.8 
Y is the coefficient for stress redistribution due 
to creep at higher temperatures (Perry, Table 6-6) 


Y= 0.4 


Thus, the minimum wall thickness = 0.0031 in (0.079 mm). 
The actual wall thickness is 0.00425 in (0.11 mm). 
Other tubing sizes used in the process fluid line were 


checked in a similar fashion. 
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Due to the flow rates actually used experimentally, the 
velocities in the coil ranged from approx. 10 to 30 ft/s 
(0.3 to 0.9 m/s). Higher velocities could be achieved with 
the use of larger pressure drops, if necessary to ensure 
entrainment of condensation occurring during cooling. 

The exit line tubing was sized so that, as nearly as 
possible, the velocity of the process stream out of the coil 
was maintained until out of the whole calorimeter assembly. 
Pubingmot,0.0355/in O.D.-by 0.00625hin wall (0.90 mm O7D. by 
0.159 mm ward provided an inside diameter of 0.023 in (0.58 
mm), allowing the insertion of a 0.020 in (0.51 mm) dia. 
thermocouple, leaving an annular area similar to the inside 


Gross_sécrronaltiarea of ithemheat éxchange, coil. 


C. Heat transfer calculations. 

The process fluid line entrance to the calorimeter was 
reconstructed based on the following calculations. The 
calculations below are done in British units with S.I. units 


given where appropriate. 


Thermal conductivity 


Stainless steel: k_. = 113 els = 16.3 W/me°C 
hreft’<°F 
Nylon: k, = 1.68 Bturin = 9.24 W/me°C 
hre£t7-°F 
Temperatures 
let line temp. be T, = 320°F = 160°C 
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Inlet constructions 


A) "Capped Cylinder" - well insulated and evacuated 


Dimensions: (in inches) 
Oel25 
OV375 
Disk Cylinder 
r) : aT e at 
q eee BESTE dr q as Ay dx 
; 7% : pm T. 
q° is = “i, r2etetef a1 qe dx s ak oA -f ax 
- + L ss x r 
c c b 
detingy = iniry) Ss he ee 
Kos oe c b 


= “ko 2°Te te (CT, ~ T.) 
(L) q = 0.280499(320 - T,) Btu/hr (2) q°31.1523 = T, - 77°F 
Substitute (1) into (2) to solve for Tos then solve for q. 

+ T, = 295.05 OF. q = 7.0 Btu/hr = 2.05 Watt 


If poorly insulated, this construction may be analysed as a fin with 


uniform cross section and uniform k and he — T = T(x) only 


For this case, 4er D D h oP Pa= pn(OeD a) 
gee ites Bech a Mat AGA | Abra (Oe ne oie T Dee) 
dx x x 4 
General sol'n: de Cees + Coeu 
For T= 77% and T, , = 320% with h, = 2=“—— ; 
om ia hr ft* + OF 
=-1 
the constants are: m = 0.34 in A C, = -365 F » and C, = 365% 
Thus, T-T_= 65 (ete a1) OF 
7.5 na) 
: a ez 2 fe 0.34x _ 
Then, dee, bP T dx hea 033375.) 365 Ce 1)dx 
; 0 Boe ibe ape 
Ge enée 6.0 G0 13h mist OMe S600 Watt = 0.8 Watt 
For a more severe estimate, use he = 20 AE : 
hreft~° OF 


1.08x 


- = =_ 0 : = 
Then T- T_ = 60(e 1)*F = dein 6.6 Watt 


However, this value is a high estimate since the construction is insulated. 
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Inlet constructions (cont'd) 


B) Bulkhead Reducing Union in Nylon Bushing in Lid 


Dimensions: | 


(in inches) 


pmies: 


St rots UnLOn Nylon Bushin 


: dT s = 8 a @ © onl 
q=- Kos Zemeret dr q = k SIZ BO A 
re . T "s fi 
° dr h ° dr : 
af = -k aenevef dt q°J, — = -k aemeeef at 
is ss . rt n T 
Fs Ss b b 
q°(in Ea In rs) -q°(ln ie ie tr) hast 
Kkee2en ¢f Ss 
n 


= -k j2emete(T, - 7.) 


(1) q = 2.6675(320 - T,) Btu/hr (2) 4°17.098 = T, - 77°F 


Substitute, G) into (2)) to .solve for T. and q. 


ale =i 7nr; q = 13.91 Btu/hr = 4.07 Watt 


C) Bulkhead Union in Lid 


Dimensions: 


~ le 0.375 


Ow3)5 


(in inches) 


AN 
Siro teen LOR only ee 


° dT 
q = ha emer et dr 


Similar to the above, 
q°(ln tees In vr) = “ko 2° Tete (T, - Ty) 
q = 1.683°(243) Btu/hr = 409 Btu/hr = 120 Watt 


Thus, of the inlet constructions examined, construction A appears 


to provide the minimum heat loss to the calorimeter lid. 
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Some inlet line heat leak calculations 


The estimated radiative heat losses from a section of the 
inlet line inside the calorimeter are compared to the 
corresponding conductive heat losses for two relatively low 
temperature cases. It is found that even for these cases 
radiation is more important than conduction, and this will be 


even more true at higher operating temperatures. 


OF065.s1'n OF 1651 em 
0.047 in = 0.1194 cm 


Dimensions: Tube O.D. 


Dh 
Tube section length (distance between 
measuring thermocouples) = 2 in = 5.1 cm 
Physical Constants and Properties: 


keg = 113 Beulinyneeee se oF = 16.3 W/meK 


Q 
| 


and view factor) 


A. Conductive heat transfer 


1) 1, = 3d2 Ke nam 08.5. K 3) (At sg-3 7k 


1 2 
= gies ee 
Shae K gies dx aga a Ax 
. Me el eR So OOheaeps (2406uG. 050 luack 
0 = 1.14 x 1073 Watt 


HS Wr 4333.05K 07, = 832958 K 3) AD meta 00K 


Similar to the above, 


Q = 1.31 x 10me Watt 


B. Radiative heat transfer (without shield) 


i) let average tube em temp be T, = 310 K ey = 295 K 
e 4 Surr 
Q = TeDeAxeerFeoe g - diate) 


habe igeration Gh 03.088, Gg) Watt 


0 = 2.49 x Wee Watt 


17) Let TY = 330 K; similar to the above, Q = 6.43 x ioe Watt 


= 1.355 x lene Peies cuss = 5.669 x 107° Wee 


1; F=1 £4(conservative values for emissivity 
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Appendix 2 - Instrument Calibrations 


A. Calibration of Thermocouples. 


The thermocouples were calibrated in the Instrument Shop and 


the calibration results for each thermocouple were 


correlated by least squares fit to an equation of the 


following 


where Y is one 


A summary 


Sw. TAG 
No. No. 
1 1 
Z 6 
3 8 
4 n/a 
5 23 
6 4 
a. 5 
8 a 
9 9 
10 = 
al n/a 
t2 n/a 


form: 


laeve =ael 
27 ee 
Si) Ye =e LN els) 
Seo =i) 


of 7/C locations and 


Location 


Mass Flow Meter 
Freon Return 
Proceekhluldn supply 
Cal. Inlet -*Upper 
Cal. Inlet - Lower 
Cal. Outlet 

Proc. Fluid Prht. 
Not in use 
Refrigeration Unit 
Not in use 

M.F.M. Signal 

P.D. Trans. Signal 


Y= A + BeXM + C/M + D*ln(M) 


Cheatin cy 


CT i eek) 


CT nec) 


CT oie eke) 
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correlations follows: 


G32) 


of the following functions of temperature: 


(3n3a) 
3630) 
(shagekey, 
ei aeishy, 


Correlation Coefficients 


A 


4.093 


eke 


0.866 
3059 
3.004 
Osise 
2.806 


6.038 


B 


21.998 
18.0541 
I7%8)7 
02003 
0.002 
17.856 
18.044 
n/a 
OT Ba 
n/a 
n/a 
n/a 


C 


-6.543 
=O 050 

0.939 
-0.064 
70.026 

Ve0de/ 
LE AST =! 


-0.064 


D 


SCZ. 6 
7 13 
2.662 
0.936 
On 200 
PE he, 
Te ctERG 


0.942 
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The maximum uncertainties in temperature arise from the 
correlation and are at most +0.1 K in the range of 
temperatures used. 

B. Calibration of Pressure Gauges. 
In the calibration of the two Heise gauges used, any 
deviations greater than 0.1% of full scale were noted. In 
the range of pressures used in this work, there were no such 
deviations noted. Thus, inlet pressures could be reported 
with a maximum uncertainty of 20 kPa and the outlet 
pressures could be reported with a maximum uncertainty of 10 
kPa. 

C. Calibration of Flow Meters. 
The mass flow meter calibration curves appear in Figures 
A2.1 and A2.2. The linear least squares correlation of the 


calibration data for each meter are as follows: 


BROOKS Model No. 5810 BIF2E5A - (Flow meter A) 
Flow rate (g/min) = 0.0062 + 0.2019*V 
BROOKS Model No. 5810 B1iJ2E5A - (Flow meter B) 


Flow rate (g/min) = 0.0611 + 1.6988*V 


In the ranges used, the maximum lnc®eteint es due to the 
correlation are 25 mg/min (0.18 mmol/min) for flow meter A, 
and 100 mg/min (0.73 mmol/min) for flow meter B. 

The Precision Scientific wet test meter was calibrated and 
checked periodically. The accuracy is 0.5% or 0.0005 cu ft 


per 0.1 cu ft of gas passed through the meter during a test. 
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Appendix 3 - Calculation Routines 


C CONTENTS OF FILE CALOR 
c 
C CALOR READS DATA FOR TEST GAS (IF AWY) AND FOR F-11 
C AS WELL AS EXPERIMENTAL DATA FOR THE CALORIMETER. 
C ENTHALPY CHANGES ARE CALCULATED FOR EACH TEST INTERVAL AND 
C THE RESULTS AND EXPERIMENTAL INFO STORED IN ARRAY “DHCAL”. 
C MAIN PROGRAM READS DATA AND CALLS SUBROUTINE CALC TO DO 
C MECESSARY CALCULATIONS AFTER WHICH MAIN WRITES RESULTS 
C INTO GCUTPUT FILE OR TO THE TERMINAL. 
C LOGICAL UMIT WOS. ARE: 1 - M2 ENTHALPY DATA 
is 2 - €O2 ENMTHALPY DATA 
is 3 ° PREOM DATA 
c 4 - EXPERIMENTAL DATA 
c S - TERMIWAL 
e€ G@ ° TERMIMAL OR RESULTS FILE 
c 
is 
REAL HT(20),HP(&0),HH(20,40) ,FRD(10,4) 
REAL EX(10,20) ,DHCAL(10,10, 20). 
INTEGER GAS(3)/‘°N2’ ,°CO2’,’TEST’/ ,MS$G(20) 
COMMON HT,HP,HH,FRO,EX,DHCAL ,WPTS ,NPF 
600 13 181,10 
00 13 v#1,10 
ODO 13 K#1,17 
13 OHCAL(1,45,K)80.0 
READ(3,101) NPF 
READ(3,103)((FRO(I,J),J381,4),181, MPF) 
READ(4,101) KGAS 
1F (KGAS.EO.3) GO TO 10 x 
READ(KGAS,101) WNTH,NPR 
READ(KGAS ,102) (HT(1),181,WTH) 
READ(KGAS ,102) (HP(J),J81,NPH) 
READ(KGAS,102) ((HH(1,J3),181,NTH), S81, BPH) 
10 READ(4,101) NRUNS 
00 100 KRUN=!1, WRUNS 
READ(4,101) MPTS,KTYP 
00 15 J#i,17 
READ(4,102) (EX(I,J3),121,NPTS) 
c WRITE(6,102) (O2x(1,J4),121,NPTS) 
DO 14 18#1,NPTS 
IF (E€X(1,J5).80.0.0.AND.1.6T.1) Ex(1,JSdPEX(Teoi,d) 
14 COWTINUE 
18 CONTINUE 
is READ(4,130)(MSG(J),J21,20) 
Ci3O FORMAT (2044) 
CALL CALC(NTH,NPH,KRUN ,KGAS ,KTYP) 
NINTSNPTS-1 
IF(KTYP.EO.0) GO TO 90 
TF(KGAS.EQO.3) GS TO BO 
WRITE(6,105) KRUN,GAS(KGAS) ,(J.J#1,NINT) 
WRITE(6, 107) (DMCAL(KRUN,J,8),J21,NINT) 
WRITE(6,108)(OMCAL(KRUN,J,9),JSt1,MINT) 
WRITE(6,109) (DWCAL (KRUN,J,10) ,J81,NINT) 
WRITE(6,110) (DMCAL(KRUN,.J,11),J5S21,NINT) 
WRITE(6, 111) (DOMCAL(KRUN,J,13),J81,MINT) 
WRITE(6, 124) (DHCAL(KRUN,J,16),J3=1,NINT) 
WRITE(6, 129) (DHCAL(KRUN,J,17),J381,NINT) 
WRITE(6,. 131) (OMCAL(KRUN,J,12),J381,NMINT) 
WRITE(6,112)(DHCAL(KRUN,J,1),J381,NINT) 
WRITE(6, 113) (DHCAL(KRUN,J,2), 981, MINT) 
WRITE(6, 114) (OMCAL(KRUN,J,3),J581,NIWNT) 
WRITE(6, 118) (DOWCAL(KRUN,J,6),J5E1,NINT) 
WRITE(6,1186)(OHCAL(KRUN,J,4),J81,NMINT) 
WRITE(6,128) 
WRITE(6,.117) GAS(KGAS), (OHCAL(KRUW,J,5),J381,NINT) 
WRITE(6,118) (DOMCAL(KRUN,J,7),J81,NINT) 
WRITE(6,119)(OMCAL(KRUN,J,.15),JS81,NIMT) 
WRITE(6,120)(DMCAL(KRUN,J.14),J81,MINT) 
WRITE(6, 128) (DNCAL(KRUN,J,18),J81,NINT) 
WRITE(6,126)(DMCAL(KRUN,J,19),J81,NINT) 
WRITE(6, 127) (DHCAL(KRUN,J,20) ,J21,MINT) 
ce WRITE(G6,132)(MSG(4N),581,20) 
GO TO 100 
80 WRITE(6,105) KRUN,GAS(KGAS) ,(J,J81,WINT) 
WRITE(6, 107) (DMCAL(KRUN,J,8),JE1,NINT) 
WRITE(6,. 108) (OHCAL(KRUN,J,9),J81,NINT) 
WRITE(6, 109) (DHCAL(KRUN,J,10),J21,NIWNT) 
WRITE(6,110)(DHCAL(KRUN,J,11),Je1,MINT) 
WRITE(6,111)(OWHCAL(KRUN,J,13),J381,NINT) 
WRITE(6,124)(DHCAL(KRUN,J,16),J81,MIWT) 
WRITE(6, 133) (OHCAL(KRUN,J,2),J81,MINT) 
WRITE(6, 134) (OMCAL(KRUN,J,6),J81,MINT) 
WRITE(6,116)(OHCAL(KRUN,J,4),J81,NMINT) 
WRITE(6, 138) 
WRITE(6,117) GAS(KGAS) , (DHCAL(KRUN,J,S),J3t1,NINT) 
WRITE(6, 118) (DOMCAL(KRUW,J,7),J81,NINT) 
WRITE(6,118)(DHCAL(KRUN,.J,18) ,J#1,NINT) 
WRITE(6,120)(DHCAL(KRUN,.J,14),381,NINT) 
WRITE(6,125)(DHCAL(KRUN,J,18),J381,NINT) 
WRITE(6,126)(OMCAL(KRUN,.J,19),J381,NINT) 
WRITE(6,127) (OHCAL(KRUN,J,20),J81,NINT) 
c WRITE(6, 132) (MSG(J),J81,20) 
GO TO 100 
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WRITE(6,121)KRUN, (J,J=1,NINT) 
WRITE(6,122)(DHCAL(KRUN,J,1),J=1 
WRITE(6,107) (DHCAL(KRUN,J,5S),J=1 
WRITE(6,108) (DHCAL(KRUN,J,6),J=1,NINT) 
WRITE(6,123)(OHCAL(KRUN,J,9),J=1 
WRITE(6,119)(DHCAL(KRUN,J,8),J=1,NINT) 
WRITE(6,120)(DHCAL(KRUN,J,10),J=1,NINT) 
WRITE(6,125) (DHCAL(KRUN,J,18),J=1,NINT) 
WRITE(6,126)(DHCAL(KRUN,J,19),J21,NINT) 
WRITE(6,127) (DHCAL(KRUN,J,20),J=1,NINT) 
WRITE(6,132)(MSG(J),J3=1,20) 

FORMAT(/SX, ‘COMMENTS: ‘,20A4) 

CONTINUE 

FORMAT (214) 

FORMAT(10F7.0) 

FORMAT (4F10.0) 
FORMAT(//SX,’RUN’,13,716,’TEST GAS - *,A4,/ 


@€5x,’CONDITIONS FOR TEST INTERVAL’ ,9(4X,11,4X)) 


FORMAT(/5X, °TEMPERATURE’,6X, ‘Upper Inmiet’,9F9.2) 
FORMAT(/SX, ’ TEMPERATURE’ ,6X,’Lower Iniet’,9F9.2) 
FORMAT(T8,’(K)’,T23, ’Outiet’,1T34,9F9.2) 
FORMAT(S5X, ‘PRESSURE’ ,T23,’Inlet’,T34,9F9.3) 
FORMAT(T&,’(MPa)‘,T23, ‘Outiet’,T34,9F9.3) 
FORMAT(T23, ‘Atmospheric’ ,9F9.3) 

FORMAT(/,5X, ‘HEAT EXCHANGE (J/min)’,T34,9F9.3) 


FORMAT(SX, ’ENTHALPY CHANGE Predicted 2 Hat ewe ae 
FORMAT(T8&,‘(J/mol1)’,T23, ‘Measured oak) it Aa 30) 
FORMAT(/,5X, ‘FLOW RATE Rian Pred. Mol. ’,9FS.5) 
FORMAT(T&,’(1/min) or’,T23,’Meas. Mol. ‘,9FS.5) 
FORMAT(T8&,’(mmol/min)’,1T23,’Calc. Voi. ‘’,9F9.5) 
FORMAT(T17,A4,T23,’‘Meas. Vol. ’,9FS9.5) 
EORMAT(T23, “Catc.. Mol..*,9F9.5) 
FORMAT(/,5X,’Freon Sat’’n Temp. (K)‘,T34,9F9.2) 
FORMAT(5X,’Freom Meter Temp. (K)’,1T34,9F9.2,//) 


FORMAT(//5X,°RUN’,I3,T16,’BACKGROUND TEST (NO FLOW) ‘’/ 


&@Sx,’°CONDITIONS FOR TEST INTERWAL’,9(4X,11,4X)) 


FORMAT(/5X,’TEST HEAT LEAK (J/min)’,1T34,9F9.3) 


FORMAT(5xX,’‘Atmospheric Pressure (MPa)’,T34,3FS.3) 
FORMAT(SX,’Freom Return Temp. (K)‘,T3&,9FS.2) 
FORMAT(5X,’‘Heat Exch. Corr. (J/min)‘,T34,9FS.3) 
FORMAT(SX,’Ambient Air Temp. (K)’,T34,9F9.2) 
FORMAT(T32, ‘cccrccesccseeseresese?’) 

FORMAT(SX, ‘HEAT LEAK (J/min)’,734,9FS.3) 
FORMAT(S5X, ‘RATIO: LEAK OVER EXCHANGE’ ,7T34,9F9.4) 
FORMAT(SX, ‘MEAS. ENTH. CHANGE (J/mo1) OCS hray) 
FORMAT(/,5X, ‘FLOW RATE Fast Meas. Mol. ’,9F9.58) 
RORMAQ We naO ri n/m Vila 2), 5 er aawle Slee areola) sielecel= = om) 
STOP 

END 


SUBROUTINE CALC(NTH,NPH,KRUN,KGAS,KTYP) 

REAL ZN(3) 

REAL HT(20),HP(40),HH(20,40) ,FRD(10,4) 

DIMENSION EX(10,20),DHCAL(10,10,20) 

COMMON HT,HP,HH,FRD,EX,DHCAL,NPTS,NPF 

ZN(1)30.9998 

Z2N(2)20.9951 

2N(3)20,9977 

NINT=NPTS=-1 

BO 200 [21,NINT 

CTZEX. CLS Tt) Ze or | t): 

TA=(EX(1I,6)*EX(141,6))/2. 
CT1=3.0043+4 .002#TA- .O0265/TA+ .SS482ALOG(TA) 
TIZEXP(CT1)+273.15 

TB GEM Git TeX OL hy Dis a 
T2=.7321+17.8S55*TB*1.0768/TB+2.7298*ALOG(TB)+273.15 
TCS CEX( 1, 8) 4x 141,38) )/2. 

TS=2.8064+18 .0442*TC-1.2677/TC+1.1165*ALOG(TC)+273.15 
TOH(EX(1,2)+EX(1+1,2))/2. 

TFM=4.0929421.9982*TD-6 .5430/TD-8.278S5*ALOG(TD)+273.15 
TEZCEX(1, 3) +Ex¢1+1,3))72. 


TRLIQ=1.7214418,.O513*TE-0.0846/TE+1.7I3J3*ALOG(TE)+273.15 


TAMB=((EX(1,17)4EX(1%1,17))/2°32.)25/94273.15 
FYCOR=1. 

PASEX(1,15)*0.1/750.062*FTCOR 

JPROP=2 

CALL PROPi1(PA,HVAP,NPF,FRD, JPROP) 
HVAP=HVAP+24000.0 

JPROP=3 . 

CALL PROP1(PA,VSP,NPF,FRD,JPROP ) 
VSPF=VSP#0.137369 

VOLT=EX(I+1,11) 


FOR LARGE MFM B=0.0611 , M=1.6988 
FOR SMALL MFM B=0.0062 , M=0.2019 


RFMOL=( .0062+.2019#VOLT)/137.37 
RFVOLERFMOL#VSPF 

JPROP=4 

CALL PROP1(PA,TSAT,NPF,FRO, JPROP ) 
TSAT=TSAT+273.15 

FHTCAP=120.7 
SATCOR=RFMOL#(TSAT-TRLIQ) *FHTCAP 


SATCOR TO BE ADDED TO DHCAL(16) TO CORRECT FOR DIFF FROM 
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IF (KTYP.E0.0) GO TO 210 
PIZ(EX(1,13)4EX(1%1,13))/2.*0.1/14.5038 >PA 
P2=(EX(1,18)+EX(141,14))/2./1000. PA 
IF (KGAS.EO.3) GO TO 198 
CALL PROP2(T1,P1,H1,NTH,NPH,HT,HP, HH) 
CALL PROP2(T2,P2,H2,NTH,NPH,HT,HP, HH) 
DHCAL(KRUN,I,1)2H2-H1 
GO TO 199 
198 DHCAL(KRUN,I,1)2-98899. 
199 TWTME((EX(1,17)eEX(141,17))/2.-32.)25./9.%273.18 
R=8.314E-6 
VSPN=ZN(KGAS) =R*STWTM/PA 
FWTM=s1. 
RNVOL=(EX(1¢1,12)-EX(1,12))20,0283168=FWTM/DT 
RNMOLERNVOL/VSPN 
DHCAL(KRUN,1,2)2°-REMOL=HVAP/RNMOL 
IF (KGAS.NE.3) DHCAL(KRUN,1I,3)=RNMOL*(H1-°H2)/HVAP#1000. 
DHCAL(KRUN,1,4)2#RFVOL21000. 
DHCAL(KRUN,1,5)=RNVOL#1000. 
OHCAL(KRUN,I,6)=2RFMOL=1000. 
DHCAL(KRUN,1,7)=RNMOL=1000. 
DHCAL(KRUN,1,8)27T1 
DHCAL(KRUN,I,9)2T2 
DHCAL(KRUN,1,10)2P1 
DHCAL(KRUN,I,11)2P2 
DHCAL(KRUN,1,13)=PA 
DHCAL(KRUN,1,14)27FM 
DHCAL(KRUN,1,15)2T7SAT 
DMCAL(KRUN,1,16)=2RFMOL=HVAP 
IF (KGAS.NE.3) DHCAL(KRUN,1,17)2RFMOLSHVAP-RNMOL#(H1-H2) 
IF (KGAS.NE.3) DHCAL(KRUN,I,12)=0HCAL(KRUN, 1,17) /DHCAL(KRUN,1,16) 
DHCAL(KRUN,1,18)}2TRLIO ; 
DHCAL(KRUN,1,19)2SATCOR 
DHCAL(KRUN,1,20)27AMB 
GO TO 202 
210 DNCAL(KRUN, 
DHCAL(KRUN,, 
DHCAL(KRUN, 


, 1) =2REMOL*HVAP 
,2)SREMOL®1000. 
»3)=RFVOL=1000. 


I 

I 

I 
DHCAL(KRUN,1,&)2TS 
DHCAL(KRUN,1,5)27T1 
DHCAL(KRUN,1,6)2T2 
OHCAL(KRUN,1,8)2TSAT 
DHCAL(KRUN,1,9)2PA 
DHCAL(KRUN,1,10)2TFM 
DHCAL(KRUN,I,18)=TRLIO 
DHCAL(KRUN,1,19)2SATCOR 
DHCAL(KRUN,1,20)=TAMB 


202 CONTINUE 
00 201 K#i1,20 
201 DHCAL(KRUN,10,K)2DHEOAL(KRUN,10,K)*DHCAL(KRUN,1,K)20T 
200 CONTINUE 
00 208 K#ei,20 
2058 DHCAL(KRUN,10,K)2O0HCAL(KRUN,10,K)/EX(NPTS, 1) 
RETURN 
DT in. min., PT, P2,PA jin MPa., 11,72, TEM, TWIM, etc. in K 
R im MJ/molcK, YSPN,VSPF in cu.m/moi 
RFVOL,RNVOL ith cu.m/min, RFEMOL,RNMOL in meoil/min 
Hi,H2,HVAP im J/moi 
END 


SUBROUTINE PROP2(X,Y,Z,NK,NY,AX,AY,22) 
REAL AX(20),AY(40),22(20,40) 
BO 300 J121,NX 
IF (AX(1).LT.X) GO TO 360 
KeI 
GO TO 301 

300 CONTINUE 

301 06 302 J=1,NY 
EF. AAV (I) CLT. YS) GOL TO, 392 
L2J 
GO TO 363 

302 CONTINUE 

303 Re C220Ky bie 22K ba SAY Ati AY (c= 1) ))) 
AzRe(VYoAV(L-=1))4ZZ2(K,L°-1) 
REA Ked Led 2K) yy ee) AAA Yo Kalter AYICIL = s1:3)1); 
BeRs(VY-AV(L*1))422(K-4,b°1) 
R=(A°B)/(AX(K)*AK(K-1)) 
ZERE(K-AX(K-1))%8 
RETURN 
END 


SUBROUTINE PROP1(X,Y,L,YY,uSP) 
REAL YY(10,4) 
00 400 I12i,L 
TE sOYiY GF wala ble X) GO TO 400 
M=eI 
GO TO 40F 
400 CONTINUE 
401 R=(YY(M,JP)-YY(M°1,J5P))/(YY(M, 1)°VY(Me1,1)) 
YeRe(X-VYY(M,1))*YY(M, JP) 
RETURN 
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THESE SUBROUTINES ARE USED IN PLACE OF SUBROUTINE PROP2 
IN THE PROGRAM “CALOR” IF SPLINE INTERPOLATION OF THE 
ENTHALPY DATA IS DESIRED 


SUBROUTINE PROP2(X,Y,Z,NX,NY,AX,AY,AZ) 
REAL AX(20),AY(40),AZ(20,40),B(4),C(4),D(4),XX(&),YY(8&) 
REAL 22(4),S(4) 
INTEGER INDX(4),INDY(4) 
DQ 300 I=1,NX 
IF (AX(I).LT.X) GO TO 300 
Ket 
GO TO 30% 

300 CONTINUE 

301 DO 302 yei,NY 
IF (AY(J).LT.Y¥) GO TO 302 
Leu 
GO TO 303 

302 CONTINUE 

303 IF (K.LT.3) GO TO 305 
IF (K.EQ.NX) GO TO 307 
DO 304 I=1,4 

304 INDX(I)=K-°341 
GO TO 310 

305 DO 306 121,4 

306 INDX(I)s1 
GO TO 310 

307 DO 308 I=1,4 

308 INDX(I)sK-4+I 


310 IF (L.LT.3) GO TO 315 
IF (L.EQ.NY) GO TO 317 
DO 318 J=1,4 

316 INDY(J)2#L-344 
GO TO 320 

315 DO 316 J=1,8 

316 INDY(J)2J 
GO TO 320 

317 00 318 J=i,4 

318 INDY(J)sL-44u 


320 DO 340 Jti,4 
DO 330 I121,4 
KXX(1)SAX(CINDX(I)) 
ZZ(1)SAZ(INOX(T), INDY(J)) 


CORRECTION FOR TWO VALUES IN CO2 TABLES IN LIQ. PHASE 
IF (INDY(J).E90.7.AND.INOX(1I).EQO.1) XX(1)5291.55 

IF (INDY(J).EQ0.7.AND.INDX(1).E9O.1) 22(1)231849. 

IF (INDY(J).EQ.8.AND.INOX(1).EO.1) XX(1)2295.25 

IF (INOY(J).EQO.8.AND.INDX(1).EQ.1) 22(1)231096 


330 CONTINUE 
NZ@ 
CALL SPLINE(N,XX,22,B8,C,D) 
S(J)®SEVAL(N,X,XX,22,B8,€,0) 
340 YY(J)SAY(INDY(J)) 


CALL SPLINE(N,YY,S,8,C€,D0) 
ZESEVAL(N,Y,YY.S,B,€,D) 
RETURN 
END 
SPLINE INTERPOLATION FROM FORSYTNE,MALCOLM&MOLER 
MASLIYAH.... CHE. 678 JAN. 1978 
SUBROUTINE SPLINE (N,X,Y,B,€,0) 
REAL X(N),.Y(N).T,B(N),C(N),D(N) 
INTEGER IB, I 
NM1=N-7 
IF ( N .LT. 2 ) RETURN 
lf N UUr. 3) 7) GO to sso 
O04) 8 G2) XA 
e(2) = (Y¥(2) *¥(1))/0¢1) 
po 510 I = 2, NMi 
O(I) = X¢I%1) + X(T) 
B(I) = 2.8(D(Ie1) + D(I)) 
COLS ty 2 CVS 1) “= V-CE)) 7 DCT) 
CUD). =O Ul+4) es CCR) 
$10 CONTINUE 
B(1) = °D(1) 
B(N) = -D(Ne1) 
Cito sxo.. 
C(N)#O. 
IF ( N .EQ. 3 ) GO TO S15 
C(i) © C(3)/(K(8)-K(2)) © CC 2)/(K(3)-K(1)) 
C(N) = C(N21)/ (KON) °X(N-2)) © CON-2)/ (KEN 1)°K(N=3)) 
C(i) = C(1)#O0(1)#82/(K(8) °K (1)) 
C(N) & -C(N)*®D(N-1)#22/(KX(N)-X(N-3) ) 
$15 DO 520 I = 2, N 
T = D(I-1)/B(I-1) 
B(I) = B(I) - T#eO(I-1) 
CUT) se CUR) <9 TS Cire) 
520 CONTINUE 
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C(N) = C(N)/B(N) 
DO 530 IB = 1, NM1 
I = N-IB 
C(l) = (C(I) - O(I)#C(1+1))/B(T) 
530 CONTINUE 


B(N) = (Y(N) °Y(NM1))/D(NM1) + D(NM1)*(C(MM1) + 2.2E6(N)) 
bDOo S40 I = 1, NM1 
B(I) = (Y( E41) -VCI))/D(T) - DCE) F(E(T+t) + 2.2C(1)) 
D(I) & (C(1+1) - C(I))/D(1) 
C(r) & 3.#€(1) 


$840 CONTINUE 
C(N) = 3.#C(N) 
D(N) = D(N-1) 


RETURN 
550 B(i) = (CY¥(2)SVC1) VS CK(2)°K(1)) 
c(i) = ©. 
O(i1) = 6. 
B(2) = B(1) 
c(2) 2 6. 
0(2) = 0. 
RETURN 
END 


FUNCTION SEVAL(N,U,X,Y, 
REAL U,X(N),Y(N),B(N),C 
DOUBLE PRECISION OX 
DATA I/1/ 

IF ( I .GE. W ) 
TE € @ VET. C1 
TE GC Ue EES oXiG? 

410 I # 1 
J = Ni 

420 K = (14J3)/2 


~w 


) G6 TO 410 
) 


IF ( U .LT. X (KD) ) J & K 
Fe ¢ U SGE. & (KY ) Fe & 
IF ( Jd .GT. I + 1 ) GO TS 420 


830 DXsUeX(1) 
SEVAL = Y¥(1)*DxX2(B(1)*¢OX*(E(1)+DK*D(1)) >) 
DERI VSB(1)+OX2(2.08€(1)+DKX#(3.08D0(1))) 
RETURN 
END 
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Appendix 4 - Fluid Properties 


‘Table A4.1 contains the relevant properties of Freon-11, in 
the range of conditions occurring in the experiments. Values 
of enthalpy as a function of temperature and pressure for 
nitrogen and carbon dioxide (used as standards) are given in 
Tables A4.2 and A4.3. 

Enthalpy values for the test gas were estimated using 
the Equi-Phase program which employs the Peng-Robinson 


equation of state. 


i 
TABLE A&.1 : FREON-11 DATA CONVERTED FROM ASHRAE FUNDAMENTALS 18977 
STORED IN SIMILAR FORM IN FILE DFREON 


P IN MPA, HVAP (ENTHALPY OF VAP’N) IN J/MOL 
7 IN €, V (SP. VOL. VAPOUR) IN CU.M./KS 


P HVAP v T (SAT) 
0.088384 248943.326 0.19366 20.0 

0.080183 28915.4480 0.139006 20.586 
©.092010 24887.4616 0.18655 24.171 
©.093865 24859.256 0.18311 21.667 
6.095787 24830.820 0.17874 22.222 
0.097664  24803.072 0.17646 22.778 
0.099608 24775.048 0.17324 23.333 
0.10188 24746.613 0.17010 23.889 
0.10368 24718.178 0.16703 24.464 
©. 10862 24689.742 ©.16403 25.0 
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20. 


.00 


-50 


-00 


- 50 


. 00 


-50 


. 00 


.50 


.- 00 


.50 


-00 


.50 


00 


-590 


-0°0 


.50 


-00 


-50 


. 00 


-00 


eo 


-1582 


- 1626 


°-1670 


cE Yer A | 


- 1756 


me He A) 


- 1928 


- 1965 


° 2044 


- 2083 


oe Nauk 


°21588 


-2198 


2231 


°2845 


°2774 


TABLE A&.2 


(ANGUS ET AL., 


NITROGEN ENTHALPY DATA FROM I.U P.A.C. 


i979) 


ENTHALPIES GIVEN IN J/MOL 


-68 


-98 


-1268 


MH Ye f 


785 


cir Bh he} 


-241 


-269 


-296 


Me se 


°-348 


-374 


-399 


-424 


-448 


-471 


-495 


-760 


-8S5 


350 


1468 


1447 


1427 


1406 


1386 


1367 


1347 


1328 


1309 


1291 


1273 


1285 


1238 


1220 


1204 


1187 


Tigis 


1156 


1140 


1006 


905 


STORED 


TEMPERATURE IN 


400 


2942 


2928 


2914 


2900 


2887 


2873 


2860 


2848 


2835 


2823 


2811 


2799 


2788 


pa Be Ae f 


2766 


2755 


2745 


2735 


2725 


2640 


2581 


450 


4418 


44093 


4800 


4391 


4382 


48374 


4366 


4388 


4350 


4342 


4335 


4328 


4321 


8314 


8307 


4301 


42935 


4289 


4283 


4236 


4208 


500 


589s 


$894 


588g 


5884 


5879 


$874 


5870 


5866 


5861 


S8sS7 


58864 


5850 


S846 


5843 


$840 


$837 


5834 


$831 


$828 


$8ii 


s8og 


IN FILE N2DATA 


K 


550 


7391 


7388 


7386 


7384 


7382 


7381 


7378 


7378 


7376 


7375 


T7374 


7373 


7372 


7372 


7371 


7371 


T3790 


7370 


7370 


7377 


7335 


600 


8894 


8894 


8894 


88g5S 


6895 


8896 


88a7 


8898 


8839 


8900 


8901 


8302 


8904 


830s 


8907 


8909 


831i 


8912 


8sis 


8941 


8976 


6S0 


10411 


10414 


10416 


10418 


10421 


160424 


10826 


10829 


10432 


160435 


10438 


10441 


10444 


10848 


10451 


10455 


104858 


10462 


10466 


10808 


10558 


700 


11944 


11948 


11982 


11986 


11960 


11965 


11969 


11974 


11978 


11983 


11987 


11992 


11997 


12002 


12007 


12012 


12017 


12022 


12027 


12083 


12146 


dy 


aearr 


teers 


ee@e\% 


ofeir 


Teer 


eer g: 


etee 


stent 


pan 


*¢re 


faeol 


saaot 


saute 


sean 


40761 


mACHDAmMVEM=A 


Zé 


Rx 


290 


300 


320 


330 


3406 


350 


360 


370 


380 


390 


400 


410 


420 


830 


4406 


450 


860 


470 


4380 


ae 
sx 
RE 


NOTE: 


34561 


35027 


35484 


35934 


36381 


36825 


37268 


37710 


38152 


38595 


39039 


39484 


39930 


40378 


40827 


41273 


41732 


42186 


42643 


TABLE A&.3 
(ANGUS ET AL. 
VAPOUR 

3.0 3.5 
33775 33437 
34295 34009 
34791 34540 
35270 35047 
35740 35538 
36202 36018 
36660 36490 
37114 36957 
37566 37420 
38018 37881 
38469 38341 
38920 38800 
39371 39258 
39824 297 07 
40277 40176 
40732 40636 
41188 41097 
41645 4i1S59 
42104 42022 
82565 42487 


THE VALUES 


HAVE BEEN 


THE SATURATION TEMPERATURES 


OF H AT 290 K AND 5.5 MPA AND AT 290 K AND 6.06 MPA 
FROM THE APPROPRIATE VALUES AT 
<€ 300 K AT THESE PRESSURES) 


“BACK-EXTRAPOLATED" 
(290 < TCSAT) 


1973) 


34274 


34813 


35329 


35828 


36316 


36797 


37272 


37743 


38211 


38678 


39144 


39609 


80074 


405839 


41005 


41472 


41939 


42408 


32602 


33350 


33987 


34566 


35110 


35632 


36138 


36633 


37121 


37603 


38080 


38556 


393029 


39501 


33972 


80483 


40914 


41385 


418S7 


482330 


ENTHALPIES GIVEN IN J/MOL 


PRESSURE 

5.90 5.5 
32030 31258 
32955 324893 
33675 33330 
34304 34023 
34882 34643 
35429 35219 
35955 35767 
36466 36295 
36967 36810 
37460 37316 
37948 37814 
38432 38307 
38913 38796 
393382 393282 
33869 33766 
40346 40248 
40822 40730 
41298 41211 
41774 416392 
42251 42172 


STORED IN FILE CO2DATA 


IN MP4 


32982 


33721 


34390 


35000 


35573 


36120 


36651 


37169 


37679 


38181 


38678 


39172 


39662 


40151 


40637 


41123 


41603 


82094 


23734 


25069 


27909 


31598 


32848 


33748 


345064 


35181 


35809 


36405 


36978 


37535 


38078 


38612 


39138 


39658 


40174 


40686 


41195 


41702 


CARBON DIOXIDE ENTHALPY DATA FROM I.U.P.A.C. 


24975 


26392 


30930 


32489 


33458 


34267 


34379 


35632 


36246 


36834 


374802 


37956 


38498 


33032 


39559 


40081 


40598 


81112 


41624 


24896 


26620 


30126 


32032 


33152 


34022 


34772 


35451 


36085 


36688 


37269 


37833 


38385 


38926 


39460 


33988 


4osii 


41030 


41546 


28828 


26388 


29264 


31566 


32830 


33769 


34560 


38267 


35922 


36541 


37135 


37710 


382771 


38820 


393361 


39895 


80423 


40947 


41468 


94 


24426 
25510 
26716 
28090 
29604 
31058 
32296 
33334 
34232 
35036 
35777 
35471 
37131 
37765 
38379 
38378 
393563 
40139 


80705 
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Appendix 5 - Sample Experimental Data Files 


t, 

Bs 

win 
©.00 10.00 20.00 
1.201 1.208 
1.105 
1 Uae 1.179 


~ 
we 
e 
~ 
- 
w 
@ 
oe 
- 
uw 
e 
o 


-063 

1.38 1.99 1.98 
2.836 2.763 2.9895 
16003 1002. 1003. 
6605. 
703.2 

23. 

72.8 

6,1, 

©.00 10.060 26.0600 
1.217 1.214 
1.103 

1.181 

7.083 7.080 7.108 
1.460 
12.701 12.704 12.721 

-064 

2.46 2.47 2.46 
5.875 6.185 6.813 
1000 

6805s. 6810. 6505S. 
702.7 

23.5 

73. 

Ss, 


1.292 1.296 
1 045 12.002 12.038 
o73 
ibe ys 
3.031 3.1835 3.3575 
1387 i39s. 
$820. $$2s. 9510. 
702.4 
23.5 
73. 
CCeCccC N2 DATA JULY 21, 


DATA FILE N2.2 


€ KGAS ®& 


C KRUNS & WO. 


1 FOR W2, 


C mPTsS = wo. 


30.00 


1.180 
6.828 


1.385 
12.436 


1.98 
3.216 


6610. 


30.00 


23 T/C 


2 FOR CO2, 


3 FOR OTHER GASES 
OF EXPT’L RUNS IN DATA FILE 


OF DATA PTS IW RUN; 
40.00 $0.00 61.00 c 
1.208 1.208 
1.108 1.108 
1.182 1.183 e 
¢c 
6.852 6.5498 6.880 
1.383 1.388 1.388 € 
12.4853 12.470 12.888 € 
€ 
Ou c 
3.438 3.664 3.911 c 
1002. 
80.00 80.00 
1.210 
1.107 9.1032 
1.180 
7.104 7.187 
1.480 1.483 
12.738 12.785 
2.858 2.46 
6.9518 7.223 
6510. 6600 
40.00 
1.220 
1.188 
5.503 
12.000 
3.687 
1396. 
HTR (300 MA) 
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KTVP © 1 FOR WORMAL CALC’N 
TIME OF READINGS 


EXPERIMENTAL DATA : T/C, PRESSURE, FLOWMETER. 
AND OTHER READINGS GIVEN HERE 


BLAWK ENTRIES AT ANY TIME INDICATE THAT 
READING IS UNCHANGED FROM PREVIOUS VALUE 


TME TWO BLAWK LINES 16% BACK RUN DATA SET 
ARE AVAILABLE FOR STORAGE OF OTHER 
EXPERIMENTAL DATA, IF REQUIRED 


< 


eelit stat 


*' 3482 wee stats 
—caeionen 


Sorte ss  taweeens OAT |. ata seragmeongee ¥ 
¥eie wiers esureaka ave if 


Tae? Seeger Sart eee Fe apeerar oe 
OUISY CUCL VEER HER OF ERSWAED BF 


'¥2 «790 @4a O2eS Ot Qyeis wees owe 
aguyYO 3G BOEWETS 965 Prenat ata. 
ete(uewe St TAO satepe) oe 3 


2, 
3, C SEE EXPLANATORY COMMENTS ON PREVIOUS PAGE 
2,0, C KTYP = © HERE FOR "BACKGROUND" HEAT LEAK CALC’N 


089 
es alee 
1.238 1.240 
1.198 1.201 
8.850 852 
230 .129 
148 .148 
698.1 
22. 
72. 
Sei 
0.00 10.00 25.00 35.00 485.00 
1.180 TS 7s e 
1.083 


6.715 6.704 G50 7 hg, 6.709 6.714 
1.383 1.382 1.381 
12.707 12.657 12.639 12.604 12.589 


Z2nhe 2.05 25 ey 2.06- 2.06 
§.485 5.647 5.890 6.083 6.215 


494. ass. 494. 
3030. 3035. 3035. 3025. 3025. 
697.93 

I ee 

ano 

Lin Aa 

°o.00 12.00 20.00 40.00 56.00 
tr ts he 182 

TOF 1.076 1.078 
1.187 


6.106 6.086 6.067 5.330 5.881 


jae oh 1.315 1.313 1.302 tie? S47) 
12. 800. 12.332 12.8900 12.716 12. 746 
-068 
1.65 1.64 1 Giz 1252 1.50 
yee he) Tia AS he 380) ho S32 TS'3 
500. s0O1. $62. 
3210. 3220. 3230. 3280. 
687.7 
23. 
73. 


ccecec C02 DATA JUNE 30 
CCCCC DATA FILE €CO2.1 


wer ates 


218g 
a at a 


- 
a) 
- 
w 


= 
on 
a 
~ 


C SEE COMMENTS ON PREVIOUS PAGES 


11.00 20.00 30.00 40.00 $0.00 60.900 


1.182 UPS CFE} 


~ 
wo 
- 
@ 
—_ 
wi 
wv 
N 
= 


1.186 
1.091 
ct [4 


7.838 Tae Sine 7.999 8.000 
-S83 1.586 
.230 13.264 


.O73 074 
- 44 -44 
.299 2.5615 


13.198 1 176 13 
o70 O70 o71 O72 
33 -40 -41 - 42 
1.308 1.8265 ae 26 2.037 z 
1O0o. 1008. 1002. 1008. 1006. 
6665. 6690. 6645. 6665. 6680. 
700.9 700.8 
10.00 20.00 30.00 20.00 
Tr. 1950 
S61 &.872 8.893 8.616 
1.664 1,667 1.671 1.674 
13.376 13.3856 13.446 13.536 
078 O73 o8i 
53 54 $3 $4 
8.137 4.482 4.767 §.083 
6SiS§ 6510 
10.006 36.00 40.00 
1.178 1.186 
1.083 
ube ee Uy tS 
247 8.205 8.200 
1.641 1.606 
13.264 13.318 13.313 
101 -O83 O54 
47 a7 475 
6.835 7.405 7.687 
1005S. 1006 1005 
6620. 6625 
10.00 20.00 30.00 489.00 
lise UK 
1,093 
StS 1.182 
702 7.740 Gf Ae} 7.743 
hes 25 Pe Si2i7) 1.831 135:31 
13.188 13.1856 13.078 13.077 


11.00 20.00 30.00 40.00 


1.451 
12.840 12.861 12. 
O60 
34 345 : 


345 .34 


2.649 2.8480 3.0525 3.265 


MIX DATA AUG 30,31 
OATA FILE MIX. 2L 


(USING LMFM) 


Sus 


2¢e8. © @€af.5 
.8ae6 


+r 
ret 


‘eta. . 
2) Ta Te a 
269.8 YR a) ed i 


5 a ry 
sta 
ih) 


re. 


@o @8. 9a;8¢ 
par,¢ <8Pt,t 
coo + a 
avy. 


eos 4 i) 
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era ty PSever 
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avs so. 
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Appendix 6 - Experimental Results 


The following pages contain all of the results of the 
later experimental work. The results of earlier experimental 
runs have not been included because changes in equipment 
configuration and operation have caused at least some 
improvement in the later performance of the calorimeter over 
that of earlier work. 

The results as shown are in a form Similar to the 
actual output of the program CALOR. Minor additional notes 


are given for explanatory purposes. 
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aiicsen £ 


eft to ativess eft Io tie ach 
Isinemiteqxe tsiliuse To asin edt wa ow 
tnemgiups at eepnedty oeuased nisi 
emoe jeeei #5 beeyves ovat nolaarige 8 ie 


aa 
ssvo tetemizoles ef? to agammiotiveq resal wats 


430 Val : te 


ay ; “ii 
eit of talimie m0? 5 Ai ste fnwode <a wast 8 


aa 


seton Lenottibbe zonié sehaiies antyouy eey 305: 
_aeadques riSvanaiqne on 


DATA FROM FILE N2 


INTERVALS MARKED WITH ** WERE 


~1 


MOST STABLE 


RUN 1 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL ) 
TEMPERATURE Lower Iniet 370.98 
«K) Outiet 298.50 
PRESSURE Iniet 9.70601 
(MPa) Outiet 8.611 
Atmospheric ©.094 
HEAT EXCHANGE (J/min) 39.128 
HEAT LEAK (J/min) 5.082 
RATIO: LEAK OVER EXCHANGE OF. P29 
ENTHALPY CHANGE Predicted 2:3 4010. 26 
(J/mol) Measured °2688.393 
FLOW RATE Lea ha | Pred. Mo}. 1.37048 
(i/min) of Meas. Mol. 1.57368 
(mmol /min) Caic. Vo}. ©.03970 
N2 Meas. Voi. 0.38227 
Caic. Mo). 14.55428 
Freon Sat’n Temp. (K) 294.73 
Freon Meter Temp. (K) 296.83 
Freon Return Temp. (K) 294.74 
Heat Exch. Corr. (J/min) °0.002 
Ambient Air Temp. (K) 295.65 
RUN 2 TEST GAS - N2 
CONBITIGNS FOR TEST INTERVAL i 
TEMPERATURE Lower Iniet 401.21 
(K) Ouzriet 301.76 
PRESSURE Inlet 8.739 
(MPa) Outies $.499 
Atmospheric 0.094 
HEAT EXCHANGE (J/min) 78.229 
HEAT LEAK (J/min) -10.788 
RATIO: LEAK OVER EXCHANGE 0). 13.7.9 
ENTHALPY CHANGE Predicted 318i). 6:7 
(J/mol) Measured °2796.09 
FLOW RATE Festi Pred. Mol. 3.580198 
(l/min) er Meas. Moi. 3.14631 
(ramol/min) Caic. Voi. ©.07838 
NZ Meas. Vo). 0.73624 
Caic. Mol. 27.837809 
Freon Sat‘n Temp. (K) 294.73 
Freon Meter Temp. (K) 297.48 
Freon Return Temp. (K) 294.74 
Heat Exch. Corr. (Jd/min) 9.003 
Ambient Air Temp. (K) 296.21 
RUN 3 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL i 
TEMPERATURE Lower Iniet 411.14 
(K) Outlet 302.61 
PRESSURE Iniet 3.743 
(MPa) Outiet 9.384 
- Atmospheric 0.094 
HEAT EXCHANGE (J/min) 104.904 
HEAT LEAK (J/min) 18.447 
RATIC: LEAK OVER EXCHANGE On att 
ENTHALPY CHANGE Predicted -3457.49 
(J3/mo1) Measured -3038.98 
FLOW RATE F113 Pred. Mol. 4.80028 
(i/min) or Meas. Moi. 4.21923 
(mmol/min) Cale. Vo). ©.10642 
N2 Meas. Vol. ©o©.3930887 
Calc. Mo}. 34.51961 
Freon Sat’n Temp. (K) 294.74 
Freon Meter Temp. (K) 297.88 
Freon Return Temp. (K) 294.73 
Heat Exch. Corr. (J/min) ©.005 
Ambient Air Temp. (K>) 296.48 


3 
3 
° 
° 
28 


.57268 
- 10222 
-07827 
»7&332 
. 24706 


294.73 
237.60 
294.75 
-0,.007 
296.21 


-15 


IN THE RUW 


sz 
3 


369.21 
238.35 
8.694 
$.609 
0.094 


37.666 
4.821 
©.1280 
°2290.62 
-2626 .85 


1.32098 
1.51489 
0.03822 


0.37661 
14.33878 


294.73 
296.92 
284.72 

0.002 
238.65 


76. 
ah a 
0). 

21/114. S'3 
2 72ar 19 


-$3833 
-094687 
-07808 


OuwwW 


©.74332 
- 24706 


294.73 
297.76 
294.75 
“0.010 
296.21 


3 


812.02 
302.398 
3.7683 
9.394 
0.094 


106.732 
“13.165 
-0.1233 
-3473.28 
-3091.91 


4.822193 
4.29272 
0.10828 


0.908397 
-51961 


rig 
2 


1 
1 
° 


c 
14 


1 


rj 
-3 


4. 
4. 
°. 


3. 
34). 


zs 
4 


369.07 
288.31 
§.691 
8.606 
0.094 


37.300 

4.795 
©.1285 
287.30 
624.70 


- 30734 
.50019 
-O3785 


-37327 
-21128 


294.73 
286.99 
294.73 

©.000 
235.65 


as 
4 


412.25 
302.98 
$.760 
9.399 
0.094 


Oi 
12.681 
©.1180 
480.43 


113.08 


83213 
32211 
10902 


90897 
$1961 


462 


368.35 
2938.28 
2.6391 
9.617 
©.094 


36.208 
4.717 
©.1303 
°-2266.01 
°2608.45 


1.266806 
1.45610 
6.03674 


0.36497 
.- 89554 


294.73 
297.06 
294.74 
“0.001 
295.65 


ss 


5 


at2.39 
302.97 
3.760 
9.394 
©.094 


106.732 
°-12.633 
“0.1184 
-3485.03 
-3116.19 


4.80082 
4.28272 
°o.10828 


©.90189 
.- 25070 


367.389 
298.26 
9.691 
8.611 
0.094 


36.204 
4.8858 
6.1343 
°2251.94 
°2603.21 


1.259671 
1.45610 
0.03674 


0.366293 
- 90746 


294.73 
297.12 
294.74 
-0.001 
295.65 


ze 


401.36 
301.80 
8.739 
$.494 
0.094 


78.595 
°11.536 
-6.1468 

°3184.72 
= 2777. 32 


3.62496 
3.16101 
0.07875 


0.74873 
. 390083 


294.73 
eee 
284.75 
°-0.011 
296.21 
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RUN 4 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 380.79 
(K ) Outlet 298.99 
PRESSURE Iniet 6.981 
(MPa) Outiet 6.781 
Atmospheric ©.9093 
HEAT EXCHANGE (J/min) 54.481 
HEAT LEAK (J/min) 7.205 
RATIO: LEAK OVER EXCHANGE ©.1322 
ENTHALPY CHANGE Predicted °-2563.06 
(J/mol) Measured °-2953.68 
FLOW RATE anit Pred. Mei. 1.90122 
Ci/min) or ,Mea@s. Mol. 2.19097 
(mmo 1 /min) Caic. Vol. 0.05536 
N2 Meas. Voi. 0.48705 
Calc. Mol. 18.488502 
Freon Sat’n Temp. (K) 284.68 
Freon Meter Temp. (K) 237.39 
Freon Return Temp. (kK) 294.60 
Heat Exch. Corr. (J/min) ©,.0623 
Amoient Air Temp. (K) 286.76 
RUN 5 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 3936.90 
(K) Outiet 300.59 
PRESSURE Iniet 6.3985 
(MPa) Outlet 6.694 
Atmospheric 0.094 
HEAT EXCHANGE (J/min) 55.200 
REAT LEAK (J/min) °-16.8393 
RATIO: LEAK OVER EXCHANGE -0.1968 
ENTHALPY CHANGE Predicted -2819.86 
(3/moi) Measured =-2357.03 
FLOW RATE Feit Pred. Mol. 2.65635 
(l/min) or Meas. Moy. Za 20ST 
(mmol /min) Caic. Voi. 6.05591 
NZ Meas. Vol. 0.61447 
Caic. Moi. 23.41945 
Freen Sat‘n Temp. (K) 2934.78 
Freon Meter Temp. (Kk) 296.98 
Freon Return Temp. (K) 294.81 
Heat Exch. Corr. (J/min) -0.608 
Ambient Air Temp. (K) 285.93 
RUN 6 TEST GAS - NZ 
CONDITIONS FOR TEST INTERVAL i 
TEMPERATURE Lower Initet 403.25 
(K) Outiet 302.33 
PRESSURE Iniet 6.988 
(MPa) Outiet 6.594 
Atmospheric °o.094 
HEAT EXCHANGE (J/min) 70.915 
HEAT LEAK (J3/min) 32 tae ue 
RATIO: LEAK GVER EXCHANGE -0.3031 
ENTHALPY CHANGE Predicted -3141.24 
(J/mo1) Measured -2410.53 
FLOW RATE P= Pred. Mol. Se Weil aS fe Bo 9] 
(l/min) or Meas. Moi. 2.85236 
(mmol /min) Calc. Vol. 0.07187 
NZ Meas. Vol. 0.77305 
Calc. Mol. 29 .41876 
Freon Sat‘n Temp. (K) 294.77 
Freon Meter Temp. (K) 297.05 
Freon Return Temp. (K) 294.79 
Heat Excn. Corr. (J/min) -@.009 
Ambient Air Temp. (K) 296.21 


2807.02 


be i 3 le 
-2428.31 


z. 
2 


° 


0.76739 
29 .20337 


294.77 
297.05 
2834.73 
-0.009 
296.21 


Or} 
OF3 


1 
2. 
° 


x 
3 


3792.29 
298.88 

6.978 
6.783 
.093 
§1.182 

$.783 
0.1130 
$20.46 
841.85 


-82613 
OS869 
-08202 


-47872 
.01589 


294.68 
297.45 
294.60 

0.021 
296.76 


° 


e 
22 


-OS5S1 


- 60032 
- 87984 


294.78 
296.39 
294.81 
-6.008 
235.33 


2 
°2 


1 
2. 
° 


}. 


V7 


379.46 
288.80 
6.374 
6.783 
6.093 


§1.%382 

6.186 
0.1208 
$28.18 
875.69 


» 80982 
OS869 
05202 


47006 
- 80148 


294.68 
297.53 
294.60 

©.021 
296.76 


e4 
“2 


1 
725 
te] 


100 


378.88 
2838.78 

6.974 
6.783 
©o.033 


§1.587 

6.530 
0.1267 
$23.39 
896.21 


.-81678 
O7339 
-05239 


- 47006 
- 80154 


294.68 
237.49 
294.60 

©.021 
296.76 
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RUN 9 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outiet 
PRESSURE Iniet 
(MPa) Outiet 


Atmospheric 


HEAT EXCHANGE 
HEAT LEAK 


(J/min) 
(J/min) 


RATIO: LEAK OVER EXCHANGE 
ERNTHALPY CHANGE Predicted 
(J3/mo)) Measured 
FLOW RATE Foie Pred. Mol. 

(T/min) or Meas. Mol. 
(mmol /min) Calc. Vol 
N2 Meas. Vol 
Caic. Mol 
Freen Sat’n Temp. (K) 
Freon Meter Temp. (K) 


Freon Return Temp. (K) 
Heat Excn. Corr. (J/min) 
Ambient Air Temp. (K) 


RUN 8 TEST GAS - 
CONDITIONS FOR TEST 


N2 
INTERVAL 


TEMPERATURE 
(K) 

PRESSURE 
(MPa) 


Lower Iniet 
Outies 
Iniet 
Outlet 


Atmospheric 


HEAT EXCHANGE (J/min) 

HEAT LEAK (JU/min) 

RATIG: LEAK OVER EXCHANGE 

ENTHALPY CHANGE Predicted 
(J3/moi) Measured 


FLOW RATE Fei 
(l/min) or 
(mmo1/min) 


Pred. 
Meas. 
Caic. 


Mol. 
Mol. 
Vol. 
N2 Meas. 
Caitc. 
Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (Kk) 
Heat Excn. Corr. (J/min) 
Ambient Air Temp. (K) 


N2 
INTERVAL 


RUN E) TEST GAS - 
CONDITIONS FOR TEST 
TEMPERATURE Inlet 
(K) 
PRESSURE 
(MPa) 


Lower 
Outiet 
Iniet 
Outiet 
Atmospheric 


HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 


RATIO: LEAK OVER EXCHANGE 
ENTHALPY CHANGE Predicted 
(J/mo1) Measured 
FLOW RATE Lo Fi | Pred. Mol. 

(l/min) or Meas. Mo}. 
(mmol/min) Caic. Vol. 
N2 Meas. Yo!. 
Calc. Mol. 
Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 


1 2 
361.30 361.26 
ZOE ls 297.76 

3.520 Sie Sti: 

3.309 3.306 

0.098 o.094 
31.085 31.085 

8.872 9.102 
0.2854 0.2928 

= WBZ. S26 =) 119i 20r. 9.0 
-2690.90 -2716.28 
©.88349 ©.88421 
1.28033 1.25033 
©.03150 ©6.03150 
0.30299 ©.30016 
11.55204 11.44814 
294.77 294.77 
296.64 296.66 
2G tel, 294.77 
-0.001 -0,001 
235.65 2935.65 
xz 

1 2 
371.82 37k 39'S 
298.47 2398.50 

3. 527 2 S207, 

3.194 3.196 

0.094 0.094 
35.470 35.470 

1.895 1.424 
©.0534 0.0401 

2 Sin - er ~eieila. Sia 
-2321.83 -23058.54 
1.35047 1.36942 
1.82670 1.42670 
©.03535 ©.03595 
0.80068 ©.40351 
15.27680 15.38481 
294.77 294.77 
296.86 296.91 
294.898 294.89 
-0.022 0.022 
295.65 2985.65 
zz 

1 2 
380.68 380.55 
299.32 299.34 

3.541 3.541 

3.094 3.094 

0.094 0.094 
43.510 43.875 
75 TeX) i Are) 
-©.0621 -0.0510 

-2444.11 -2439.67 
22301 30) =23 24). 2:3 
1.85865 1.85479 
1.75005 1.76474 
0.04411 0.04448 
©.49696 0.49683 
18.90669 18.90178 
294.76 294.76 

297.24 239:7.723 

Freon Return Temp. (K) 294.83 294.83 
Heat Exch. Corr. (J/min) “0.016 -0.016 
Ambient Air Temp. (K) 296.21 236.21 


-2 


1 
1 


°. 


3°. 
WS 


343.58 


~ 3533 
-43405 
03613 


393901 
21303 


294.77 
296.93 
294.89 
-0,.022 
295.65 


Ors 


401.21 


-891593 
-118085 
.02817 


- 30362 
.-57615 


294.77 
296.72 
294.77 
-0.001 
295.65 


Pr: 


LOE 


428.91 


0.88205 


- 11805 
02817 


. 30016 
- 844818 


294.77 
296.78 
298 17, 
-0.001 
295.65 
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DATA FROM FILE N2.2 


INTERVALS MARKED WITH ** WERE MOST STABLE IN RUN 


RUN 1 TEST GAS - N2 ze 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 335.30 
(K) Outiet 300.34 
PRESSURE Inlet 7.006 
(MPa) Outiet 6.699 
Atmospheric 6.094 
HEAT EXCHANGE (J/min) 73.836 
HEAT LEAK (J/min) Cee fh 
RATIO: LEAK OVER EXCHANGE ©.c0Si 
ENTHALPY CHANGE Predicted 2982.51 
(J/mol) Measured °-2987.83 
FLOW RATE Feit Pred. Mol. 2.98476 
(l/min) or ‘Meas. Mol. 2.96994 
(mmo l/min) Cate. Vol. 0.07473 
N2 Meas. Vol. 0.64562 
Caic. Mol. 24,62967 
Freon Set’n Temp. (K) 294.78 
Freon Meter Temp. (K) 296.70 
Freon Return Temp. (K) 298.91 
Heat Exch. Corr. (J/min) °0.046 
Ambient Air Temp. (K) 2935.65 
RUN z TEST GAS - NZ2 
CONDITIONS FOR TEST INTERVAL i 
TEMPERATURE Lower Iniez? 408.36 
(K) Outtet 301.53 
PRESSURE Iniet 6.3988 
(MPa) Outtet 6.601 
Atmospheric ©.084 
HEAT EXCHANGE (J/min) $1.378 
HEAT LEAK (J/min) Oy ag ee! 
RATIO: LEAK OVER EXCHANGE -0.0298 
ENTHALPY CHANGE Predicted * 3231.63 
(J/mol) Measured °31338.13 
FLOW RATE Fgentct Pred. Mol. 3.78494 
(l/min) or Meas. Mol. 3.67842 
(ramoi/min) Catc. Vol. 6.08262 
N2 Meas Vol. 6.76455 
Calc. Mol. 29.11873 
Freon Sat’n Temp. (K) 294.76 
Freon Meter Temp. (K) 297.01 
Freon Return Temp. (K) 294.87 
Heat Execn. Corr. (J/min) °0.049 
Ambient Air Temp. (K) 295.93 
RUN 3 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 376.12 
(K) Outiet 298.48 
PRESSURE Iniet? 9.726 
(MPa) Outiet 9.616 
Atmospheric 0.094 
HEAT EXCHANGE (J/min) 47.530 
HEAT LEAK (J/min) 6.364 
RATIO: LEAK OVER EXCHANGE 0.13398 
ENTHALPY CHANGE Predicted -28504.06 
(J/mol) Measured -2891.20 
FLOW RATE Lot it Pred. Mol. 1.68574 
(l/min) or Meas. Mol. 1... G:.A,72 
(mmo l/min) Calc. Vol. °o.04818 
N2 Meas. Vol. 9.43183 
Caic. Moi. 16.43965 
Freon Sat’n Temp. (K) 294.75 
Freon Meter Temp. (K) 296.393 
Freon Return Temp. (K) 294.33 
Heat Exch. Corr. (J/min) -0.042 
Ambient Air Temp. (K) 295.93 


2 


395.94 
300.35 
7.006 
6.699 
0.094 


73.470 
©.475 
0.0065 
-2983.34 
-3002.74 


2.93616 
2.985524 
0.07442 
0.64138 
24.46773 


234.78 
296.73 
294.91 
-0.046 
295.65 


-3148.91 


3.76045 
3.66072 
©o.08225 
o.758898 
28.902395 


294.76 
296.98 
294.87 
-0.049 
295.93 


zs 


3 


405.74 
301.54 
6.988 
6.599 
0.0394 


90.647 
-3.606 
"9.0398 
°3242.893 
-3118.81 


3.79168 
3.64603 
0.09188 


©.76314 
06473 


294.76 
296.35 
294.87 
“0.048 
295.83 


=z 


3 


376.85 
298.56 


ze 
5 


395.63 
300.25 
7.009 
6.704 
0.094 


73.105 
©.74&7 
0.0102 
-2976.89 
-3007.74 


-91088 
- 94085 
-974805 


OoOnwW 


©.63713 
24.30557 


294.78 
296.858 
294.97 
°-0.067 
2958.65 


3 


806.07 
301.56 
6.388 
6.5998 
0.098 


91.013 
-6.216 
-0.0463 
-3282.29 
-3108.32 


3.830268 
3.66072 
©.08225 


6.76886 
29 .28080 


294.76 
296.91 
294.391 
0.066 
295.93 
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zs 
6 


395.61 
300.27 
7.006 
6.704 
©.094 


73.105 
©.3920 
0.0126 
-2878.90 
-3013.81 


- 90356 
-940S5 
.07405 


ONWN 


0.63584 
- 25661 


294.78 
296.86 
294.93 
°-0.083 
295.65 
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DATA FROM FILE N2.3 


INTERVALS MARKED WITH ** WERE MOST STABLE IN RUN 


RUN 1 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(kK) Outiet 
PRESSURE Inlet 
(MPa) Out iet 


Atmospheric 


HEAT EXCHANGE (J/min) 

HEAT LEAK (J/min) 

RATIO: LEAK OVER EXCHANGE 

ENTHALPY CHANGE Predicted 
(J/mol) Measured 


FLOW RATE (oy ten Pred. Mo}. 
(l/min) or ‘Meas. Mol. 
(mmot/min) Caic. Vol. 


ux 
1 


338.77 


297.88 


6.8383 
6.794 
@.0984 


21.876 
6.318 


@.2928 
°1302.75 
-1882.24 


° 
° 
© 


-61395 
- 86819 
-O2173 


2 


338.73 


oO. 


° 
°o 


297.59 
6.833 
6.794 
0.094 


21.981 


-%1301.29 
°1873.45 


61325 
- 88289 
02210 


N2 Meas. Voi. 
Calc. Mol. 


Freon Sat’n Temp. ( 
Freon Meter Temp. ( 
Freon Return Temp. (K) 
Heat Exch. Corr. (J/min) 
Ambient Air Temp. (K) 


RUR 2 TEST GAS - W2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Cutie? 
PRESSURE Iniet 
(MPa) Outlet 


Atmospheric 


HEAT EXCHANGE (J/min) 

HEAT LEAK (J/min) 

RATIO: LEAK OVER EXCHANGE 
ENTHALPY CHANGE Predicted 


(Jd/moei) Measured 
FLOW RATE Feil Pred. Mol. 
(l/min) oF Meas. Moi. 
(mmol /min) Catc. Voi. 
N2 Meas. Voi. 


Caic. Mol. 


Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (K) 
Heat Exch. Corr. (Jd/min) 
Ambient Air Temp. (K) 


RUN 3 TEST GAS - N2 
CONDITIONS FOR TEST INTERVAL 


TEMPERATURE Lower Iniet 
(K) Outlet 
PRESSURE Iniet 
(MPa) Outiet 
Atmospheric 


HEAT EXCHANGE (J/min) 

HEAT LEAK (J/min) 

RATIO: LEAK OVER EXCHANGE 

ENTHALPY CHANGE Predicted 
(J/mot) Measured 


FLOW RATE Fr qt Pred. Mol. 
(Ci/min) or Meas. Mol. 
(mmol/min) Calc. Voi. 


° 
11 


oat 
=a 


1 
1 
° 


o 
20 


irs 
-1 


- 30882 
71176 


i 


359.40 
2939.06 
6.837 
6.604 
0.094 


36.552 
-1.666 
0.0456 
892.67 
810.18 


-$3782 
.487078 
-O3682 


Oe 
- 19246 


298.96 
298.11 
294.89 

0.012 
296.65 


300.72 
6.968 
6.634 
0.094 


51.894 
10.037 
0.1934 
317.26 
941.71 


° 


Tt. 
lig 
©}. 


aa) 
= 


- 30582 


¥1.71166 


294.97 
2c Orr 
294.97 
-0.000 
296.48 


ae 
2 


359.07 
2993.05 
6.837 
6.604 
©.094 


36.552 
°0.848 
0.0232 


-1882.77 
-1840.10 


=e 
2 


374.87 
300.69 
6.968 
6.634 
0.094 


51.894 
10.100 
©.1946 
317.04 
939.57 


-49447 
-08809 
.05229 


N2 Meas. Vol. 
Calc. Mol. 


Freon Sat‘n Temp. (K) 
Freom Meter Temp. (K) 
Freen Return Temp. (K) 
Heat Exch. Corr. (J/min) 


Ambient Air Temp. (K) 


-6€9942 


°. 


70020 


-72606 26.75560 


294.95 
298.48 
294.37 
-0.006 
296.93 


294.95 
298.50 
294.97 
-0.006 
296.93 
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RUN 4 TEST GAS - N2 
CONDITIONS FOR TEST INT 


TEMPERATURE Lower 
(K > Outier 

PRESSURE Inte 
(MPa) Outlier 
Atmosp 


HEAT EXCHANGE (J/min?} 
HEAT LEAK (J/min) 


RATIO: LEAK OVER EXCRAN 
ENTHALPY CHANGE Predic 
(J/moi) Measur 
FLOW RATE Fai Pred. 
(l/min) or ‘Meas. 
(mmol /min) Canc. 
N2 Meas. 
Calc. 
Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (K) 
Heat Excn. Corr. (J/min 


Ambient Air Temp. (K) 


RUR 5 TEST GAS - N2 
CONDITIONS FOR TEST INT 


TEMPERATURE Lower 
(K) Outiet 

PRESSURE Iniet 
(MPa) Outiet 
Atmesp 


HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 
RATIO: LEAK OVER EXCHAN 


ENTHALPY CHANGE Predic 
(J/mot) Measur 
FLOW RATE oe Vi | Pred. 
(i/min) or Meas. 
(mmol /min) Calc. 
N2 Meas. 
Cale 


Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (K 
Heat Exch. Corr. (J/mi 
Ambient Air Temp. (K) 


RUN 6 TEST GAS - N2 
CONDITIONS FOR TEST INT 


ERVAL 1 2 


Intet 395.47 395.32 


301.32 301.28 
6.389 6.9893 
6.6884 6.684 
heric 0.084 0.094 


§8,.200 58.200 
°-12.789 -°41.866 


GE -0.20593 -0.1837 
ted -2937.62 -2834.62 
ed °2836.04 -2857.95 
Mol. 2.87249 2.84341 
Mol. 2.38208 2.38204 
Vol. ©.0S8966 ©.OS866 
Voi. ©.63571 0.63005 
Mol. 24.30177 24.08519 


294.34 294.94 
298.18 238.19 
294.87 298.87 
) 0.021 0©.621 
296.76 296.76 


=f 


ERVAL 1 2 


Iniet 417.83 417.84 


304.78 304.84 

6.996 6.996 

6.387 6.384 

heric ©.094 0.0984 


90.614 91.344 
-36.317 -36.284 


GE -0, 4008 °0.3972 
ted -3800.27 -3498.62 
ed -2488.78 -2503.398 
Mol. §.10732 5.13839 
Mol. 3.64603 3.67542 
Vel. 0.08132 0.09206 
Vol. 0.94861 &©.95428 
Mol. 36.26314 36.473961 


294.94 294.394 
298.32 288.37 
294.87 294.87 
) 0.032 0.032 
296.76 296.76 


ERVAL 1 2 


TEMPERATURE Lower Iniet 394.37 394.76 
(K) Outlet 303.43 303.48 
PRESSURE Iniet 6.986 6.996 
(MPa) Outiet 6.394 6.384 


Atmospheric 0.094 °o.094 


HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 


RATIO: LEAK OVER EXCHANGE 
ENTHALPY CHANGE Predicted 


76.001 76.732 
-26.839 -26.566 
°-0.3545 -0.3462 

-2822.15 -2831.96 


(J/mot) Measured -2083.61 °-2103.64 
FLOW RATE | ia I | Pred. Mol. 4.14209 8.15648 
(i/min) or Meas. Mol. 3.05813 3.08752 
(mmol/min) Calc. Vol. ©.07658 0.07731 
N2 Meas. Vol. ©o.9S5S869 ©o.35569 

Calc. Mol. 36.87586 36 47577 

Freon Sat‘n Temp. (kK) 294.95 294.95 
Freom Meter Temp. (K) 298.96 2939.01 
Freon Return Temp. (K) 294.89 294.89 
Heat Exch. Corr. (J/min) ©.022 0.022 


Ambient Air Temp. (K) 


237.32 297.32 


=2 
2 


“3 
2 
° 


3 


394.62 
301.21 
6.989 
6.694 
0.094 


$8 .470 
11.243 
0.1823 
SS aed 
445.16 


. 86504 
35264 
.05893 


-62854 
.812388 


294.94 
298.22 
294.87 

6.920 
296.76 
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DAT& FROM FILE €62.1 


INTERVALS MARKED WITH *3 WERE MOST STABLE 


RUN 1 BACKGROUND TEST (NO FLOW) 
CONODITIGNS FOR TEST JNTERYVAL i 
TEST REAT LEAK (J/min) 6.533 
TEMPERATURE Lower Iniet 297.44 
(kK) OQutiet 2936.70 
Atmospheric Pressure (MPa) 9.033 
Freon Sat’nm Femp. (K) 294.58 
Freon Meter Temp. (Kk) 296.25 
Freon Return Temp. (K) 284.66 
Heat Exch. Corr. (J/min) 0.001 
Ambient Air Temp. (K) 2935). 30 
RUN 2 TEST GAS - E€602 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 398.52 
(kK) Outiet 300.23 
PRESSURE Iniet 3.502 
(MPa) Outlier? 3.126 
Atmospheric 6.093 
HEAT EXCHANGE (J/min) 76.061 
HEAT LEAK (J/min) -2.405 
RATIO: LEAK OVER EXCHANGE °0.0316 
ENTHALPY CHANGE Predicted -8496 .53 
(J/moi) Measured °4358.71 
FLOW RATE POU Pred. Mol. 3.15482 
(i/min) or Meas. Moi. 3.08813 
(mmo i/min) Calc. Vol. ©.07756 
Co2 Meas. Voi. ©.45873 
Caic. Mol. 17.886026 
Freon Sat’n Temp. (K) 294.857 
Freon Meter Temp. (KK) 296.29 
Freon Return Temp. (K) 294.48 
Heat Excn. Corr. (J/min) ©.034 
Ambient Air Temp. (K) 295.65 
RUN 3 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Initet 387.36 
(K) Outiet 298.35 
PRESSURE Iniet 3.540 
(MPa) Outiet 3.303 
Atmospheric 0.093 
HEAT EXCHANGE (J/min) 61.074 
HEAT LEAK (J/min) 3.096 
RATIO: LEAK GVER EXCHANGE 6.0807 
ENTHALPY CHANGE Predicted -8145 .41 
(J/moet) Measured -8366.77 
FLOW RATE LEM Pred. Mol. 2.33108 
(l/min) or Meas. Mol. 2.45553 
(mmoi/min) Calc. Vol. ©.06229 
co2 Meas. Vol. 0.36812 
Caic. Mol. 13.98611 
Freon Sat‘n Temp. (K) 294.56 
Freon Meter Temp. (K) 296.25 
Freon Return Temp. (K) 294.36 
Heat Excn. Corr. (J/min) ©.060 


Ambient Air Temp. (K) 2985.93 


-4@ 
-& 


3 
3 
° 
° 
les 


IN RUN 
. ze 
2 3 
388.54 3938.58 
300.22 300.21 
3.502 3.439 
3.128 3.123 
©.03s3 6.033 
76.732 78.426 
hoe by 473 °2.858i1 
©.60224 °-0.0338 
499.29 -44983.84 
800.63 -8352.8i 
- 18674 3.17662 
-08752 3.07282 
07831 0.867793 
-45873 0.46156 
45018 17.58791 
294.87 294.57 
296.29 296.35 
294.48 298.48 
6.9034 0.034 
295.65 298.65 
ax 
2 3 
387.01 385.59 
298.91 298.78 
3.544 3.547 
3.308 3.318 
©.093 0©.03s3 
60.343 $6 .687 
2.536 1.385 
0.04206 0.0244 
133.20 -4079 .64 
314.48 -4181.80 
32419 2.22348 
42613 2.27916 
06155 0.08782 
36812 0.35679 
$8611 13.55575 
294.56 294.56 
296.293 296.32 
294.36 294.35 
o.0s9 ©.058 
295.33 295.93 
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DATA FROM PILE CO2.2 


INTERVALS MARKED WITH ** WERE MOST STABLE 


RUN 1 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outiet 
PRESSURE Iniet 
(MPa) Outiet 
Atmospheric 
HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 


RATIO: LEAK OVER EXCHANGE 
ENTHALPY CHANGE Predicted 
(Jd/mo) Measured 
FLOW RATE | petty 5] Pred. Mol. 

(l/min) or ‘Meas. Mol. 
(mmol /min) Calc. Vor. 
C02 Mees. Voi. 
Calc. Mol. 
Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 


Freon Return Temp. (K) 
Heat Excn. Corr. (J/min) 
Ambient Air Temp. (kK) 


RUN 2 TEST GAS - C02 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outiet 
PRESSURE Iniet 
(MPa) Outiet 
Atmospheric 
HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 


RATIO: LEAK OVER EXCHANGE 
ENTHALPY CHANGE Predicted 
(J/mo}) Measured 
FLOW RATE Fost Pred. Mol. 

(i/min) or Meas. Moi. 
(mmol /min) Calc. Vol. 
co2 Meas. Voi. 
Calc. Mol. 
Freon Sat‘n Temp. (Kk) 
Freen Meter Temp. (K) 


Freon Return Temp. (K) 
Heat Exch. Corr. (J/min) 
Ambient Air Temp. (K) 


RUN 3 TEST GAS - C62 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Iniet 
(K) 
PRESSURE 
(MPa) 


Lower 
Outiet 
Iniet 
Outiet 
Atmospheric 


HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 


IN RUN 

1 2 3 
406.08 405.73 40S .80 
301.80 301.73 301.77 
2... 520 3.530 3.530 
3.098 3.108 3.0393 
©.033 0.083 0.033 
82.628 81.531 80.800 
7.343 “10.814 On i S51 fie ak} 
-6.0889 °-0.1326 °O.1507 
-8746.54 -4734.88 °-4727.88 
°&4357.32 °4180.40 -4108.77 
3.674795 3.713458 cS ae fee hi Aare 
3.32268 3.27859 3.249198 
©.08406 0.08294 0.08220 
0.49696 ©.51112 Ors Ts3 7 
18.962987 19.80323 19.66533 
288.65 284.65 294.65 
296.64 286.64 296.76 
294.70 294.70 2934.70 
“0.013 °-0.018 0.0198 
295.884 298.54 295.54 

a2 

1 2 3 
407.20 407.26 407.74 
301.73 361.78 301.80 
Soa. 3.520 KER I 
3.091 3.093 3.083 
9.0s3 0.0s3 0.093 
81.895 82.260 82.626 
-160.896 -10.174 °$.819 
°0.1282 “O77 1237) -6.11406 
°-4798&.51 °&78986.78 °4817.06 
-4249.858 -4268.82 -43246.12 
3.71T835 hoe elie 2a hi} 3.760146 
3.293329 3.30798 3.32268 
°o.08328 ©.08365 ©.684802 
0.50546 ©.50546 0.50121 
18.27007 19.27007 19.108i11 
284.66 294.66 2948.66 
rh Reo as 297.26 297.28 
284.70 284.76 294.70 
-0.015 -0.015 °0.015 
298.93 298.33 2938.83 

at 

i 2 3 
384.83 385.64 385.15 
299.41 299.48 2393.53 
ee 3 | 3.530 3.530 
Sy Se 3.318 3h a2 
©.03933 °0.0s3 ©.093 
§1.325 §3.621 §3.782 
OT 4.01 “2.044 =O. 131 
RATIO: LEAK OVER EXCHANGE -0,.0091 -0.0385 -0.0024 
ENTRALPY CHANGE Predicted 4002.88 -4048 43 +4026 .02 
(J/mol) Measured -3966 .83 -3898.19 -4016.24 
FLOW RATE Fett Pred. Mol. 2.10703 2.21436 2.16684 
(i/min) or Meas. Mo). 2.088093 2513218 2.16158 
(mmol/min) Care. VO 0.05280 0.08392 0.05466 
co2 Meas. Vol. 0.34366 0.35711 0.38139 
Calc. Mol. 13.08959 13.60158 13.38377 
Freon Sat‘n Temp. (K) 294.66 294.66 294.66 
Freon Meter Temp. (K) 296.94 29:25:42: 237 223 
Freon Return Temp. (K) 294.70 294.70 294.70 
Heat Exch. Corr. (J/min) -0.009 -0.009 -02.010 
Ambient Air Temp. (K) 2936.21 296.21 296.21 


<8 
°4 


Oww 


8 
-4 


Ow WwW 


12) 
ig 
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403.94 
301.68 
3.520 
3.103 
0.083 


79.338 
°-8.683 
0.1094 
654.98 
195.80 


. $3955 
. 19040 
- 98071 


49554 
90838 


294.65 
296.96 
294.70 
-0.018 
295.584 


& 


408.21 
301.86 
3.523 
3.086 
©@.0933 


83.357 
11.082 
©.1326 
831.53 
268.88 


- 79681 
.- 35208 
-O8 477 


-§12583 
- 53990 


284.66 
237.32 
234.70 
°0.018 
298.33 


-4 
°4 


Owuw 


5 


403.07 
301.57 
3.813 
3.108 
0.093 


SLOSS fea] 
°-93.584 
0.31231 
625.32 
118.49 


-51699 
-13161 
.07922 


- 49554 
- 90898 


284.65 
287.18 
294.70 
0.018 
295.54 


°-3 
-4 


i 
2 
° 


3.520 
3.321 
.@33 


-$258 

3.136 
0.0608 
$68.05 
223.12 


.96197 
o8s8os 
05280 


32281 
29541 


284.66 
287.32 
294.70 
-0.008 
296.21 
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° 4680.23 
-4167.51 


3.59941 
3.20510 
-08108 


-50121 
- 12506 


294.65 
287.28 
294.70 
-0.018 
235.54 
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RUN 6&8 TEST GAS - CO2 
COMDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outiet 
PRESSURE Iniet 
(MPa) Outiet 


Atmospheric 


HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 


RATIO: LEAK OVER EXCHANGE 


ted 
ed 


Mol. 
Mol. 
Vo). 


Vo). 
Mo}. 


ENTHALPY CHANGE Predic 
(J/moi) Measur 
FLOW RATE Feii Prec. 
(¥/min) or -Meas. 
(mmo l/min) Caic. 
coz Meas. 
Caitc. 
Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (K) 
Heat Exch. Corr. (Jd/min) 


Ambient Air Temp. (K) 


SAS sare) = 
-44483.96 - 


©.67S06 


©.15574 
§.932730 


294.66 
296.38 
234.70 
-0.0084 
296.48 


2 


355.86 
297.46 
3.510 
3.453 
©.083 


28.975 
8.663 
©.3335 
2868.70 
4304.09 


°o.68621 
1.08456 
0.02641 


0.18887 
6.03501 


294.66 
286.89 
294.70 
°0.004 
296.48 


3 


355.31 
297.46 
3.513 
3.456 
0.033 


25.610 
8.907 
0.3478 
-2843.76 
-4360.33 


0.67167 
1.02987 
©.02608 


©.18433 
5.87334 


294.66 
296.99 
294.70 
-0.004 
296.48 
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DATA FROM FILE CO02.3 


INTERVALS MARKED WITH ** WERE MOST STABLE IN RUN 


RUN 1 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 1 2 
TEMPERATURE Lower Iniet 383.88 386.07 
(K) Gutiet 298.59 298.74 
PRESSURE Iniet 3.547 3.564 
(MPa) Outiet 3.378 3.313 
Atmospheric ©.093 0.093 
HEAT EXCHANGE (J/min) 74.6098 78.265 
HEAT LEAK (J/min) 19.331 18.287 
RATIO: LEAK OVER EXCHANGE 0.2891 0.2337 


ENTHALPY CHANGE Predicted 
(J/mol) Measured 


-6011.41 +8096 43 
-5414.27 -§345.43 


FLOW RATE Fe 99 Pred. Mol. 
(l/min) or Meas. Mol. 2.99934 3.14631 
(amo l/min) Calc. Voi. 0.07624 0.073988 

€02 Meas. Voi. 6.36245 ©.38511 
Calc. Moi. 13.78005 14.64145 


2.22220 2.41115 


Freon Set‘n Temp. (K) 294.56 294.50 
Freon Meter Temp. (K) 296.26 296.19 
Freon Return Temp. (K) 298.58 294,54 
Heat Exch. Corr. (J/min) “0.013 =0.013 
Ambient Air Temp. (K) 235.09 295.09 
RUN 2 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL i 2 
TEMPERATURE Lower Iniet 386.86 386.67 
(K) Outiet 298.85 298.85 
PRESSURE Iniet 3. 530 3.526 
(MPe) Cuties 32208 3.2:78 
Atmospheric ©.03s3 6.033 
HEAT EXCHANGE (J/min) 79.360 783.360 
HEAT LEAK (J/min). 19.526 18.640- 
RATIO: LEAK OVER EXCHANGE 0.2460 0.2475 


°4115.77 °4107.88% 
“5458.85 ~5458.85 


ENTHALPY CHANGE Predicted 
(J/mo1) Measured 
FLOW RATE Fe11 Pred. Mol. 2.40546 2.80084 
(t/min) or Meas. Mel. 3.19040 3.19040 
(mmo i /min) Calc. Vo. ©.08108 ©.08108 
€e2 Meas. Voi. o,.38228 ©.38228 

Caic. Mol. 14,.83783 14.537983 


Freon Sat’n Temp. (K) 294.51 294.81 
Freon Meter Temp. (K) 296.13 296.17 
Freon Return Temp. (K) 294.62 294.62 
Heat Exch. Corr. (J/min) -0.041 -0.0861 
Amoient Air Temp. (K) 295.08 298.09 
RUN 3 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 1 2 
TEMPERATURE Lower Iniet 393.35 392.79 
(K) Outiet 300.26 300.22 
PRESSURE Intet 3.862 3.562 
(MPa) Outiet 3.323 3.323 
Atmospheric 0.093 ©.093 
HEAT EXCHANGE (J/min) 57 769 57.769 
HEAT LEAK (J/min) -3.282 Cy ar fas th} 


RATIO: LEAK OVER EXCHANGE 
ENTHALPY CHANGE Predicted 
(J/mol) Measured -4121.45 -4137.41 
FLOW RATE Feit Pred. Mo}. 2.48522 2.43255 
(l/min) or Meas. Mol. 2.32325 2.32325 
({mmol/min) Calc. Vor. 0.05868 0.05868 
co2 Meas. Voi. 0.36670 ©.36529 
Gaic. Mol. 14.01656 13.86247 


Freon Sat‘n Temp. (K) 294.70 294.70 
Freon Meter Temp. (K) 296.36 2936.35 
Freon Return Temp. (K) 294.78 294.79 
Heat Exch. Corr. (J/min) -0.028 -0.028 
Ambient Air Temp. (K) 295.54 295.54 


3 


387.81 
298.91 
3.557 
3.288 
0.093 


80.824 
18.204 
0.2282 
°4154.73 
-S362.53 


-S1i738 
- 24919 
08259 


oun 


6.39644 
15.07201 


294.50 
296.15 
294.54 
-0.014 
235.09 


3 


392.55 
300.19 
3.562 
3.323 
©.093 


$7 .403 
-2.947 
-0.0513 
43225 ai) 
=O atid 


- 42707 
. 30855 
.05831 


-365293 
-96258 


wo ONN 


294.70 
296.32 
294.79 
-0.028 
295.54 
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387.74 
298.97 
3.537 
3.270 
0.093 


81.190 
18.740 
o.2308 
143.48 
386.81 


-$1083 
. 26388 
-08297 


- 39643 
07198 


284.50 
296.13 
294.54 
°0.0148 
295.09 


zn. 
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386.88 
288.86 
3.526 
3.280 
0.083 


78.895 
18.678 
©.2364 
118.46 
383.77 


- 482483 
-t7871 
-O8071 


-38Sii 
. 64586 


234.51 
296.39 
294.62 
-0.04)1 
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RUN 4 TEST GAS - C02 

CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 365.79 
(kK) Outiet 2937.98 
PRESSURE Iniet 3.568 
(MPa) Ourjiet 3.461 
Atmospheric ©.093 
HEAT EXCHANGE (J/min) 35.111 
HEAT LEAK (J/min) 8.804 
RATIO: LEAK OVER EXCHANGE 0.2508 
ENTHALPY CHANGE Predicted °-3287.60 
(J/mo1) Measured °4387.88 
FLOW RATE Fea Pred. Moi. 1.08793 
Ci/min) or ,Meas. Mol. 1.41200 
(ramol/min) Caine. Voi 0.03568 
C02 Meas. Voi. 0.20984 
Calc. Mol. 8.00186 
Freon Sat’n Temp. (K) 294.68 
Freon Meter Temp. (K) 286.55 
Freon Return Temp. (K) 294.70 
Heat Exch. Corr. (J/min) “0.003 
Ambient Air Temp. (K) 2938.65 

RUN § TEST GAS - CA2 

CONDITIONS FOR TEST INTERVAL i 
TEMPERATURE Lewer Iniet 407.12 
(K) Outiet 302). 10 
PRESSURE Iniet 3.548 
(MPa) Outie? eee 
Atmospheric 0.033 
HEAT EXCHANGE (J/min) 76.410 
HEAT LEAK (J/min) °-14.967 
RATIO: LEAK OVER EXCHANGE °0.1959 
ENTHALPY CHANGE Predictead -47798 .36 
(J/mel) Measured -3896.55 
FLOW RATE Beit Pred. Mol. 3.674706 
(l/min) oF Meas. Moi. 3.07282 
(mmol /min) Caic. Voi. ©.07766 
ce2 Meas. Voi. 0.50121 
Caic. Moi. 19.11963 
Freen Sat’n Temp. (K) 294.68 
Freon Meter Yemp. (K) 237.05 
Freon Return Temp. (K) 294.72 
Heat Exch. Corr. (J/min) 0.015 
Ambient Air Temp. (K) 295.33 


SES 
-8 


1 
1 
° 
° 
7 


2 


366.11 
298.02 
3.558 
3.448 
0.0393 


35.477 

$.232 
©.2602 
296.02 
455.39 


-O8545 
- 42670 
.O3605 
» 20852 
-96263 


294.68 
296.61 
294.70 
-0.003 
295.65 


2 


407.13 
3602.18 


-67168 
-13161 
.-O7915 


11803 


294.68 
2937.06 
2934.72 
°0.015 
295.93 


-4 


419.46 


07421 
. 44180 
-03643 


.21238 
- 11007 


294.68 
296.65 
294.70 
°0.003 
295.65 
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DATA FROM FILE CO 


2.4 


INTERVALS MARKED WITH ** WERE MOST STABLE IN RUN 


RUN 1 TEST GA 
CONDITIONS FOR TE 


TEMPERATURE 
(K) 

PRESSURE 
(MPa) 


HEAT EXCHANGE (J/ 

HEAT LEAK (J/min) 

RATIO: LEAK OVER 

ENTHALPY CHANGE 
(J/mol) 


FLOW RATE Fe-i1 
(i/min) ofr 
(mmol /min) 


co2 


Freon Sat’n Temp. 
Freon Meter Temp. 


S$ <2 002 =2 
ST INTERVAL 1 2 
Lower Iniet 374.50 375.28 
Outiet 298.16 298.17 
Inlet Spas §.119 
Outiet 5.033 5.020 
Atmospheric 0.033 ©o.093 
min) 61.086 61.811 
12,940 13.847 
EXCHANGE On 2109) 0.2240 
Predicted -4398.54 -4425.86 
Measured °-5580.84 -5703.59 
Pred. Mo). 12 93S3'2 1.92824 
‘Meas. Mol. 2.48553 2.48492 
Calc. Vol. 0.06239 0.06314 
Meas. VYo!. 0.28883 ©.28606 
Caic. Moi. 160.84459 10.83725 
(K) 2948.52 284.82 
(K) 297.79 29 72,33 


Freon Return Temp. (Kk) 234.56 294.56 
Heat Exch. Corr (J/min) °-0.013 -0,013 
Ambient Air Temp. (K) 296.21 A hee Sh 
RUN 2 TEST GAS - C02 

CONDITIONS FOR TEST INTERVAL j 2 


TEMPERATURE Lower Iniet 32 Sha7 7 385.84 
(K) Outiet 2998..19 299.20 
PRESSURE Iniet 5.095 §.082 
(MPa) Gutiet 4.945 4.943 
Atmospheric ©.033 0.083 
HEAT EXCHANGE (J/min) 78.373 73.273 
HEAT LEAK (J/min) NOS 7241 196.869 
RATIO: LEAK OVER EXCHANGE 0.1368 0.1387 
ENTHALPY CHANGE Predicted -473930.65 -8&4792.23 
(J/mol) Measured -5548.81 -5563.87 
FLOW RATE Feit Pred. Mol. 2e71973 2.771376 
(i/min) or Meas. Mol. 3.15072, 3.15072 
(mmoi/min) Caic. Vol. 0.08007 0.08007 
C02 Meas. Voi. ©.37378 0.37284 
Caic. Mel. 14.12179 14.08612 
Freon Sat’n Temp. (K) 294.81 294.51 
Freon Meter Temp. (K) 298.38 298.40 
Freon Return Temp. (K) 294.52 294.52 
Heat Exch. Corr. (J/min) 9.003 0.003 
Ambient Air Temp. (K) 237.08 237.04 
RUN 3 TEST GAS = C02 x2 
CONDITIONS FOR TEST INTERVAL 1 2 
TEMPERATURE Lower Iniet 404.72 4204.43 
(K) Outiet 301.35 301.32 
PRESSURE Iniet S.2106 Saiki: 2 
(MPa) Outlet 4.863 4.858 
Atmospheric 0.093 0.083 
HEAT EXCHANGE (J/min) 106.789 106.789 
HEAT LEAK (J/min) 1.837 1726 
RATIO: LEAK OVER EXCHANGE 0.0135 ©.0162 
ENTHALPY CHANGE Predicted -5473.61 -5458.58 
(J/mo1) Measured -5548.26 -5548.24 
FLOW RATE F-11 Pred. Mo}. 4.23496 4.22335 
(1/min) or Meas. Mo). 4.28272 4.29272 
(mmol/min) Cate; Vote 0.10924 ©.10924 
co2 Meas. Vo!. ©.50970 ©.50970 
Caic. Mol: 19.24728 18.28734 
Freom Sat’n Temp. (K) 294.47 294.47 
Freon Meter Temp. (K) 237.58 237.87 
Freon Return Temp. (K) 294.56 294.56 
Heat Exch. Corr. (J/min) -0.045 -0.045 
Ambient Air Temp. (K) 296.76 296.76 


-& 
cit} 


776.61 
s939.78 


-69634 
.- 16101 
-08033 


- 37166 
-04149 


294.51 
298.42 
294.52 
°0.003 
297.04 


as 
& 


374.32 
288.19 
$.102 
$.008 
©.0s3 


61.811 
13.902 
0.2249 
°83939 .06 
-S6€7S8.51 


4.92602 
2.48832 
0.06314 


©.28741 
10.89086 


2948.52 
298.01 
284.56 
=0. 013. 
286.21 


2.69846 
3.17571 
0.08071 


©.37095 
14.01486 


294.51 
288.42 
294.52 
°-©.003 
297.04 
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67.370 
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-11812.13°11899.31 
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DATA FROM FILE CO2.5L 
INTERVALS MARKED WITH ** WERE MOST STABLE IN RUN 
RUN 1 TEST GAS - COZ es ed 
CONDITIONS FOR TEST INTERVAL 1 2 3 4 
TEMPERATURE Lower Iniet 386.26 395.15 393.72 3982.45 
(K) Outiet 302.73 302.59 302.36 302.16 
PRESSURE Iniet 9.735 8.728 S733 S734 
(MPa) Outlet 8.673 9.673 8.683 9.683 
Atmospheric 0.0393 0.093 0.0383 0.0Ss3 
HEAT EXCHANGE (J/min) 180.276 177.199 171.046 169.508 
HEAT LEAK (J/min) 11.185 18.817 11.859 14.085 
RATIO: LEAK OVER EXCHANGE 0.0620 ©.0876 ©.0693 0.0831 
ENTHALPY CHANGE Predicted -11658.53-11619.39-11573.61-11532.46 
(JS/mol) Measured °124828.73°12734.49-12435.79-12577.54 
FLOW RATE Feit Pred. Mo? 6.796860 6.49904 6.39873 6.24743 
(i/min) or Meas. Mol 7.24641 7.12274 6.87541 6.81358 
(mmo l/min) Calc. Vol. 0.18455 0.18140 0.17510 0.17383 
co2 Meas. Vol 0.38369 0.36812 0.36387 0.38683 
Calc. Mol 14.50360 13.91481 13.75435 13.47703 
Freon Sat‘nmn Temp. (K) 294.45 294.45 294.45 294.45 
Freon Meter Temp. (K) 298.12 298.09 298.13 298.17 
Freon Return Temp. (Kk) 294.36 294.36 294.36 2948.36 
Heat Exch. Corr. (Jd/min) 0.076 0.074 6.072 0.071 
Ambient Air Temp. (K) 236.21 296.21 236.21 296.21 
RUN 2 TEST GAS - C02 zs 
CONCITIONS FOR TEST INTERVAL 1 2 3 4 
TEMPERATURE Lower Iniet 448.04 448.64 449.64 450.58 
(K) Outiet 312.89 Si 2 .. 37, 33) 4S 313.29 
PRESSURE Intet 9.728 $.742 9.745 8.732 
(MPa) Outlet 8.528 $.545 9.545 SS 23 
Atmospheric 6.0383 0.033 0.0833 0.033 
HEAT EXCHANGE (J/min) 387.912 390.988 397.141 392.526 
HEAT LEAK (J/min) °-12.8855 12), 917.5 - 13). 20s coy Wr Ley oy fe 
RATIG: LEAK OVER EXCHANGE °0.0331 -0.0332 -90.0345 °-6.08583 
ENTHALPY CHANGE Predicted ©“12229'.69- 42247. 19:9 12237. 75<12221.23 
(J/mo1) Measured -131837.41°118S53.80-11829.13-11691.65 
FLOW RATE Le ft Pred. Moi. 16.116058 16.23909 16.51631 16.49405 
(i/min) or Meas. Moi. 18.59385 15.71751 15.96485 15.77935 
(mmo l/min) Caic. Vo). 0.39661 0.39975 0.40604 6.40133 
C62 Meas. Vol. 0.86649 0.87216 0.88773 6.88773 
Calc. Mol. 32.76996 32.98415 33.57307 33.57318 
Freon Sat‘n Temp. (K) 294.49 294.49 294.48 294.49 
Freon Meter Temp. (K) 298.80 298.83 2938.90 238.37 
Freon Return Temp. (K) 294.42 294.42 294.42 294.42 
Heat Exch. Corr. (J/min) 0.126 OF B27, 0.129 0.128 
Ambient Air Temp. (K) 296.48 296.48 296.48 296.48 
RUN 3 TEST GAS - CO2 zB ae 
CONOXTYONS FOR TEST INTERVAL 1 2 3 4 
TEMPERATURE Lower Iniet 458.33 458.90 458.62 458.40 
(K ) Outiet 314.95 314.84 314.73 314.63 
PRESSURE Inlet 9.763 9.763 9.759 9.759 
(MPa) Outiet 9.405 9.415 9,408 $.403 
Atmospheric °.03s3 0.033 °©.093 0.0s3 
HEAT EXCHANGE (J/min) 430.955 429.417 426.341 426.341 - 
HEAT LEAK (J/min) °22.440 -23.388 "22.784 °-23.123 
RATIO: LEAK OVER EXCHANGE -0.0821 “0.0545 -0.0538 “0.0542 
ENTHALPY CHANGE Predicted °11964.58-11990.62-12008 .50-12025.24 
(J/mo1) Measured -11372.40-11371.28-11399.32-11406.59 
FLOW RATE Fe-i1 Pred. Mol. 18.22729 18.20357 18.05560 18.06926 
(l/min) or Meas. Mo}. 17.32516 17.26334 17.13968 17.13968 
(mmo 1/min) Saic. Vet. 0.44023 0.43866 6.43552 0.43552 
ca2 Meas. Vol. 1.00100 0.99752 0.98794 0.98731 
Calc. Mol. 37.89482 37.76331 37.40059 37.37674 
Freon Sat’n Temp. (K) 294.52 294.82 294.52 294.52 
Freon Meter Temp. (Kk) 298.83 298.85 298.86 298.82 
Freon Return Temp. (K) 294.38 294.37 294.36 294.36 
Heat Exch. Corr. (J/min) 0.281 0.301 0.319 0.313 
Ambient Air Temp. (K) 296.48 296.48 296.48 296.48 


.- 16389 
- 78175 
-17195 


.- 35239 
.- 32031 


294.45 
298.17 
294.36 

0.070 
296.21 


94. 
14. 
0.0364 
223.61 
733.90 


-41832 
84118 
40290 


88348 
81252 


294.49 
298.93 
294.43 

0.7110 
296.48 
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-O1132 
751758 
-17195 


- 34405 
-00508 


294.85 
298.19 
294.36 

0.070 
296.21 
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RUN 4a TEST GAS - CO2 bed 


CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 438.43 
(K) Outiet 310.45 
PRESSURE Inlet §$.749 
(MPa) Outiet 9.565 
Atmospheric 0.033 
HEAT EXCHANGE (J/min) 341.135 
HEAT LEAK (J/min) -18.756 
RATIG: LEAK OVER EXCHANGE -0.0550 
ENTHALPY CHANGE Predicted °12582.74 
(J/moi) Measured -11926.98 
FLOW RATE Fim) t Pred. Mol. 14.468715 
(i/min) or Meas. Moi. 13.71413 
(mmo i/min) Calc. Vol. 0.34852 
€o2 Meas. Yol. 0.75606 
Calc. Mol. 28.601S7 
Freon Sat‘’n Temp. (K) 294.51 
Freon Meter Temp. (K) 298.64 
Freon Return Temp. (K) 294.24 
Heat Exch. Corr. (J/min) 0.445 
Ambient Air Temp. (K) 296.65 

RUN Ss TEST GAS - COZ 

CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 431.25 
(K) OQutiet 309.34 
PRESSURE Iniet 9.741 
(MPa) Outiet 9.567 
Atmospheric 0.092 
HEAT EXCHANGE (J/min) 311.120 
HEAT LEAK (J/min) -8.4&86 
RATIO: LEAK OVER EXCHANGE -90.0273 
ENTRALPY CHANGE Predicted -12468.94 
(J/mol) Measured -12134.93 
FLOW RATE Foi Pred. Mol. 12.84321 
(i/min) or Meas. Mol. 12.50220 
(mmol/min) Eaic. Vol. 0.31993 
€02 Meas. Voi. 0.88243 
€aic. Mol. 25 .63840 
Freon Sat’n Temp. (K) 294.30 
Freon Meter Temp. (K) 2993.09 
Freon Return Temp. (K) 294.32 
Heat Exch. Corr. (J/min) °-0.027 
Ambient Air Temp. (K) 296.48 

RUN 6 TEST GAS - CO2 

CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 433.32 
(K) Outiet 3098.80 
PRESSURE Inlet 9.755 
(MPa) Outiet 9.567 
Atmospheric 0.092 
HEAT EXCHANGE (J/min) 303.437 
HEAT LEAK (J/min) -26.512 
RATIO: LEAK OVER EXCHANGE -0.0874 
ENTHALPY CHANGE Predicted -12495.20 
(J/mol) Measured -11491.16 
FLOW RATE E= 44 Pred. Mol. 13.25838 
(i/min) or Meas. Mol. 12.19304 
(mmol/min) Calc. Vol. 0.31218 
co2 Meas. Vol. 0.70367 
Calc. Mol. 26 .40607 
Freon Sat‘’n Temp. (K) 294.29 
Freon Meter Temp. (K) 239.13 
Freon Return Temp. (K) 294.38 
Heat Excn. Corr. (J/min) -0.137 


Ambient Air Temp. (K) 296.65 
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$33.42 


- 40695 
- 73886 
- 34915 
- 75323 
- 49481 


294.51 
298.62 
294.22 

0.478 
236.65 


2 


431.66 
303.41 
8.745 
9.570 
0.092 


14.198 
-$.629 
0.0306 
475.29 
104.32 


-01280 


-63093 
-35747 


294.30 
299.06 
234.32 
-0.028 
296.48 


2 


433.10 
309.77 
3.752 
9.562 
0.092 


303.437 


“12 
a He | 


24.953 
0.0822 
486.34 
§37.55 
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28 
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-12 


13 
12 
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° 
25 
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=31.2' 
=20t 


13 
12 
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°. 
26. 


$88.46 


- 35278 
-67702 
- 34758 


7$040 
- 38774 


294.51 
2938.60 
294.22 

0.476 
296.65 


18.737 
-8.336 
©.0264 
484.69 
163.55 


-O2266 
- 68770 
- 32468 


.-69093 
-957858 


284.30 
288.99 
284.32 
-0.028 
296.48 


309.71 
9.752 
9.567 
0,092 


00.389 
24.547 
0.0817 
473.98 
$37.10 


-OS573 
06937 
- 30902 


693376 
03416 


294.28 
299.14 
294.38 
°-0.136 
296.65 
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432.76 
309.63 
9.752 
9.8567 
9.092 


18.814 
°-9.510 
0.0298 
494.89 
132.96 


.- 18381 
81137 
32784 


69942 
27666 


298,30 
298.97 
294.32 
-0.028 
296.48 


94.204 


0.1049 
486.36 
300.67 


-06241 
. 82204 


0.30269 


.69376 
03418 


294.29 
2939.18 
294.38 
°-0.133 
296.65 
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833.07 
308.69 
Sir si2 
3.870 
©.092 


18.814 
-$.7290 
0.90305 
$02.83 
132.92 


20194 
81137 
32784 


69943 
27675 


294.30 
293.02 
294.32 
-0.028 
296.48 
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RUN 7 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outlet 
PRESSURE Iniet 
(MPa) Outiet 


Atmospheric 


HEAT EXCHANGE (J/min) 

HEAT LEAK (J/min) 

RATIO: LEAK OVER EXCHANGE 

ENTHALPY CHANGE Predicted 
(J/mol) Measured 


FLOW RATE Ley | | Pred. Mol. 
(l/min) or Meas. Mol. 
(mmol/min) Calc. Voi. 


co2 Meas. Voi. 
Caic. Mol. 


Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (K) 
Heat Exch. Corr. (J/min) 
Ambient Air Temp. (K) 


RUN & TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outiet 
PRESSURE Intet 
(MPa) Outiet 


1 2 
4S1.493 451.78 
313.63 313.68 

$.721 $.724 

3.407 3.405 

0.092 0.092 

400.390 401.928 
-11.643 - 14.300 
°0.0291 -0.0356 
-12033.04-12028.43 
-11693.00-11615.18 
16.5S637 16.72493 
16.08852 16.15034 
0.41209 0.41367 
0.91322 0.92287 
34.24182 34.60371 
294.28 294.28 
293.34 299.36 
294.42 294.42 
°0.281 -0.282 
296.76 296.76 

1 2 
430.98 4830.77 
303.25 309.22 

9.732 9.732 

9.418 $.418 

o.03g33 0.093 


Atmospheric 


HEAT EXCHANGE (J/min) 


HEAT LEAK (J/min) 


RATIO: LEAK OVER EXCHANGE 


ENTHALPY CHANGE 
(J/mol) 


Predicted 
Measured 


°12380.89-12374.76 
-11882.20-11819.34 


FLOW RATE Liteon i Pred. Mol. 
(i/min) or Meas. Mol. 
(mmol /min) Calc. Yo. 


c€02 Meas. Vo). 
Calc. Mol. 


Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (K) 


Heat Exch. Corr. (J/min) 
Ambient Air Temp. (K) 


RUN 9 TEST GAS - C02 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outlet 
PRESSURE Iniet 
(MPa) Outiet 


Atmospheric 


HEAT EXCHANGE (J/min) 
HEAT LEAK (J/min) 
RATIO: LEAK OVER EXCHANGE 


ENTHALPY CHANGE 


Predicted 


(J/mol) Measured 
FLOW RATE F-11 Pred. Moi. 

(i/min) or Meas. Mol. 

(mmol/min) Caics Vel. 


co2 Meas. Vol. 
Calc. Mol. 


Freon Sat’n Temp. (K) 


Freon Meter Temp. (K) 
Freon Return Temp. (K) 
Heat Exch. Corr. (J/min) 


12.18934 12.18336 
11.69838 11.63654 
0.293757 0.23600 
0.64862 6.648562 
24.49129 24.49139 
294.48 294.48 
29757 F 297.80 
294.54 294.54 
-0.078 -0.077 
295.54 295.84 

=e 

1 2 
350.51 350.397 
236.65 296.63 
3.564 3.575 
3.494 3.504 
°.0s3 ©.0s3 
$4.133 $2.903 
41.652 40.013 
©.7694 0.7564 
-2671.07 -2696.49 
-11585.02-11067.14 
0.80173 0.51818 
rye ae A 2.12664 
0.08535 0.05410 
0.12318 0.12601 
4.67268 4.78015 
294.48 294.48 
296.80 296.86 
294.53 294.54 
“0.012 -0.014 
295.54 295.54 


Ambient Air Temp. %“K) 


401.928 
°11.587 
-0.0288 
-12023.70 
-11692.62 


16.61592 
16.18034 
0.41367 


0.91676 
34.374854 


294.28 
293.34 
294.42 
°0.282 
296.76 


287.3934 
-18.S86 
°0.0646 
-12382.92 
°11603.88 


12.32188 
11.574&70 
0.29443 


0.65412 
24.81358 


294.468 
297.71 
284.54 
-0.077 
295.54 
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DATA FROM FILE €02.6 


INTERVALS MARKED WITH *2 WERE 


RUN 1 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Out iet 
PRESSURE Iniet 
(MPa) Outies 


Atmospheric 


HEAT EXCHANGE (J/min) 

HEAT LEAK (J/min) 

RATIO: LEAK OVER EXCHANGE 

ENTHALPY CHANGE Predicted 
(J/mol) Measured 


FLOW RATE Fient't Pred. Moi. 
(1/min) or ‘Meas. Mol. 
(ramo l/min) ‘Catic., Vol. 

C02 Meas. Voi. 
Caic. Mol. 


Freon Sat’n Temp. (K) 
Freon Meter Yemp. (K) 
Freon Return Temp. (K) 
Heat Exch. Corr. (J/min) 
Ambient Air Temp. (K) 


RUN 2 TEST GAS - C€O2 
CONDITIONS FOR TEST INTERVAL 


TEMPERATURE Lower Iniet 
(K) Outiet 
PRESSURE Iniet 
(MPa) Outiet 
Atmospheric 


HEAT EXCRANGE (J/min) 

HEAT LEAK (J/min) 

RATIO: LEAK OVER EXCHANGE 

ENTHALPY CHANGE Predicted 
(J/mol) Measured 


FLOW RATE Ft Pred. Mo}. 
(i/min) or Meas. Mol. 
(mmoi/min) Caitic. Vol; 

coz Meas. Vol. 
Caic. Mol. 


Freon Sat’‘n Temp. (K) 
Freon Meter Temp. (K) 
Freon Return Temp. (K) 
Heat Exch. Corr. (J/mi 
Ambient Air Temp. (K) 


RUN 3 TEST GAS - CO2 
CONDITIONS FOR TEST INTERVAL 
TEMPERATURE Lower Iniet 
(K) Outiet 
PRESSURE Iniet 
(MPa) Outiet 


Atmospheric 


HEAT EXCHANGE (J/min) 

HEAT LEAK (J/min) 

RATIO: LEAK OVER EXCHANGE 

ENTHALPY ChANGE Predicted 
(J/mot) Measured 


FLOW RATE F-11 Pred. Mol. 
(i/min) or Meas. Mol. 
(ramo l/min) Catic. Vol. 


co2 Meas. Vol. 
Calc. Mo}. 


Freon Sat’n Temp. (K) 
Freon Meter Temp. (K) 
Freom Return Temp. (K) 
Heat Exch. Corr. (J/min) 


Ambient Air Temp. (K) 


MOST STABLE IN RUN 


=2 


1 2 


391.26 391.53 
299.31 289.35 


3.536 3.836 
3.298 3.295 
©.095 0.098 


73.401 73.884 
12.822 12.802 
0.1747 o.1688 
°4302.57 -4311.54 
-$213.28 -5185.82 


2.43800 2.46326 
2.95524 2.96258 
0.07328 0.07348 
©.36104 6.36387 
14.07872 14.19019 


295.24 235.24 
237.42 237.42 
295.28 298.29 
=O. O1 7 =O 007 
238.37 295.37 


1 2 


390.57 330.47 
298.406 2939.38 


3.8532 3.515 
3.268 3.260 
°o.o0ss ©.0ss 


2.48827 2.40981 
3.014803 2.98464 
0.07483 9.07410 
0.36698 ©.36151 
14.30907 14.08569 


288.21 288.21 
297.14 236.87 
298.21 295.21 
©.000 ©.000 
295.09 295.09 


i 2 
380.21 379.84 
298.56 298.49 

3.635 3.628 

3.475 3.467 

0.0gss ©o.095 


°o.2855893 0.28411 
11.08180 11.02335 


295.10 295.10 
297.22 297.18 
298.21 295.21 
-0.032 -0.032 
295.37 2358.37 
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DATA FROM FILE MIX. i 


INTERVALS MARKED WITH *2 WERE MOST STABLE IW RUN 


RUN 1 TEST GAS - TEST 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 384.73 
(K) Outiet 301.89 
PRESSURE Iniet 7.005 
(MPa) Outiet 6.663 
Atmospheric ©.093 
HEAT EXCHANGE (J/min) 112.637 
MEAS. ENTH. CHANGE (J/mol1) -34938.39 
FLOW RATE Feit Meas. Mol. 4.52788 
(mmoil/min) Caic. Voi. o.11s19 
or (1i/min) tee e coerce 
Tes? Meas. Voi. ©.8as58i6 
"Calc. Mol. 32.18764 
Freon Sat‘n Temp. (Kk) 294.48 
Frean Meter Temp. (K) 237.83 
Freon Return Temp. (K) 294.56 
Heat Exch. Corr. (J/min) -0,.043 
Ambient Air Temp. (K) 297.08 
RUN 2 TEST GAS - TEST 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 372.53 
(K) Outie? 300.47 
PRESSURE Tniet 7.067 
(MPa) Outiet 6.828 
Atmospheric 0.0693 
HEAT EXCHANGE (J3/min) 8&.180 
MEAS. ENTH. CHANGE (J/moi1) °3272.08 
FLOW RATE F-i1 Meas. Mol. 3.54314 
(mmol /min) Ceaic. ¥or. ©.09014 
or (l/min) SORE CNIS ae 
TEST Meas. Voi. ©.71500 
Catc. Mo}. 26.93716 
Freon Sat‘n Temp. (K) 294.48 
Freon Meter Temp. (K) 297.03 
Freon Return Temp. (K) 2946.57 
Heat Exch. Corr. (J3/min) “60.038 
Ambient Air Temp. (K) 296.76 
RUN 3 TEST (GAS =) (TEST ze 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 419.13 
(K) Outiet 303.07 
PRESSURE Iniet 7.029 
(MPa) Outiet 6.780 
Atmospheric ©.093 
HEAT EXCHANGE (J/min) 137.134 
MEAS. ENTH. CHANGE (J/mo1) -S080. 8&2 
FLOW RATE F-11 Meas. Mol. 5.51261 
(mmoi/min) Cates Vers 0.14024 
or (1/min) ER CIO i 
TEST Meas. Vo). O.71642 
Calc. Mo}. 26 .99086 
Freon Saet’nm Temp. (K) 294.48 
Freon Meter Temp. (K) 289320 
Freon Return Temp. (K) 294.56 
Heat Exch. Corr. (J/min) -0.052 


Ambient Air Temp. (K) 296.76 


2 


384.87 
30; .$3 
7,012 
6.665 
0.083 


113.738 
-3833.45 


4.57197 
©.17631 
©.8S517 
32.18774 


294.48 
237.94 
294.56 
-0.084 
297.08 
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372.78 
300.487 
7.6076 
6.833 
©.083 


87.044 
°3231.37 


1 
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°. 


fe} 
32 
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1 
“4 
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384.39 
301.95 
7.015 
6.670 
©.0393 


13.734 
$845.19 


-§7187 
11631 


- 85234 
-08122 


294.48 
287.94 
294.56 
°@.088 
2937.08 
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373.85 
300.53 
7.070 
6.828 
o.083 


gs .968 
268.75 


-61663 
-09201 


- 73087 
-§23380 


294.48 
297.38 
294.84 
°0.026 
296.76 
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418.17 
303.098 
7.032 
6.778 
0.093 


35.671 
910.19 


- 48382 
.13875 


- 73340 
-63052 


294.48 
227 a 
294.56 
-0.082 
296.76 


ze 
4 


385.09 
301.95 
7.015 
6.673 
0.093 


113.734 
-3845.20 


4.87197 
0.173631 


©.88233 
32.08112 


294.48 
297.98 
294.56 
-0.04844 
287.04 
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378.1% 
300.66 
7.078 
6.823 
©.0$s3 


$1.431 
-3283.71 


3.67542 
©.08350 


0.73807 
27.84387 


294.48 
297.48 
294.54 
20.626 
296.76 
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8138.97 
303.23 
7.039 
6.778 
©.093 


138:231 
-49485 .83 


§.S5671 
0.14137 


©.74190 
27.35059 


294 88 
297.07 
2946.56 
-0.0583 
296.76 
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385.17 
301.97 
7.018 
6.673 
6.093 


14.465 
567.98 


-6€0137 
- 11706 


- 85234 
08122 


294.48 
297.90 
294.56 
-0,0448 
297.04 
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RUN &@ TEST GAS - TEST 

CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 421.93 
(K) Outiet 303.65 
PRESSURE Iniet 6 S74 
(MPa) Outiet 6.683 
Atmospheric ©.,093 
HEAT EXCHANGE (J/min) 158.7638 
MEAS. ENTH. CHANGE (J/mol) -§243.89 
FLOW RATE F-it Meas. Mol. 6.26218 
(mmol /min) Calc. Vol. ©.15910 
or (1/min) weet ec ceee 
TEST Meas. Vol. 0.78438 
“Calc. Mol. 29.70883 
Freon Sat’n Temp. (K) 294.82 
Freon Meter Temp. (K) 296.47 
Freon Return Temp. (K) 294.62 
Heat Exch. Corr. (J/min) °0.074 
Ambient Air Temp. (K) 298.65 

RUN 5 TEST GAS - TEST 

CONDITIONS FOR TEST INTERVAL i 
TEMPERATURE Lower Iniet 398.73 
(kK) Outiet 302.50 
PRESSURE Iniet 7.001 
(MPa) Outiet 6.713 
Atmospheric 0.0393 
HEAT EXCHANGE (J/min) T2277 TI4. 
MEAS. ENTHK. CHANGE (J/mol) -4114.55 
FLOW RATE F-11 Meas. Mol. 4.314001 
(mmol /min) Caic. Vol. 0.12476 
or (1/min) 7oe CS E> se'ee 
TEST Meas. Voi. 0.78837 
Calc. Mo}. 23 .68350 
Freon Sat‘n Temp. (K) 238.82 
Freon Meter Temp. (K) 296.84 
Freon Return Temp. (K) 298.5856 
Heat Exch. Corr. (J/min) -0.025 
Ambient Air Temp. (K) 295.82 


as 
2 


422.42 
303.71 
6.377 
6.685 
0.093 


147.725 
-4988.39 


5.93884 
©.18089 
©o.78198 
29.61368 


294.52 
296.50 
294.62 
-0.070 
295.65 
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399.65 
302.81 
7.001 
Uo WAS 
o.0s3 


122.134 
°-4114.54 


4.91001 
0.12476 
©.78438 
23 .68356 


294.582 
296.92 
284.586 
0.025 
235.82 


399.67 
302.51 


122.500 
°4481.81 


4.9247t 
6.12814 


0.78184 
28 .57637 


294.52 
296.32 
294.56 
°0.026 
295.82 
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399.65 


123.231 
-41$1.50 


4.985411 
o.12588 


0.78437 
23 .68346 


298.52 
296.91 
294.56 
-0.026 
295.82 
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399.65 
3602.51 
7.001 
6.708 
0.033 


21.769 


°$102.22 
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-89531 
.124398 


78438 
- 68356 


294.82 
286.91 
294.56 
°60.025 
295.82 
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DATA FROM FILE MIX.2L 


INTERVALS MARKED WITH ** WERE MOST STABLE IN RUN 


RUN 1 TEST GAS - TEST xs 
CONDITIONS FOR TEST INTERVAL 1 2 
TEMPERATURE Lower Iniet 417.67 418.03 
(K) Outiet 302.81 302.87 
PRESSURE Inlet 6.981 6.985 
(MPa) Outiet 6.743 6.781 
Atmospheric 0.0933 0.093 
HEAT EXCHANGE (J/min) 130.985 134.060 
MEAS. ENTH. CHANGE (J/mo1) -49335.77 -5038.19 
FLOW RATE Fett Meas. Mol. §.26775 §.39142 
(mmoi/min) Caic. Vol. 0.13305 0.13617 
or (l/min) cower ts ec ecco cceeoce 
TEST Meas. Vol. 0.69505 6.69691 
Calc. Mol. 26.53786 26.60872 
Freon Sat’n Temp. (K) 294.70 294.706 
Freon Meter Temp. (K) 296.32 296.33 
Freon Return Temp. (K) 2948.68 294.68 
Heat Exch. Corr. (J/min) Qo.012 O02 
Ambient Air Temp. (K) 295.09 295.09 
RUN 2 TEST GAS - TEST aE 
CONDITIONS FOR TEST INTERVAL 1 2 
TEMPERATURE Lower Iniet 632.45 432.18 
(K) Outlet 305.66 305.66 
PRESSURE Inlet? 6.388 6.988 
(MPa) Outiet 6.613 6.611 
Atmospheric 0.0s3 0.093 
REAT EXCHANGE (J/min) 174.038 177 -4K:3 
MEAS. ENTH. CHANGE (J/mol) -5160.81 -5201.89 
FLOW RATE Fiesta Meas. Moi. 6.393907 Wi ieee 
(mmo i/min) Calc. Voi. 0.17682 0.17995 
or (1/min) croc sece eeccc cc ececes 
TEST Meas. Vo}. 6.88348 °o.89198 
Caic. Mol. 33.72290 34.08710 
Freon Sat’n Temp. (K) 294.69 294.69 
Freon Meter Temp. (K) 296.83 296.50 
Freon Return Temp. (K) 238.68 294.68 
Heat Exch. Corr. (d/min) ©.0093 0.009 
Ambient Air Temp. (K) 295.08 295.09 
RUN 3 TEST GAS - TEST 
CONDITIONS FOR TEST INTERVAL i 2 
TEMPERATURE Lower Iniet 426.36 426.02 
(K) Outiet 304.63 304.57 
PRESSURE Iniet 7.026 7.026 
(MPa) Outiet 6.716 6.716 
Atmospheric 0.093 0.093 
HEAT EXCHANGE (J/min) 158.885 185.585 
MEAS. ENTH. CHANGE (J/moi1) -§049.27 -5048.27 
FLOW RATE | fo 4) Meas. Moi. 6.25708 6.28708 
(mmol/min) Calc. Vol. 0.15804 0.15804 
ha (MIVA Ya 8 Cris. oP IT Ot 
TEST Meas. Voi. 0.80703 0.80703 
Catc. Mol. 30.81332 30.81332 
Freon Sat’n Temp. (K) 284.70 294.706 
Freon Meter Temp. (K) 286.38 296.31 
Freon Return Temp. (K) 294.66 294.66 
Heat Exch. Corr. (J/min) ©.0293 0.029 


Ambient Air Temp. (K) 295.09 295.09 
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418.62 


137.136 
-5084.50 


5.S1S09 
9.13931 


0.70650 
26 .97136 


294.69 
296.34 
294.68 

0.010 
295.09 
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432.47 
305.73 
6.9388 
6.608 
0.033 


174.038 
-5111.66 


6.99907 
0.17682 


©.8sig98 
34.98720 


» 294.69 
286.50 
294.68 

0.009 
295.09 
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425.59 
304.52 
7,026 
6.718 
0.093 


187.122 
°5183.41 


6.31882 
©.1583960 


0.79853 
30.48904 


294.70 
296.27 
284.67 

0.022 
295.09 
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420.27 
303.18 
7.009 
6.761 
°o.03s3 


140.211 
-4969 .47 


§.63875 
0.14244 


0.73907 
22.21446 


294.69 
296.37 
294.66 

0.9028 
295.09 


Re 
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4832.87 
305.78 
6.988 
6.606 
©.093 


479-003 
-§1858.52 


7.122784 
6.17995 


0.89481 
34.15523 


294.69 
296.50 
294.68 

©.009 
285.09 
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421.73 
303.42 
7.008 
6.748 
0.033 


46.362 
168.42 


. 88608 
-14871 


- 74190 
-31887 


294.69 
296.40 
294.64 

0.036 
2935.09 
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421.92 
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7.023 
6.766 
©.033 


46.362 
188.38 


. 88608 
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RUN 4 TEST GAS - TEST 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 416.32 
(K) Outlet 302.93 
PRESSURE Inlet 7.002 
(MPa) Outiet 6.763 
Atmospheric 0.9093 
HEAT EXCHANGE (J/min) 134.063 
MEAS. ENTH. CHANGE (J/mol) °4937.10 
FLOW RATE Fett Meas. Mol. §.39142 
(mmoi/min) Caice ivoilic ©.13623 
or (i/min) cor cccccces 
TEST Meas. Voi. Ora 27, 
Calc. Mol. 27.15424 
Freon Sat’nm Temp. (K) 294.68 
Freon Meter Temp. (K) 286.72 
Freon Return Temp. (K) 294.64 
Heat Exch. Corr. (J/min) 0.030 
Ambient Air Temp. (K) 298.37 
RUN 5 MGST GaAs = TEST ae 
CONDITIONS FOR TEST INTERVAL 1 
TEMPERATURE Lower Iniet 406.95 
(K >) Gutiet 301.53 
PRESSURE Inlet 7.036 
(MPa) Outiet 6.853 
Atmospheric 0.093 
HEAT EXCHANGE (J/min) 118.615 
MEAS. ENTH. CHANGE (J/mol1) -§038.76 
FLOW RATE F-11 Meas. Mol. 4.64942 
(mmol /min) Calcm Vet. 0.11783 
or (l/min) wecrceccceceres 
TEST Meas. Vol. 0.60238 
Caic. Mol. 22.94§822 
Freon Sat’n Temp. (K) 294.67 
Freon Meter Temp. (K) 297.085 
Freon Return Temp. (K) 294.64 
Heat Exch. Corr. (J/min) ©.019 
Ambient Air Temp. (K) 295.54 


2 


416.87 
302.39 
7.002 
6.763 
0,033 


137.138 
-§030.384 


s.515098 
0.138335 
0.71500 
27.26225 


298.68 
296.72 
294.64 

0.031 
288.37 


2 


407.02 
301.54 
7.036 
6.853 
°©.093 


117.183 
°-5117.93 


4.71126 
6.113909 
0.60095 
22.89072 


294.67 
29 7te 0 
294.64 
@.9020 
295.54 


1 


417.42 
303.05 
6.399 
6.761 
0,093 


37.138 


°§020.42 


vo on 


1 


-51509 
.13935 


- 71641 
- 31609 


294.68 
296.8848 
294.66 

0.018 
295.37 


3 


4o7.11 
301.55 
7.036 
6.883 
0.093 


17.183 


-S111.25 


2 


4 
° 


°. 
2 


-71126 
.11309 


60173 
.92064 


294.67 
287.09 
294.64 

0.020 
298.54 


=e 
4 


417.48 
303.08 
6.995 
6.761 
0.093 


140.213 
°-S$§112.75 


-6387S8 
- 14248 


- 719258 
42421 


“NO ow 


294.68 
296.95 
294.68 

0.0058 
238.37 


ze 
4 


407.16 
301.55 
7.036 
6.8853 
0.033 


115.615 
-§084&.14 


4.64942 
0.11783 


0.60173 
22.92075 


294.67 
287.08 
294.63 

©.028 
288.54 
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Appendix 7 - Error Analysis 


The precision of the results and the calculated heat leaks are 
estimated using the accuracies of the quantities which enter into these 
calculations. A summary of the primary equations involved in this 


estimation follows: 


44) sak) 3 Seared: - d(Ah, was? (A7.1) 
d(Ah og) = 4peag = 4(-A,*m,/m,) (A7.2) 
d(Ah 4) * (-m,/m,)da_ - (A,/m,)dm_ + (A m_/m2)dm, (4.1) 
d Cy pealeces satan is (A7.3) 
dh = (j+£P)dT + (k+£T)dP (A7.4) 


[f£,j, and k are well defined in terms of the] 
[tabular values used in the interpolation ] 


The differential quantities in the equations above are 
approximated by small values representing the uncertainties in the 
corresponding variables. It should be noted that the uncertainties 
are "plus or minus" and use must be made of the case which gives the 
larger (i.e., more conservative) estimate of the uncertainty of the 
quantity of interest. Calculation of these uncertainties for several 


cases is presented below. 
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A. Determination of the uncertainty in the experimental value 
of enthalpy change - d(Ah ) 


meas 


Examine equation (4.1) term by term. 


E =| } ; 
d (P -P)+X 
pt EL oe sat ie 


1/ dd. = 
ae 
=a KdP where K = | ————— | = 15180 J/moleMPa 
sat pt apr 
for ambient conditions 
Since dP = 6.7 x 10° MPa (0205 2torr 4 


day. a-t0iet-al/mol 


This interpolation is at least as accurate as linear 
interpolation of A_ with the reciprocal of saturation 


pressure in the narrow range of ambient conditions. 
An approximate value for h- is 24850 J/mol. 
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2h dm is equal to d(atbV) = bdV or to the maximum correlation 


error, whichever is larger. dv = +0.01 volt. 


For flow meter A, b = 0.2019 giminenevel eke 


Thus, bdV = +0.002 g/min, but the maximum correlation 


error in the range used is about +0.025 g/min. 
Hence, for flow meter A, 


dm. = +0.025 g/min = +1.8 x toseroi iia 


For flow meter B, b = 1.6988 Gere ee eee 


Thus, bdV = +0.017 g/min, but the maximum correlation 


error is about +0.1 g/min in the range used. 
Hence, for flow meter B, 


dm. mit0Q,1°C/ Mineman) wonx 108 ot in 
3/ dm. = d(PAV/2RTt) 
_ Pd(AV) + AVdP - PAVdT/T - PAVdt/t 
zeReTet 
d(AUy) =c20, 00 fe? eh 288 a 10e-m Ay es. 5: x 10D om 
By tay) oo ler ONS = 0.0933 MPa 
dT = +0.25°F = +0.14 K 2296. Louk 
dt = +0.01 min C= "L070 ymin 
R = 8.314 x 107°MI/mol K z = 0.998 


Thus, dm, = +1.4 x Tom mo nin 
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Having the above uncertainty estimates, d(Ah as) can be 


calculated from equation (4.1) for a number of cases: 


a) m. = 10 mmol/min ; m = 1 mmol/min (using N, and flow meter A) 


Z 


= + 
d(Ah as +482 J/mol 


bo. *m =635 mmol /min: m = 5 mmol/min (using N, and flow meter A) 


2 


d(Ah as? = +142 J/mol 


c) m- = 10 mmol/min ; mn = 1 mmol/min (using CO, and meter A) 


£ 2 


= + 
Oe eae ae 4482 J/mol 


d) m. = 35 mmol/min ; mn = 20 mmol/min (using CO, and meter B) 


ie 


d (Ah Je =ot5/5y Ji/mol 
meas 


It is clear that the use of higher process fluid flow rates results 
in lower uncertainties of the measured enthalpy changes. Case d) 
shows that this is negated by the use of flow meter B, made necessary 
by the higher flow rates of the reference fluid. Flow meter B has 


a higher range but also a higher correlation uncertainty. 


B. Determination of the uncertainty in the predicted value 


areas 
Examine equations (A7.3) and (A7.4) term by term after deriving 


of enthalpy change - d(Ah 


equation (A7.4). This equation arises as a result of a two - 


dimensional linear interpolation as follows: 
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Continuing from the previous page, two more convenience terms 
may be defined: 


Let f = —= and g = sist ° 
e e 


The enthalpy at some condition may then be expressed: 


nas (£(P-P ) + g)(T-T ) + a(P-P ) +b 


Upon rearrangement, 


i ae (g-£P )T + (a-fT )P + fPT + (fP T - gT - aP +b) 
9 


With another three convenience terms defined: 
4 = (g-£P) 0k = (a-fT ). Me 1 =2(EPNT) ='egT°’- ae +b) 
the predicted value of enthalpy is written 


hy p SST + ke + fpr + 1 


and the derivative may be taken: 


dh = (j+£P)dT + (k+fT) dP (A7.4) 


Applying this to both inlet and outlet conditions as given in 
equation (A7.3), 


d(‘n ed = cheer - dh, | CAT. 3) 


Thus d(4h ed) is expressed as a function of the measurement 


uncertainties dT and dP. Several cases are now examined. 


Case l. Nitrogen Inlet: qT, = 401.36 K ; Pp. Se 0) Seba 
Outlet: Ty = 301.80 K ; P, = 9.494 MPa 
Now (j+£P) = g + £(P-P) = SSS) +o 


and (k+fT) 


a+f(T-T ) = ae) +a 


ane a pe] yb 
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The calculations for case 1 follow: 
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Inlet 
sy APRS get ae OW ah es Jymole™ | 5 
= aE = -20 Mpa Ae? D = 2/00/)J (m0) 
_ 1428302 4209 Sayed mol a 
c= 0.5 = -12 “Pa? d = 4289 J/mol 
and e = 450 - 400 = 50 K 
Thus, (j+£P) = 93101. J3/molsk ands (kf£T) = -19.8 dimol 
» J in . in : MPa 
Outlet 
HOE (46 8) ry Se MOD a 
gare CECT Fa aam PR 46 are. b 448 J/mol 
Se LS Ghali Le, J/mol : L 
c= a = 30 “pa tal 1171 J/mo1 
and e = 350 = 300 = 50 K 
fe : BY, J/mol 
Thus, (TEP ee = 32.5 J/mol*K and (K+£T) ot 45.4 Mpa 
Hence, d(Ah ed? = £(31.1432.5)dT + (19.8+45,.4) dP 
=e vod l, Oo e2ecr 
Since dT = t0slek. -and *dPS=, +0.01sMPa = (SeeSAppra2) ; 
d (Ah ) = +7.0 J/mol for this case. 
pred 
Case 2. Nitrogen Inlet : Ty ashen oN) AGH Pp. = 6.975 MPa 
Outlet : T, = 3005.7 Kes P. = 6.694 MPa 
Following the method of Case l, 
. eS oh mol 
(j+£P) = Sf 0) Jf mol *Ki2 (k+fT) = -26.5 Pa? 
; ve oe any J/mol 
(pte) = 931,6 J/mol«K: and (K+fT) out = -51.8 Daa 
Thus, d(Ah___,) = +62.6 dT + 78.3 dP = +7.0 J/mol. 
pred 
Case 3. Nitrogen Inket™: qT, =e Je Olas Pr = 3.524 MPa 
Outlet : T, = 298.54 K ; P, = 3.199 MPa 
Similar to the above, 
; = aes aon SMO s. 
(j+£P) so 30;. L.. Senos Kees (k+£T) . = -33.9 pee 3 
; a reas ae J/mol 
GtiP)ns = 31.0 J/mol*K ; and (k+£T) ot = -58.8 SMoES 


Thus, d(Ah | = +61.1 dT + 92.7 dP = +7.0 J/mol. 
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The following cases are for CO, : 


Case 4, inlets: Ty aU 458).40 0K: PL = 9.759 MPa 
Outlets: T, mte3 14) Oakes P, = 9,403 MPa 
Similar to the above, 
; (ktfT), = -177.6 ——— “ized 
(j+EP) = 52,60 J/mol skis. ( Nps = ane 
oe 


CEFR) te = 359 J/mol*K ; and (K+£T) it = -1144 —— 


Thus, d (Ah, 34) = +411.6 dT + 1321.6 dP = +54.8 J/mol. 
Case 5. Inter=: qT, = 374,92 K 3 Pp. = 5.102 MPa 
Ourclet: qT, = 298.19 K ; P, = 5.008 MPa 
As above, 
5 2 oe J/mol 
(j+£P) = 49.6 J/moleK ; (k+fT) . = —301.2 ‘Pa? 
eae 


igetP) = 93.0 J/mol*K ; and (kt+£T) vt = -1043 


Thus, d(Ah, 4) = +142.6 dT + 1344 dP = +27.7 J/mol. 
Case 6. Inlet: T) = 398.56 K. ; Py = 3.499 MPa 
Outlet : T, = 300.20 K ; P, = 3.118 MPa 
As above, 
(G+£P) sO = 45.9 J/mol*K ; (k+fT), ==242,3 Amel ; 
“ nf Pech te mel 
Cohn th yes = 50.4 J/mol*K ; and (k+f£T) out = -570.6 
Thus, d(Ah, 4) = +96.3 dT + 812.9 dP = +17.8 J/mol. 


Summarizing the above, it appears that the uncertainty in the prediction 
of the enthalpy changes is relatively constant for the cases using Nos 
DUE Lor CO,, the uncertainty increases with pressure of the run. 

The appropriateness of the use of (two-dimensional) linear interpolation 
can be judged by comparing (j+fP) with 7 and (k+fT) with @ for a 


number of conditions T and P. Some comparisons are given below. 


Nitrogen : Conditions (j+f£P) fo (k+£T) 0) 
[J/mol °K] [J/mol+MPa] 

Case 1: Inlet Se cb lpes: -19.8 -19.7 

Case 23, Outler (slo SyAMAL -51.8 -51.8 


Case 3: Inlet 300% BO? -33.9 =33.0 
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For carbon dioxide, the I.U.P.A.C. tables (Angus et al., 1973) 


give only the values of ee and comparisons are made below. 


Carbon see Condi tione (j+£P) C_(from tables by 
dioxide faye ek] interpolation) 
Case 4: Inlet 52.6 DZS 
Case 5: Outlet 9370 abi! 
Case 6: inlet 45.9 45.9 


It appears that linear interpolation is generally accurate for 
both No and CO, 
comparing the result of a cubic spline fit with the results of 


linear interpolation. Differences were negligible. 


C. Determination of the uncertainty in the calculated values 


of the heat leak term - dda cak: 
Use equation (A7.1) and examine several cases. 
qq oak = NG ae - d(Ah cas? CAT<1) 
For Nos low flow rate : qq oak F 7 ¢2a482 : +490 J/mol 
aoe ghicterie ee ! : eit] 291 42"=.2150.J/mol 
For CO, > low flow rate : qq) oak 2 ~30) et) 482 F +519 J/mol 
ys Ueehiche 3 4 me 5r 575, eto 30 mol: 


It should be noted that ane different F-1l1 flow meters are 


used for the two co, cases given above. This, unfortunately, 


obscures the trend seen in the Ny cases showing that the 


uncertainty in this calculation decreases as process fluid 


flow rate increases. 


in the regions of use. This was also tested by 
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